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Magnetic clusters, i.e., assemblies of a finite number (between two or three and several hundred) of 
interacting spin centers which are magnetically decoupled from their environment, can be found in 
many materials ranging from inorganic compounds, magnetic molecules, artificial metal structures 
formed on surfaces to metalloproteins. The magnetic excitation spectra in them are determined 
by the nature of the spin centers, the nature of the magnetic interactions, and the particular 
arrangement of the mutual interaction paths between the spin centers. Small clusters of up to 
four magnetic ions are ideal model systems to examine the fundamental magnetic interactions 
which are usually dominated by Heisenberg exchange, but often complemented by anisotropic 
and/or higher-order interactions. In large magnetic clusters which may potentially deal with a 
dozen or more spin centers, the possibility of novel many-body quantum states and quantum 
phenomena are in focus. In this review the necessary theoretical concepts and experimental 
techniques to study the magnetic cluster excitations and the resulting characteristic magnetic 
properties are introduced, followed by examples of small clusters demonstrating the enormous 
amount of detailed physical information which can be retrieved. The current understanding of 
the excitations and their physical interpretation in the molecular nanomagnets which represent 
large magnetic clusters is then presented, with an own section devoted to the subclass of the single- 
molecule magnets which are distinguished by displaying quantum tunneling of the magnetization. 
Finally, some quantum many-body states are summarized which evolve in magnetic insulators 
characterized by built-in or field-induced magnetic clusters. The review concludes addressing 
future perspectives in the field of magnetic cluster excitations. 
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I. INTRODUCTION 

Magnetic clusters are defined as systems which consist 
of a finite number N of interacting spins that are magnet- 
ically isolated from the environment, where the number 
N may be as small as two or can be several hundreds. 
Larger entities of thousands to a few millions of atoms, 
which are usually referred to as nanoscale or ultrafine 
particles, are excluded from the present review. Mag- 
netic clusters either occur naturally in pure compounds, 
where the magnetic isolation of the clusters is provided 
by nonmagnetic ligands, or are formed artificially, e.g., in 
solid solutions of magnetic and nonmagnetic compounds. 
Prototypes of pure compounds are molecular nanomag- 
nets in which a polynuclear magnetic metal core is em- 
bedded in a diamagnetic ligand matrix. Examples as- 
sociated with cooperative systems are diluted magnetic 
compounds, in which the magnetic ions are randomly 
distributed, so that different types of magnetic clusters 
(7V-mers, TV = 1,2, 3,...) are simultaneously present. 

Many of the characteristic physical properties of mag- 
netic clusters are determined by the magnetic excitation 
spectra which reflect the nature of the fundamental mag- 
netic interactions between the spins in a cluster. The lat- 
ter can be quantitatively accounted for by a spin Hamil- 
tonian in which a bilinear Heisenberg-type exchange in- 
teraction is usually the dominant term, often comple- 
mented by additional terms describing the anisotropic 
and/or higher-order interactions. As long as the num- 
ber N of spins is reasonably small, exact analytical solu- 
tions of the spin Hamiltonian can be obtained. Accord- 
ingly, small magnetic clusters are ideal model systems 
to explore experimentally the limitations of the theoreti- 
cal models. However, the synthesis strategies to produce 
molecular magnetic materials have advanced enormously 
in recent years, making available bounded molecular 
magnetic clusters with the number N of magnetic cen- 
ters varying from one to several dozens; the record cur- 
rently stands at N = 84 in the torus- like molecule Mng4 
( Tasiopoulos et a/1 12004 ). In this new class of magnetic 
materials, now c o mmonl y called molecular nanomagnets 
(iGatteschi et al. I, l2QQ6bh . clusters with more than four 



metal ions are rather the rule than the exception, shift- 
ing the main scientific challenge away from that of study- 
ing the nature of the basic interactions towards that of 
exploring the possible consequences of having such inter- 
actions in a lattice of exchange-coupled spin centers. 

In the past decade, magnetic cluster systems have 
become a topic of increasing interest and relevance in 
condensed matter science. They are not only inter- 
esting in themselves for determining the origin and 
the size of the fundamental magnetic interactions, but 
they have also important applications in both tech- 
nology and science. Molecular nanomagnets are cur- 
rently considered among the most promising candidates 
as the smallest nanomagnetic units ca pable of stor- 
ing a nd processing quantum information (jTroiani et all 
l2QQ5h . In particular, magnetic molecules with large 
spin ground states and negative axial anisotropy or 
single-molecule magnets were found to exhibit outstand- 
ing properties such as slow relaxation of the magneti- 
zation and stepped magnetic hysteresis curves due to 
quantum tunneling of the ma g netization at low temper- 
atures dChristou et all l2000l: IChudnovskv and Tejadal 
119981 : IGatteschi and SessoliL l2003h .~ 

Another emerging field concerns transition-metal per- 
ovskites in which magnetic polarons evolve upon hole 
doping and behave like magn etic nanopart i cles e mbed- 
ded in a nonmagnetic matrix ( Phelan et all 120061 ) . Fur- 
thermore, we mention the remarkable observations of 
quantum phase transitions, Bose-Einstein condensation 
and field-induced three-dimensional magnetic ordering in 
weakly interacting antiferromagnetic dimer compounds 
(jGiamarchi et all 120081 : ISachdevL Il999h as well as the at- 
tractive phenomenon of spin-Peierls d imerization i n bot h 
organic and inorganic compounds ( Brav et~al\ , Il975h . 
Very recently, magnetic cluster systems could also be fab- 
ricated dire ctly on surfaces using sc anning microscope 
techniques (|Hirjibehedin et all , l2QQ6h . Even in biology, 
magnetic metal clusters are important subunits. Polynu- 
clear iron clusters are contained in proteins like adren- 
odoxin and ferredoxin, which are involved in the photo- 
synthetic process to convert light energy into chemica l 
energy based on electron transfer mechanisms (Griffith, 
Il972l ). All these systems have a common property, 
namely the presence of magnetic clusters, whose char- 
acterization is therefore a key issue in both theoretical 
and experimental investigations. 

An important class of magnetic cluster systems con- 
cerns magnetic nanoparticles containing a few hundreds 
of atoms produced by e.g. sputtering, inert-gas conden- 
sation techniques, direct ball milling or microemulsion- 
based syntheses. In principle, the spin dynamics ex- 
hibits many features similar to the molecular nanomag- 
nets, such as superparamagnetic relaxation and quan- 
tized spin-wave states. Pioneering neutron scattering 
experiments i n this area have be en performed for Fe 
nanoparticles ([Hennion et aZ.L fl994l) and for nanocrysta l- 
lites of hematite (jHansen et all 1 1 9971 : iKuhn et all l2006h . 
In practice, the interpretation of experimental data on 
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magnetic nanoparticles is made difficult by line broad- 
ening effects due to variations in size and form, and the 
presence of a variety of surface spin states. Any discus- 
sion of the magnetic properties needs considering these 
aspects, and for that reason the magnetic nanoparticles 
are excluded from the present review, in particular as 
excellent re views of the theoretical and exper i menta l as- 
pects exist (jHendriksen et all 1 1993c iKodamal 1 19991 ). 



II. BASICS 

A. Spin Hamiltonian 

The form of the appropriate spin Hamiltonian of mag- 
netic cluster systems depends on two essentially indepen- 
dent terms: the nature of the interacting systems in the 
absence of interactions, and the physical n ature of th e 
mechanisms responsible for the interactions ( Wolifl97lh . 
The total Hamiltonian describing interacting spins can 
generally be written as 



(i) 



(i) 



Our main objective in this review is to introduce the 
theoretical concepts and the experimental techniques ap- 
plied to the study of magnetic cluster excitations as well 
as to give a snapshot of recent developments achieved 
for the most important classes of materials whose prop- 
erties are largely governed by the presence of magnetic 
clusters. The outline of this review is as follows: Sec. [XT] 
starts with a summary of the basic terms of the underly- 
ing spin Hamiltonians and with a short description of the 
most powerful experimental techniques, followed by rep- 
resentative examples on small magnetic clusters (dimers, 
trimers, and tetramers) in Sec. IIIII to demonstrate the 
enormous amount of detailed physical information result- 
ing from the study of magnetic cluster excitations. Clus- 
ters with a small number N of coupled magnetic ions are 
excellent choices for this purpose, as the basic interac- 
tions are all present in them and can be treated exactly, 
hence a straightforward comparison between theory and 
experiment is possible with little ambiguity. This is no 
longer the case for the emerging field of large magnetic 
clusters discussed in Sec. US in which the number N of 
magnetic centers can be as large as several dozens, so 
that the interpretation and analysis of the experimental 
results require the use of sophisticated tools. In addition, 
novel physical aspects such as complex many-body quan- 
tum states come into play as a consequence of the large 
cluster size. An important subclass of the large magnetic 
clusters will be addressed in Sec. El namely the single- 
molecule magnets in which at low temperatures mag- 
netic hysteresis or slow magnetic relaxation and quantum 
tunneling of the magnetization can be observed. Many 
methodologies applied to the large magnetic clusters have 
their origin in the field of quantum spin systems which 
are summarized in Sec. IVII for cases where the presence 
of magnetic clusters is the most important ingredient to 
understand their quantum spin properties. We conclude 



where ^> and the terms H\ L) may or may not be 
comparable with . The classification of different cases 
depends on the eigenfunctions of which provide a 

basis for the description of and . Since is 
by definition large, we generally consider only the ground 
state at a time, and this may be characterized as one of 
following types: 

Type S: The system has negligible orbital admixtures. 
The ionic spins s$ are good quantum numbers to define 
the spin Hamiltonian. 

Type Q: The orbital angular momentum is quenched, 
thus the spin Hamiltonian can be expressed as for type 
S. However, orbital effects like the ligand field have to be 
considered in . 

Type L: There is orbital degeneracy (or near degener- 
acy) with 21 i + 1 states. The (weak) spin-orbit coupling 

H = A % • Si has to be considered in 

Type J: The spin-orbit coupling is large, thus the to- 
tal angular momentum = 1^ + is a good quantum 
number. The spin Hamiltonian is expressed as for type 
S systems (replacing s$ by 

It is mentioned that type Q and S are also known as 
systems with "orbitally non-degenerate ground terms" 
and type L as systems with "orbitally degenerate ground 
terms" or "first-order orbital angular momentum". In 
the following we write down and for types S 

and Q which constitute the large majority of magnetic 
clusters studied so far. 

A widely used approach to describe the spin inter- 
actions Hij is the Heisenberg-Dirac-Van Vleck (HDVV 
Hamiltonian (jDiracl 119261 : iHeisenberd . 11926c IVan Vied 
[l932|), 



with a brief outlook in Sec. IVIII Given the abundant 
literature on the topic of magnetic cluster excitations, 
the present review is necessarily bound to be incomplete 
in terms of both materials covered and references cited. 
The experimental results were chosen according to their 
didactical suitability to illustrate the concepts, including 
both historical data from the pioneering time and data 
of today's research. 



(2) 



where is the spin operator of the zth ion in the clus- 
ter and Jij the exchange parameter which couples the 
magnetic ions at sites i and j. In the literature the 
HDVV Hamiltonian is often described as — Si • Sj or 
+ ^2 JijSi-Sj, in contrast to the convention adopted here. 
Hence the exchange parameters given in the present work 
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are always adjusted to be in agreement with Eq. (j2j). H 
commutes with the total spin S = s^, thus S and M 
are good quantum numbers and the eigen functions can 
be written as \tSM), where —S<M<S and r stands 
for any other quantum numbers required for distinguish- 
ing the spin multiplets unambiguously. Often labels are 
omitted for convenience. Any anisotropic term added to 
the HDVV Hamiltonian Eq. (|2]) lifts the M-degeneracy 
of the spin states \SM). 

The exchange coupling may not always be isotropic, 
thus we extend Eq. (j2j) to include exchange anisotropy, 

H = — 2 y ^ ( Si X Sj X + J?j ' SiySjy + J i j SizSjz} . (3) 

i<j 

A special case of anisotropic spin-spin coupling is pro- 
vided by the dipole-dipole interaction which is always 
present in addition to the exchange interaction Eq. (j2j), 



2 2 



, (§i * Rij)(Sj • Rij) 



(4) 



where g is the Lande splitting factor, \ib the Bohr magne- 
ton, and Hij = Ri — Rj the vector defining the distance 
between the spins at sites i and j located at R^ and Rj . 
In most magnetic clusters the typical distance between 
metal ions is 3.0-3.5 A yielding g 2 ^ 2 B lR% ^ 0.1 K. The 
dipole-dipole interaction is hence normally much smaller 
than the HDVV interaction and can often be disregarded, 
though in some cases it cont ributes appre c iably to the 
overall magnetic anisotropy ([Abbati et all |2001[ ). An- 
other type of exchange anisotrop y is described b y the 
Dzyaloshinski -Moriya interaction ([Dzvaloshinskil Il958t 
iMorival [i960) which, however, vanishes for the case of 
inversion symmetry, 



(5) 



where d^j is known as the Dzyaloshinski-Moriya vector. 

The HDVV Hamiltonian Eq. (13 is based on the bilin- 
ear spin permutation operator (iHerrind . Il966t ) 



p ij = g 0- + ^ ' § i) ' 



(6) 



A more complete Hamiltonian takes permutations of 
more than two spins into account. The relevant terms 
up to second order (biquadratic terms) are defined by 



P 2 



PijPjk 
PijPkl 



1 

4 



1 + 2si • Sj + (Si • SjY 



(7) 



S ? ; • So 



Sj • s k + (Si • Sj) (sj • s k )] , (8) 
s k • si + (Si • Sj) (gfe • s^)] , (9) 



which refer to two-spin, three-spin, and four-spin inter- 
actions, respectively. 



The Hamiltonian has to be introduced in es- 

sentially two cases. For ^ 1/2 systems, single- ion 
anisotropy must be considered, which for the case of axial 
anisotropy reads 



1 



Si(Si + 1) 



and for planar anisotropy we have 



(10) 



(ii) 



where often Di = D and/or E{ = E for all sites. In cases 
with Si > 2 also higher-order anisotropy terms 

H = ^B'iOlis^ + BlOl^ + BiMsi) (12) 



may have to be included, where 0™(si) are Stevens op- 
erato r equivalents bu ilt up by fourth-order spin opera- 
tors (iHutchingsl . fl964) . Finally, the action of an external 
magnetic field B is described by 



H 



Si = g/i^B • S. 



(13) 



For type S and Q systems, the HDVV Hamiltonian 
dominates usually over the anisotropic terms in the total 
spin Hamiltonian, and a first-order perturbation treat- 
ment of the anisotropy is an excellent startin g poin t 
(strong-exchange limit) ([Bencini and Gatteschil |l990). 
The energy spectrum is then structured into spin mul- 
tiplets with a definite value of S for each of them, and 
the eigen functions are well described by the \rSM) spin 
functions. The energies of the spin multiplets are gov- 
erned by the HDVV Hamiltonian (exchange splitting), 
and each spin multiplet is further split by the magnetic 
anisotropy (anisotropy splitting or zero-field splitting, 
ZFS). The possible transitions may be distinguished into 
inter- multiplet (AS ^ 0) and intra- multiplet (AS = 0). 
Unless stated otherwise, the strong-exchange limit is al- 
ways assumed. 



B. Experimental Techniques 

The spin interactions discussed in the preceding Sec- 
tion give rise to discrete energy levels and wave functions 
which can be determined by a variety of experimental 
methods. The most powerful techniques are certainly 
spectroscopic methods such as inelastic neutron scatter- 
ing (INS) and optical spectroscopies, which allow a di- 
rect determination of the spin states. Some information 
about the spin states can also be obtained by resonance 
techniques, e.g., by paramagnetic resonance experiments 
(EPR). Information on the spin states is also contained 
intrinsically in the thermodynamic magnetic properties, 
however, an extraction of reliable parameters is not al- 
ways possible due to the integral nature of these proper- 
ties. 
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In the following the two spectroscopic methods mainly 
applied to the study of magnetic cluster systems are 
briefly introduced. These include INS and optical spec- 
troscopies which both have their merits and should be 
considered as complementary methods. Optical spectro- 
scopies, on the one hand, can be applied to very small 
samples of the order of 10 /xm 3 ; they provide highly re- 
solved spectra so that small line shifts and splittings can 
be detected, and they cover a large energy range so that 
inter-multiplet transitions can easily be observed. Neu- 
tron scattering, on the other hand, is not restricted to 
particular points in reciprocal space, i.e., interactions be- 
tween the spins can be observed through the wave vector 
dependence, the peak intensities can easily be interpreted 
on the basis of the wave functions of the spin states, and 
data can be taken over a wide temperature range which 
is important when studying linewidth phenomena. INS 
as the most widely used spectroscopic technique is de- 
scribed below in detail, followed by short descriptions of 
optical spectroscopies and EPR techniques as well as by 
a summary of the thermodynamic magnetic properties. 



1. Inelastic Neutron Scattering 

The principal aim of an INS experiment is the determi- 
nation of the probability that a neutron which is incident 
on the sample with wave vector k is scattered into the 
state with wave vector k'. The intensity of the scattered 
neutrons is thus measured as a function of the momen- 
tum transfer 



ftQ = ft(k-k'), 



(14) 



where Q is known as the scattering vector, and the cor- 
responding energy transfer is given by 



2m V 



(15) 



where m is the mass of the neutron. Eqs. (fT4|) and ([T5]) 
describe the momentum and energy conservation of the 
neutron scattering process, respectively. For |k| = |k'| we 
have from Eq. ([T5]) huj = 0, i.e., elastic scattering. For 
inelastic scattering, Q can be decomposed according to 
Q = G + q, with a reciprocal lattice vector G and a wave 
vector q. INS experiments thus allow us to measure the 
magnetic excitation energy at any predetermined point 
in reciprocal space , most conveniently by triple-axis crys- 
tal spectrometry (iBrockhousd . 119551 ). In extended sys- 
tems this yields the dispersion relation huj(q). In mag- 
netic clusters the excitations are dispersion-less but the 
scattering intensity shows a characteristic dependence on 
momentum transfer (Q dependence, vide infra). For INS 
experiments on polycrystalline samples various types of 
time-of-fli ght (TOF) spectrom eters are usually more ap- 
propriate ( Furrer et a/.Ll2QQ9l ). 

The neutron scattering probability for magnetic cluster 
excitations can be derived from the master formula for 



magnetic scattering ( Lovesevl 1 19871 ): 
d 2 cr ^ / Q a Q^ 

a(3 



dftduj 



= C(Q)£(^-^)s Q/3 (Q, w ), (16) 



where 



C(Q) = ( 7 ro) 2 yF 2 (Q)e-™ (17) 
and S a p(Q,u)) is the magnetic scattering function 
S aP (Q,u) = ^e lQ - R -5> A (A|s iQ |A'> 

xiX'lsj^Sihu + Ex-Ex,). (18) 

Herein is 7 = -1.91, r = 0.282xl0 -12 cm the classi- 
cal electron radius [(jro) 2 = 0.29 barn], F(Q) the di- 
mensionless magnetic form factor defined as the Fourier 
transform of the normalized spin density associated with 
the magnetic ions, exp[— 2W(Q)] the Debye- Waller fac- 
tor, and a, (3 = x,y, z. |A) denotes the initial state of the 
scatterer, with energy E\ and thermal population factor 
Px [Eq. dUl)]; and |A') its final state with energy Ey . 

The essential factor in the cross section is the magnetic 
scattering function S a p(Q,uj) which will be discussed in 
more detail below. There are two further factors which 
govern the cross section for magnetic neutron scattering 
in a characteristic way: Firstly, the magnetic form fac- 
tor F(Q) which usually falls off with increasing modulus 
of the scattering vector Q. Secondly, the polarization 
factor (5 a p — QaQp/Q 2 ) tells us that neutrons can only 
couple to magnetic moments or spin fluctuations perpen- 
dicular to Q which unambiguously allows to distinguish 
between different polarizations (transverse and longitu- 
dinal) of spin excitations. 

The magnetic scattering function S a p(Q,uj) contains 
two important terms: Firstly, the structure factor 
exp(zQ • Hij) which directly reflects the geometry of the 
cluster; secondly, the matrix elements (\\si a \\'} which 
determine the strength of the transition |A) — ^ |A') as 
well as corresponding selection rules. 

For magnetic clusters we describe the eigen state |A) by 
\tSM). The matrix elements can then be calculated by 
introducing irreducible tensor operators (ITOs) T^{si) of 
rank 1, which are related to the spin operators Si a : 

To(3i) = 8 iz , fl^Si) = {Six ± S iy ) • (19) 

In the HDVV model the states \SM) are degenerate with 
respect to the magnetic quantum number M, so that 
Eq. (jTHJ) has to be summed over both M and M' . Using 
the Wigner-Eckart theorem 



{SM\fi( Si )\S'M') = (-1) 



S-M 



S 1 S' 

-M q M' 

x {S\\f\ Si )\\S') 



(20) 



6 



we find 



J2(SM\f^ Si )\S'M')(S'M'\f q \( Sj )\SM) 



MM' 



= -{S\\f\s i )\\S')(S'\\f 1 (s j )\\S). (21) 

The two-row bracket ( ...) in Eq. (J^UJ) is a Wigner-3j sym- 
bol (jRotenberg et aZX Il959h . It vanishes unless 



AS = S'-S = 0,±1, 
AM = M' — M = 0,±1, 



(22) 
(23) 



which establish the selection rule for INS in spin clus- 
ters. Thus, INS experiments allow us to detect not only 
splittings of individual spin multiplets (AS = 0), simi- 
lar to EPR experiments (Sec. III.B.3|L but also splittings 
produced by magnetic interactions (AS = ±1). The eval- 
uation of the reduced matrix elements on the right-hand 
side of Eq. (|2T|) depends on the details or many-body 
structure of the spin functions \rSM). 

Equations ([T8]) -(|2T ]) strictly apply to magnetic cluster 
systems of type S and Q. A theoretical treatm ent of the 
scatt ering by L- and J-type ions was given in IJohnstonL 
1966. However, the calculation is complicated, and we 
simply quote the result for Q — >• 0. In this case the 
cross section measures the magnetization, jli = —/jLb(U + 
2§i), i.e., a combination of spin and orbital moments that 
does not allow their separation. This clearly contrasts to 
magnetic scattering by X-rays. For INS an approximate 
result can be obtained for modest values of Q. We replace 
the spin operator by 



where 



1 



+ + !) + *(* + !) 



2j;(ji + l) 



(24) 



(25) 



is the Lande splitting factor. 

If uj is a positive quantity in the scattering function 
S a p(Qi the neutron loses energy in the scattering pro- 
cess and the system is excited from the initial state A 
which has energy Tiuj less than the final state A' . Con- 
sider now the function S a p(Q, —uj) where uj is the same 
positive quantity. This represents a process in which the 
neutron gains energy. The transitions of the system are 
between the same states as for the previous process, but 
now A' is the initial state and A is the final state. The 
probability of the system being initially in the higher 
state is smaller by the factor exp(— fiu/hsT) as compared 
to its probability of being in the lower energy state, hence 

S a p(Q,-uj) = exp(--^-)S a/3 (Q,uj), (26) 

which is known as the principle of detailed balance. 
Eq. ([26]) has to be fulfilled in experimental data taken 



in both energy-gain and energy-loss configurations which 
correspond to the so-called Stokes and anti-Stokes pro- 
cesses, respectively. 

Using the integral representation of the S function the 
scattering function S a p(Q,uj), Eq. (fT8]h transforms into 
a physically transparent form: 



*3 

+oo 



x / ((s ia (O)s j0 (t)))e- wt dt, (27) 

J — oo 

where ((si a (0)sjp(t)}} is the thermal average of time- 
dependent spin operators, or th e van Hove pair corre- 
lation function ([Van Hovel Il954f ) for spins. A neutron 
scattering experiment measures the Fourier transform of 
the pair correlation function in space and time, which is 
clearly what is needed to describe a magnetic system on 
an atomic scale. 

The van Hove representation of the cross section in 
terms of pair correlation func tions is relate d to the 
fluctuation-dissipation theorem (|Lovesevl 1 19871 ): 



S a p(Q,uj) 



1 — exp(- 



Tiuj 
kB~T J 



(28) 



Physically speaking, the neutron may be considered as a 
magnetic probe which effectively establishes a frequency- 
and wave vector-dependent magnetic field B(Q,cj) in 
the sample, and detects its response, the magnetization 
M(Q,co>), to this field given by 

M a (Q,uj) = J2xMQ:0j)Bp(Q : uj), (29) 



where Xap(QiW) is the generalized magnetic susceptibil- 
ity tensor. This is really the outstanding property of 
the neutron in a magnetic scattering measurement, and 
no other experimental technique is able to provide such 
detailed microscopic information about magnetic com- 
pounds. 

For polycrystalline material Eq. (fl6|) has to be aver- 
aged in Q space, whic h in zero magnet ic field can be 
performed analytically (Wald mannl . |2QQ3[ ): 



d 2 a 
dftdu 



C(Q) ^ J2k j ^ QRi ^ ' sj+MQRij) 



XX' 



^^Tq* (Rij)Tq (siSj)]S(huj + E\ — E\>). (30) 



th 



jo(x) and j2(x) are the spherical Bessel functions of 
and 2 nd order, and s, = ((X\s ix \X'), <A|^|A'), (\\s iz \\')). 
For an isotropic spin cluster described by only the HDVV 
Hamiltonian the 2 nd order term vanishes: 



sm(QRij) 



dVtduj 



XX' 



QRij 

x(S"||T 1 (s i )l|5)<5(?iw + J E A 



(silicons') 

~Ey). (31) 
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The Bessel function sm(QRij)/ (QRij) is responsible 
for a characteristic oscillatory Q dependence of the 
INS intensity, wh i ch is often very helpfu l in analysis 
( Furrer and GudelL 119771 : IWaldmannl . [2003) . Also, a use- 
ful rule of thumb is inferred: For Q — » the scattering 
intensity of AS = ±1 transitions drops to zero, while for 
AS = transitions it becomes maximal. 

Analytical results for the INS cross section were 
derived for so me cases, i.e., f o r dimers, t rimers , 
and tetramers ( Furrer and Glide! . Il979l : iGiidel et all 
Il979l : lHaraldsen et all l2QQ5h and a pentamer a nd hex- 
amer (jHaraldsenl . 1201 it lHaraldsen et all l2QQ9h . Ex- 
plicit expressions for Eq. (| 30l) can be found in, e.g. 



IWaldmann and GudelL I200E 

In practical applications of INS to magnetic cluster 
compounds the huge incoherent neutron scattering con- 
tribution of hydrogen can easily prevent observation of 
the magnetic cluster excitations. Removing or reducing 
the hydrogen content by e.g. deuteration or fluorination 
is of course the best solution, but this is often prohibitive, 
in particular in molecular clusters. Fortunately, at trans- 
fer energies from ca. 0.1 to 3 meV a window exists with 
comparatively small hydrogen scattering. This is rele- 
vant because otherwise most studies on the cluster exci- 
tations in the molecular nanomagnets for instance would 
not have been possible. 

Another point to be considered is the non-magnetic 
scattering from the lattice. Besides the "standard tricks" 
for identifying the nature of INS features, such as in- 
specting the temperature and Q dependencies, a Bose- 
correction analysis is often helpful. Here, INS data 
recorded at sufficiently high temperature are scaled by 
the Bose factor and then compared to the data at lower 
temperatures. At high temperatures, where a large part 
of the energy spectrum is accessed and the magnetic scat- 
tering intensity spread out over essentially all frequencies, 
the measured spectrum may reflect the lattice scatter- 
ing, whose temperature dependence is governed by the 
Bose factor [1 — exp(fkj / Tc^T)] -1 (neutron-energy loss). 
Accordingly, the Bose-scaled high-temperature data can 
estimate the lattice contribution at lower temperatures. 
Often this works well, especially in large magnetic clus- 
ters with a dense higher-lying energy spectrum, and al- 
lows an unambiguou s identification of magnetic peaks 
([Preiser et all l2010bl : lOchsenbein et a/.l . l2008h . 



2. Optical Spectroscopies 

Optical spectroscopies cover a large range of wave- 
lengths of light. Individual spectrometers are specialized 
devices that focus on particular parts of the electromag- 
netic spectrum produced by lamps, lasers or synchrotron 
sources. They therefore e xist in a wide var iety of types 
for different applications ([Tkachenkol 120061 ). One major 
type of optical spectroscopy is absorption spectroscopy, 
where the absorbance of a system is determined by mea- 
suring the photons which pass through (transmittance 



I 



FIG. 1 (Color online) Sketch of the processes relevant for 
optical spectroscopies. A: absorption; L: luminescence; R: 
radiationless transition. 



spectrum). Another important type is emission or lu- 
minescence spectroscopy. When a system is excited by 
an outside energy source such as light, it eventually re- 
turns back to the ground state by releasing the excess 
energy either as radiation-less transitions or in the form 
of photons as illustrated in Fig. [TJ 

A variant of the absorption and luminescence spec- 
troscopies, associated with initial transitions to dis- 
crete ex c ited e nergy states, is the Raman spectroscopy 
( Larkini l201lf ) where the system is excited to a vir- 
tual energy state and then quickly relaxes back to a 
ground-state level. Unlike a luminescence process, Ra- 
man scattering involves no transfer of electron pop- 
ulation to the intermediate state. Several variations 
of Raman spectroscopy have been developed in order 
to enhance the se nsitivity [e.g., surface-enha nced Ra- 
man spectroscopy (Lombardi and Birke, 2008) and res- 
onance Raman spectroscopy ([Chao et all Il976l )] as well 
as to improve the spatial resolution [Raman microscopy 
(jTurrell and CorsetL fl996l )]. 

Optical spectroscopies are governed by the energy and 
momentum conservation laws as in neutron spectroscopy, 
see Eqs. ([T4]) and ([T5]) . However, as the photon wave vec- 
tor is about 10 3 times smaller than a typical reciprocal 
lattice vector, only excitations close to the center of the 
Brillouin zone are observed. The calculation of intensities 
of the observed transitions is a non-trivial task. This is 
in contrast to neutron spectroscopy where the intensities 
of spin excitations are directly proportional to the square 
of the magnetic dipole matrix elements, Eq. (fT8|) . Since 
optical spectroscopies often involve intermediate states 
which are not known, approximate models have to be em- 
ployed for the calculation o f transition matrix elements 
(|Lovesev and Collinsl 1 19961 ). 

The polarization of light has great importance particu- 
larly when anisotropic systems are studied. The specific 
polarization of both the exciting and the emitted light 
can be exploited to obtain extra information concerning 
the line identification from the observed energy spectra. 
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More specifically, electronic states with transition dipole 
moments perpendicular to the electric field orientation 
will not be excited. 



3. Electron Paramagnetic Resonance 

In EPR spectroscopy the absorption of a radio- 
frequency (rf) magne tic field B r ^ by a mag netic sys- 
tem is measured ( Abragam and Bleanevl .fl986). Absorp- 
tion can occur whenever the energy hv of the radiation 
matches the energy difference of two eigen states |A) and 
|A'>, 



±hv, 



(32) 



and the absorbed power P is calculated in linear response 
theory to 

P = C(a;)^^^^^p A (A|^|A')(A'| S >|A) 

a/3 ij XX' 

x5(Huj + E x -E x ,), (33) 

with C(uj) = uj[1 - exp(huj/k B T)}(g 2 fi 2 B )/(Sh7r). For 
magnetic cluster systems with eigen states \tSM), 
Eq. ([33]) can be further evaluated and the EPR selection 
rules 



AS = 0, 
AM = ±1, 



(34) 
(35) 



established from Ei(A|^|A / ) = (rSM\S a \r' ' S> ' M') . 

It follows that EPR spectroscopy is a very direct 
method to determine anisotropics of the g factor by align- 
ing the external magnetic field B along different direc- 
tions. Similarly, anisotropics of the form defined by, e.g., 
Eqs. ([TO]) and ([TT]) . which lift the degeneracy of a par- 
ticular spin multiplet, can also be determined from the 
positions of the lines in the EPR spectra. On the other 
hand, the exchange splittings or parameters Jij are not 
directly attainable, but they can be estimated from the 
temperature variation of the signal intensities which fol- 
low the Boltzmann populations of the energy levels in- 
volved, or i n some fortunate cas es through the S'-mixing 
mechanism ([Wilson et all l2006f ) . Finally we point to the 
distinctive hyperfine structure superimposed on an EPR 
spectrum for systems with non-zero nuclear spin quan- 
tum numbers. 

It is instructive to compare Eq. ([33]) to the correspond- 
ing INS formula Eq. JE) with Eq. CEJ). The main dif- 
ference lies in the structure factor exp(Q • IUj) which is 
1 in case of EPR corresponding to Q = in INS. There- 
fore, EPR affords the detection of exactly those magnetic 
transitions which have INS intensity at Q — >• 0, which in 
the HDVV model are the AS = transitions. Physically 
speaking, in contrast to a neutron the applied radio fre- 
quency establishes a frequency-dependent but spatially 
homogeneous magnetic field, B(Q,cj) with Q = 0. 



Modern EPR spectroscopy techniques permit a large 
combination of frequency and magnetic field values ex- 
tending up to the THz regime and 25 T, respectively. 
In principle, EPR spectra can be generated by either 
varying the frequency v while holding the magnetic field 
constant, or doing the reverse. In commercial EPR in- 
struments it is the frequency which is kept fixed, and 
typical frequencies are X-band (10 GHz) and Q-band 
(35 GHz), but W-band (95 GHz) is also available. How- 
ever, the EPR techniques have progressed enormously, 
and multi-frequency high-field EPR and frequency do- 
main magnetic resonance spectroscopy (FDMRS) exper- 
iments are routinely undertaken in various laboratories. 
A recent development are THz EPR experiments us- 
i ng radiation from synchrotron sources. Fo r revie ws see 
(Gattesc hi et all l2006at Ivan Slageren et all l2QQ3h . 



4. Thermodynamic Magnetic Properties 

The thermodynamic magnetic properties depend ex- 
plicitly upon both the energies E\ and the eigen functions 
| A) of the spin excitations. Based on general expressions 
of statistical mechanics for the Gibbs free energy F and 
internal energy [/", 



F = -fc B TlnZ, 



(36) 
(37) 



we obtain, with the Zeeman term as in Eq. (fT3"]) . the 
magnetization M a , magnetic susceptibility Xaa, entropy 
S, and Schottky heat capacity cy. 



dF 



^ Pa (A| So I A), 



8M a (gii B y 



dB a k B T 

x 5>a(A|^|A)- (5> A (A|£ a |A) 



(38) 
(39) 



S 



cv 



f^T) =k B [\nZ 



(40) 



dU 
dT 



knT 2 



^PxEl-fepxE* 

A V A 



(41) 



Here, Z and p\ are the partition function and Boltzmann 
population factor, respectively: 
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For a system with magnetic anisotropy also the magnetic 
torque r appears as a useful thermodynamic quantity: 

r = ^, t = M x B, (43) 

where denotes the rotation angle around the torque 
axis [often the definition r = —dF/dcj) is found, our con- 
vention is consistent with the usual parametrization of 
magnetic field, e.g., B = £?(sin#, 0, cos#)]. 

For T — >• 0, the free energy reduces to the ground- 
state energy Eq and the magnetization to the field deriva- 
tive M a = —dEo/dB a . As function of field the ground 
state often undergoes level crossings at characteristic 
fields, which can be detected at very low temperatures 
as steps in the field-dependent magnetization (torque) 
curves. The characteristic fields allow insight into the 
magnetic excitation spectrum in a cluster, and low- 
temperature high-field magnetization (torque) measure- 
ments represent an importan t experimental technique 
(|Shapira and Bindilattii |2002[ ) , though the level crossing 
can also be detected by ot her technique s ^ e.g., proton nu- 
clear magnetic resonance (jjulien et a/1 Il999h - In order 
to check the reliability of the model parameters derived 
from spectroscopic data, it is however generally useful to 
compare the calculated thermodynamic magnetic prop- 
erties to corresponding experimental data. 

III. SMALL MAGNETIC CLUSTERS 

The aim of this section is to demonstrate how the var- 
ious interactions introduced in Sec. III. Al manifest them- 
selves for different experimental techniques. The pre- 
sented examples will be restricted to small clusters built 
up by N < 4 coupled magnetic ions, for which the under- 
lying models can be treated exactly, since only a small 
number of interactions are present and cooperative ef- 
fects do not occur, so that a straightforward compari- 
son between theory and experiment is possible with lit- 
tle ambiguity. The examples cover magnetic clusters 
which naturally occur in pure compounds as well as clus- 
ters which are artificially formed in solid solutions of 
magnetic and nonmagnetic compounds. Ideal examples 
of pure compounds are molecular transition-metal com- 
plexes, in which a polynuclear metal core is embedded 
in a diamagnetic ligand matrix. Information more di- 
rectly associated with cooperative systems results from 
diluted magnetic compounds, in which the magnetic ions 
are randomly distributed, so that different types of clus- 
ters (7V-mers, 7V = 1,2,3,...) are simultaneously present. 
Among the myriads of small magnetic cluster systems 
studied up to the present we will choose as a represen- 
tative of the pure compounds the dimeric chromium sys- 
tem [(NH 3 )5CrOHCr(NH3) 5 ]Cl5-H 2 0, which is the first 
magnetic cluster system investigated by INS. The class 
of magnetically diluted systems was pioneered in INS ex- 
periments carried out for the compound KMn^Zni-^Fa 
( Svensson et all Il978h and will be exemplified here by 



the compound CsMn^Mgi-^Bra. The chapter ends with 
further insights on particular physical aspects resulting 
from magnetic cluster excitations on some other com- 
pounds. 



A. The dimeric chromium compound 
[(NH 3 )5CrOHCr(NH 3 )5]CI 5 H 2 

1. Energy levels 

The simplest magnetic cluster system is the dimer (two 
coupled spins si and §2) for which the HDVV Hamilto- 
nian Eq. ([2]) simplifies to 

#=-2Jsi-s 2 . (44) 

Assuming identical magnetic ions (si = 52 = s) the 
eigenvalues of Eq. (|44|) are 

E(S) = -J [5(5 + 1) - 2s(s + 1)] , (45) 

with < 5 < 2s. The energy splittings defined by 
Eq. (f45|) satisfy the Lande interval rule 

E(S) - E(S - 1) = -2JS. (46) 

For Cr 3+ dimers with 5 = 3/2 the separation between 
the ground-state levels will be 2 J, 4J, and 6 J (5 = to 
5 = 3), with the state 5 = being the lowest in case 
of AFM exchange J < 0. Observed deviations from the 
Lande interval rule are often attributed to the presence 
of biquadratic exchange, 

H=-K(s l -s 2 ) 2 . (47) 

Combining Eqs. (|44j) and (|47j) yields the modified eigen- 
values 

E(S) = -Jt)-\k^\ (48) 
rj = 5(5 + 1) -25(5 + 1). (49) 

2. Structural and magnetic characterization 

The compound [(NH 3 )5CrOHCr(NH3) 5 ]Cl5-H 2 was 
characterized by X-ray diffractio n, EPR, and magnetic 
susceptibility measurements ([Veal et all Il973h . The ma- 
terial crystallizes in the tetragonal space group P42/mnm 
with four formula units in a cell of dimensions a = 
16.259(7) A and c = 7.411(7) A. The two Cr 3+ ions 
are coupled by superexchange via a Cr-O-Cr bridge with 
a bridging angle of 165.6(9)° and a Cr-0 distance of 
1.94(1) A. The analysis of the X-ray data requires two 
inequivalent positions of Cr 3+ dimers in the unit cell. 
EPR measurements gave a negligibly small upper limit 
of D < 0.002 meV for the single- ion anisotropy defined by 
Eq. ([TO]) . The magnetic susceptiblity was measured for a 
polycrystalline sample as shown in Fig. [2j The data were 
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J=-1 .95 meV, K=0 

J=-1.83 meV, K=0.07 meV 



50 100 150 200 250 300 

T[K] 



FIG. 2 (Color online) Temperature dependence of the 
magnetic susc eptibilty of [(NH 3 ) 5 CrOHCr (NH 3 ) 5] Cl 5 - H 2 O . 
Adapted from IVeal et all Il973l . 
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FIG. 3 Polarized luminescence spectrum of 
[(NH3)5CrOHCr(NH 3 )5]-Cl5-H 2 taken at T = 7 K. 
The correspon ding transition diagram is shown at the top. 
Adapted from iFerguson and Giidell |1973. 



analyzed according to Eq. ([39]) with g = 1.99 resulting 
from the EPR experiments. The agreement between the 
observed and calculated data is slightly improved when 
in addition to the Heisenberg exchange Eq. ([44]) a bi- 
quadratic term Eq. ([47)) is included. 



3. Optical spectroscopies 

Optical spectroscopies were appli ed to single crystals of 
i[NH3)5CrOHCr(NH 3 )5]Cl5-H 2 (iFerguson and Giidell . 
Il973l ). Both polarized absorption and polarized lumi- 
nescence spectra provided well resolved lines from which 
the ground-state level scheme could be directly deter- 
mined. Figure [3] shows a representative polarized lumi- 
nescence spectrum with two sets of transitions (A',B',C) 
and (A",B",C") reflecting the presence of two inequiva- 
lent dimer sites. The emission starts from an excited 
state with 5 = 2. The appearance of three lines for each 
of the transitions indicates that the selection rule AS = 
is not exact, but transitions also occur for AS = ±1 (with 
much smaller intensities) due to the spin-orbit interac- 
tion. The luminescence spectrum accurately determines 
the separations between the ground- state levels S = 1,2, 
and 3, and the separation between 5 = and 5 = 1 was 
taken from the absorption spectrum. The ground-state 
level scheme slightly deviates from the Lande interval 
rule, so that the data analysis was based on Eq. ([48]) . 
The resulting bilinear and biquadratic exchange param- 
eters J and K are listed in Table [I] The luminescence 
spectrum is strongly temperature dependent, and it is 
completely quenched at room temperature. 



4. Inelastic neutron scattering 

For the analysis of the neutron data we adjust the mag- 
netic scattering function Eq. (fT8|) to the dimer case. We 



start from the reduced matrix elements introduced in 
Eq. (j2T|) . Since T 1 ^) operates only on the zth ion of 
the coupled system, the reduced matrix elements can be 
further simplified: 



(S\\f\ Sl )\\S') 



(-1) 2S+S+1 VW 
JSS> 1 

y s s s 

(-lf- s (S\\f\ Sl )\\S'), 



1)(2S' + 1) 

(slllT^HIs), (50) 



(S\\T\s 2 )\\S') 
(sWlf^sMs) = v/ s ( s +l)(2 S +l). 



(51) 
(52) 



The two -row bracket {••• I in E q. ([50]) is a Wigner-6j 
symbol ( Rotenberg et all Il959[ ) which vanishes unless 
AS = 0,±1. From Eqs. ([20j) and ([50]) the INS selec- 
tion rules AM = 0, ±1 and AS = 0, ±1 are recovered. 
By making use of the symmetry properties of the reduced 
matrix elements defined by Eq. ([50]) , we find the following 
cross section for the dimer transition |5) — >• |5'): 



d 2 a 
dQduj 



= C(Q) 



exp[ 



' k B T\ 



E 



x- [l + (-l) AS cos(Q-R)] (SWfHs^WS'f 



iS[nu + E{S) - E(S% 



(53) 



where R = Ri — R2 is the vector defining the intra-dimer 
separation. The structure factor [1 + ( — 1) AS cos(Q • R)] 
is a powerful means to unambiguously identify dimer ex- 
citations from other scattering contributions due to its 
characteristic oscillating behavior. 

For a polycrystalline material Eq. ([53]) has to be aver- 
aged in Q space: 



d 2 a 
dQduj 



C(Q) 



exp[ 



k B T> 4 



As sm(QR) 



QR 

x(^||r 1 (s 1 )||^') 2 <5[^ + E {S) ~ E(S') 
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(a) 



[(ND 3 ) 5 Cr-OD-Cr(ND 3 ) 5 CI 3 ]D 2 



T = 10 K 
T = 30 K 
T = 165 K 
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FIG. 4 (Color online) (a) Energy spectra of neutrons scat- 
tered from deuterated [(NH3)5CrOHCr(NH 3 )5]Cl5-H 2 0. (b) 
Q dependence of the intensity of the |0) — > |1) transition ob- 
served at T = 4.2 K show n in panel (a). Adapted from (a) 
IGiidel et all Il98ll and (b) iFurrer and GiidelL Il977l . 



(54) 



The polarization and structure factors combine to the in- 
terference factor [1 + (-1) A5 sm(QR)/(QR)], which pro- 
duces a damped oscillatory Q dependence of the intensi- 
ties. 

Figure SJa) shows the temperature dependence 
of neutrons scattered from a polycrystalline sam- 
ple of deuterat ed [^N H^CrOHCr ( NH^k ] CI5 • H 2 O 

(jFurrer and Giidell . 119771 : IGiidel et all Il98lh . which 
demonstrates the successive appearance of the excited- 
state transitions with increasing temperature. The data 
confirm the ground-state splitting pattern sketched 
on top of Fig. [3l The resulting parameters based on 
Eq. (US)) are listed in Table HI The oscillatory behavior 
of the intensities vs the modulus of the scattering vector 
Q predicted by Eq. ([54]) is nicely verified as shown in 
Fig. mb). 



5. Comparison of different experimental techniques 

Table U lists the results obtained by different 
experimental techniques presented in the preceding 
subsections. From the EPR measurements the 
anisotropic magnetic effects associated with the com- 
pound [(NH 3 )5CrOHCr(NH3) 5 ]Cl5-H 2 were established 
to be negligibly small. This was verified in subsequent 
light and neutron spectroscopic investigations, which did 
not give evidence for any anisotropy-induced line split- 
tings. Due to the excellent energy resolution, light spec- 
troscopies provide at low temperatures rather precise 
spin coupling parameters, but information on their tem- 
perature dependence is severely hampered because of sig- 
nal quenching. This is not the case for inelastic neutron 
scattering, which gives evidence for a strong temperature 
dependence of the bilinear exchange parameter J of the 
order of 15% upon heating from 7 K to room temper- 
ature. The analysis of the magnetic susceptibility data 
thus results in a temperature-averaged parameter J, and 
it overestimates the biquadratic coupling parameter K 



TABLE I Coupling parameters J and K of 
[(NH3)5CrOHCr(NH 3 )5]Cl5-H 2 determined by different 
experimental techniques. 



Technique 


T[K] 


J [meV] 


K [meV] 


mag. suscept. a 


7-300 


-1.83 


0.07 


light spec. 6 


7 


-1.91(1) 


0.02(1) 


INS C 


30 


-1.88(5) 


0.03(2) 


INS C 


165 


-1.83(8) 


0.02(4) 


INS C 


293 


-1.74(14) 


0.03(7) 



Veal et all fl973l. 



a Ref. 
6 Ref. 

c Ref . IGiidel et al, 



Ferguson anc Giidel. ll973l . 



1981. 



by a factor of 2. 

B. Manganese iV-mers in CsMn x Mgi_ x Br3 

1. Structural and magnetic characterization 

Solid solutions of composition CsMn^Mgi-^Bra are 
ideal model systems for various reasons. Both CsMnBr3 
and CsMgBr3 crystallize in the hexagonal space group 
P6s/mmc, and their unit cell parameters are almost 
equal: a = b = 7.609(15) A, c = 6.52(5) A for CsMnBr 3 
(jGoodvear and Kennedy! [l972| ) and a = b = 7.610(2) A , 
c = 6.502(2) A for CsMgBr 3 (|McPherson et all Il980h . 
The structure consists of chains of face-sharing MBr6 oc- 
tahedra parallel to the c axis, where M is Mn 2+ (s — 
5/2) or Mg 2+ (diamagnetic). Spin- wave experiments 
gave evidence for a pronounced one-dimensional mag- 
netic behavior with the intra-chain exchange interaction 
exceeding the inter- ch ain exchange in t eract i on by three 
orders of magnitude ([Breitling et al I Il977t [Falk et al, 
Il987bh . All the Mn 2+ clusters in the mixed compound 
CsMn^Mgi-^Bra are thus linear chain fragments with 
composition MnArBr 3 ( A r+i) (N = 1, 2, 3, ...) oriented par- 
allel to the c axis. The Mn 2+ clusters are statistically 
distributed with the probability pn(x) for N-mer forma- 
tion given by 



Pn(x) = (1 — x)x 



N-l 



(55) 



For Mn 2+ concentrations x < 0.05, monomers and dimers 
dominate, but for x > 0.05 trimers, tetramers, etc. have 
to be considered. For the ground-state splitting pattern 
of dimers we refer to Sec. IIII.A.ll The energy levels of 
linear trimers and tetramers are summarized below. 



2. Energy levels of linear trimers and tetramers 



The HDVV Hamiltonian of a linear trimer is defined 



by 



H = -2 J (§1 • s 2 + s 2 • §3) - 2 J'si • s 3 . 



(56) 



It is convenient to introduce the spin quantum numbers 
Si 3 and S resulting from the spin coupling scheme de- 
fined by the vector sums S13 = §1 + S3 and S = S2 + S13 
with < S13 < 2s and |Si 3 - s\ < S < (S13 + s), re- 
spectively, assuming si = S2 = S3 = s. The trimer states 
are therefore defined by IS13SM), and their degeneracy 
is (2S + 1). With this choice of spin quantum numbers, 
the Hamiltonian Eq. ([56]) is diagonal and the eigenvalues 
can thus readily be derived as 

£(Si 3 ,S) = -j[S(S +1) - S 1S (S 13 + 1) - s(s + l)] 
-J'[S 13 (S 13 + l)-2s(s + l)]. (57) 

The HDVV Hamiltonian of a linear tetramer is given by 

H = -2J(si • s 2 + s 2 • s 3 + s 3 • s 4 ) 

-2 J' (§1 • s 3 + s 2 • s 4 ) - 2 J"si • s 4 . (58) 

To solve Eq. ([58]) . the total spin S = si + s 2 + S3 + S4 is 
still a good quantum number, but for a complete charac- 
terization of the tetramer states additional intermediate 
spin quantum numbers are needed, e.g., S12 = Si+s 2 and 
S34 = S3+S4 with < Si 2 < 2s and < S34 < 2s, respec- 
tively. The total spin is then defined by |Si 2 — S34I < S < 
(S12 + S34), and the basis states are the wave functions 
IS12S34SM). There is no spin coupling scheme which 
results in a diagonal Hamiltonian matrix, so that the 
eigenvalues of Eq. ([58]) have to be calculated numerically 
or by spin-operator techniques ([Juddl 1 19631 ) . 

3. Electron paramagnetic resonance 
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FIG. 5 Q-band EPR spectrum of Mn 2+ ions in CsMgBr3 at 
T — 77 K. Adapted from iMcPherson et all Il974l . 



o-6 - Mn 2+ in CsMgBr 3 
g=2.004, D=0.0115 meV 



0.4 - 




0.5 1.0 1.5 2.0 



B // c [T] 

FIG. 6 (Color online) Field dependence of the ground-state 
levels of Mn 2+ ions in CsMgBr3 . The calculations are based 
on Eqs. (|10f) and (|13p . with D > 0. The double arrows mark 
the resonances A to E observed in the EPR spectrum of Fig. 



Single crystals of CsMgBr3 doped with Mn 2+ ions 
(s = 5/2) were studied by EPR measurements at Q- 
and X-band frequencies with t he magnetic field paralle l 
and perpendicular to the c axis ([McPherson et all Il974t ). 
The EPR spectrum displayed in Fig. [5] shows the hyper- 
fine and fine structures expected for Mn 2+ ions in an ax- 
ial environment. The Q-band frequency of v = 35 GHz 
produces five resonances A to E whenever the spacing 
of adjacent Zeeman-split energy levels corresponds to 
AE = 0.145 meV [see Eq. (j32|)] as illustrated in Fig. H 
Each resonance is characterized by six oscillations due 
to the hyperfine interaction, since the nuclear spin of 
manganese is / = 5/2. The positions of the five reso- 
nances do not occur at equidistant spacings, which in- 
dicates the presence of a non-zero single-ion anisotropy. 
From the positions of the resonances the Lande splitting 
factor g = 2.004(1) and the axial anisotropy parameter 
\D\ = 0.0115(2) meV were obtained. Note that the sign 
of D cannot be determined from EPR experiments at 
elevated temperatures. 



4. Optical spectroscopies 

Single crystals of CsMn^Mgi^Bra (0.04 < x < 
0.20) were investigated by optical spectroscopies 



([McCarthy and Glide! Il984l ). In particular, Mn 2+ pair 
excitations were observed in the absorption spectra as 
shown in Fig. [71 The weak absorptions at T = 1.4 K can- 
not be assigned with certainty; they are either single-ion 
absorptions or due to Mn 2+ clusters with N > 3. With 
increasing temperature additional bands appear due to 
the successive population of the cluster states S = 1 to 
S = 4. The temperature dependence of the intensities is 
best described by using the HDVV Hamiltonian Eq. ([45]) 
with J = —0.88 meV as illustrated in Fig. [71 



5. Inelastic Neutron Scattering 

INS experiments performed on a single crystal of 
CsMno.28Mgo.72Br3 gave evidence for well defined Mn 2+ 
dimer transitions as shown in Fig. [8] ( Falk et ~al\ , |l984). 
The observed intensities are in excellent agreement with 
the predictions by the structure factor Eq. ([53]) ; with 
R = (0,0,1/2) the intensity has a maximum for Q = 
(0,0,1) and vanishes for Q = (0,0,2). The energies 
of the transitions |0) H> |1), |1) H> |2), |2) H> |3), and 
|3) |4) turned out to be 1.80(1), 3.60(1), 5.27(2), and 
6.74(3) meV, which considerably deviate from the Lande 
interval rule, so that the data analysis was based on 
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FIG. 7 (Color online) (a) Temperature dependence of absorp- 
tion spectra resulting from a spin-flip process observed for 
CsMno.2oMgo.8oBr3. The successive population of the Mn 2+ 
dimer states is marked by S = 1 to S = 4. (b) Observed 
intensities of the bands S = 1 to 5 = 3 as a function of tem- 
perature. The lines correspond to the Boltzmann populations 
calculated from Eq . (E E S]) wit h J = —0.88 meV . Adapted from 
iMcCarthv and GudeTll984l . 
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FIG. 8 (a) Energy spectra of neutrons scattered from Mn 2+ 
pairs in CsMn .28Mg .72Br 3 at T = 30 K. (b) Energy- 
level sequence of an antiferrom agneticalry coupled spin pair. 
Adapted from Falk et all Il984l . 



Eq. (jUJ). The resulting parameters are J = —838(5) fieV 
and K = 8.8(8) /ieV. 

Later INS experiments gave evidence f or well de- 
fined M n 3+ trimer and tetramer transitions ([Falk et all 
Il987ah . For the evaluation of the differential neutron 
cross section we refer to references in Sec. III.B.li For a 
trimer the selection rules of the transition \S13SM) —> 



TABLE II Resulting parameters from least-squares fits to 
the observed Mn trimer transitions in CsMno.2sMgo.72Br3 
(jFalk et q/.Lll98^ . 



Model 


J [fieV] 


J' [fieV] K [fieV] L [fieV] 


x 2 


a (K = L = 0) 


-870(12) 


-8(14) 


5.22 


b (K^0,L = 0) 


-786(7) 


-12(10) 14(2) 


3.13 


c (K^0,L^0) 


-777(6) 


-11(9) 8(1) 6(1) 


1.67 



\S' 13 S'M') are derived as 

AS = 0,±1, 

AS 13 = 0,±1, 

AM = 0,±1. 



(59) 
(60) 
(61) 



The smallest magnetic systems to identify three-spin 
interactions are spin trimers. The bilinear Hamiltonian 
Eq. (|56j) has to be extended in the following way: 



H 



-2 J (§1 • s 2 + s 2 • §3) - 2 J'si • s 3 



—K 



(si • s 2 ) 2 + (s 2 • S3) 2 - K' (si • S3) 



-L [(si • s 2 ) (s 2 • s 3 ) + (s 3 • s 2 ) (s 2 • si)] . (62) 

K, K', and L denote biquadratic two-spin and three- 
spin exchange parameters, respectively, which give rise 
to off-diagonal matrix elements, so that Eq. ([62]) was di- 
agonalized in first-order perturbation theory. The bi- 
quadratic K' term is neglected, since \K'\ <C \K\. The 
low-energy part of the eigenvalues E(Sis,S) is illus- 
trated in Fig. [9] for the case of Mn trimer excitations 
in CsMno. 2 8Mgo.7 2 Br3, which were identified in INS ex- 
periments according to the characteristic depe ndence of 
the c ross section Eq. (fT6|) upon Q and T (|Falk et all 
1 19861 ) . The observed transitions are marked by arrows in 
Fig. [9j Least-squares fits based on Eq. ([62]) with differ- 
ent parameter selection gave the results listed in Table ITfl 
The model including only bilinear exchange interactions 
failed as expected. The model including the bilinear and 
biquadratic terms of the Hamiltonian Eq. ([62]) resulted 
in an improved standard deviation % 2 , but only the least- 
squares fit including the three-spin interaction was able 
to reproduce the observed transitions satisfactorily. 

Recent INS experiments performed with increased 
instrumental energy resolution gave evidence for 
anisotropy-induced s plittings of Mn 2+ dimer and 
tetramer transitions (Fu rrer et aD, l2011al) . This is 
demonstrated in Fig. [10] for the dimer \S = 0) — >• \S = 1) 
transition. There are two well defined lines A and B 
which according to the approximate intensity ratio 2:1 
can be attributed to the \S = 0, M = 0) \S = 1, M = 
±1} and |5 = 0,M = 0) \S = 1,M = 0) transitions, 
respectively. A similar anisotropy-induced splitting was 
observed for the lowest tetramer |5 = 0) —> \S = 1) 
transition as well. The dimer and tetramer data could 
be rationalized by the combined action of a single-ion 
anisotropy parameter D = 0.0183(16) meV defined by 
Eq. ([TO]) and of two-ion anisotropic coupling parameters 
J = Jxx = Jyy = -0.852(3) meV and J zz /J = 0.997(1) 
defined by Eq. ([3]). The two- ion anisotropy is most likely 
due to the anisotropic part of the dipole-dipole inter- 
action Eq. ([4]). The exchange coupling J is sufficiently 
strong to keep the spins s^ antiferromagnetically aligned 
at low temperatures T < | J|/&b, but their direction with 
respect to R||c is free to rotate. Therefore, the second 
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FIG. 9 Energy level splittings of Mn trimers in 
CsMno.28Mgo.72Br3. The arrows deno te the observed tran- 
sitions. Adapted from iFalk et all Il986l . 



term of Eq. (jlj has to be averaged in space: 



(63) 



The dipole-dipole anisotropy calculated from Eq. ([63]) is 
Jzz/ J = 0.997, in agreement with the experimental find- 
ings. 



6. Comparison of different experimental techniques 

Table IIIII lists the results obtained by different exper- 
imental techniques presented in the preceding subsec- 
tions. The sign of the axial anisotropy parameter D could 
unambiguously be determined by neutron spectroscopy, 
in contrast to the EPR experiments. The parameters 
D and J exhibit a pronounced temperature dependence 
probably due to the lattice expansion with increasing T, 
whereas the parameter K remains constant. Since the 
analysis of the optical data was based on a model with 
K = 0, the resulting exchange parameter J cannot be 
compared wit h the resu l ts of the INS exper i ments . It 
was shown in Fal k et all Il984t IStrassle et aR [2004 that 
the presence of biquadratic exchange (K ^ 0) is caused 
to a ma jor ex tent by the mechanism of exchange striction 
(|Kittei1l960h . 

In extended antiferromagnets the observation of the 
spin-wave dispersion by single-crystal INS experiments 
is usually the most common approach to determine ex- 
change parameters. By applying the spin- wave formalism 
to H from Eq. ([62]) , which includes higher-order exchange 
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FIG. 10 (Color online) Energy spectra of neutrons scattered 
from CsMn^Mgi-ccBrs. The energy resolution amounts to 
55 fieV. For clarity, the data for x = 0.14 and 0.28 are shifted 
by 10 and 20 intensity units, respectively. The lines refer to 
Gaussian peak fits with equal line widths for both transitions 
A and B. Adapted from iFurrer et all l2uTIal . 



TABLE III Axial anisotropy parameter D and spin coupling 
parameters J and K of Mn 2+ dimers in CsMn x Mgi_ x Br3 
determined by different experimental techniques. 

Technique T [K] D [meV] J [meV] K [meV] 



EPR a 
Optical 6 
INS C 
INS d 
INS e 



77 
13 
30 
50 
1.5 



±0.0115(2) 



0.0183(16) 



-0.88 

-0.838(5) 0.0088(8) 

-0.823(1) 0.0087(2) 

-0.852(3) 0.0086(2) 



a Ref. [McPherson et al, lQZj 
fe Ref. McCarthy and Gudel, 1984. 
c Ref. Falk et a/., 198J. 
d Ref. Strassle et al, 2004. 
e Ref. Furrer et al, 2011a. 



terms, we find 
hu(q) - 

J eff ~- 



4s\J eff \sm(qc), 



(64) 



bs 



VJ(J-4J') + y(# + 2L)|, (65) 



where q is the wave number of the spin wave propagating 
along the c axis. From spin- wave experiments performed 
for the one-dimensional antiferromagnet CsMnBr3 (s = 
5/2) the excha nge coupling wa s determined to be J P jf = 
-0.89(2) meV ([Breitling et all 1 19771 : [Falk et aZ.LIl987bh . 
The analysis of the spin- wave dispersion just yields an ef- 
fective exchange parameter J e f f , but the individual sizes 
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FIG. 11 (Color online) High-field differential magnetization of 
{Ni4Moi2}. The circles represent the experimental data taken 
at 0.44 K, and the line a calculation based on the Hamiltoni- 
ans Eqs. (66}, ([TO]) and (fTTI) with the model par ameters listed 
in the text. Adapted from lSchnack et all . 2006. 



FIG. 12 (Color online) Energy levels of S = 1 tetramers cal- 
culated from Eq. (|67|) with J < 0. The lowest states of a given 
S value are marked by bold lines. The \S34S234S) states are 
identified at the right hand side. The insert shows the cou- 
pling parameters in a slightly distorted tetrahedron as realized 
for the compound {NuMoj^}. 



of the bilinear and biquadratic exchange parameters can- 
not be determined. This is in contrast to experiments 
on small magnetic clusters as discussed in the preceding 
sections. A numerical comparison of the J e ff values ob- 
tained from the three models listed in Table ITT1 is interest- 
ing. Using Eq. (j62j) and s = 5/2 we find J eff = -0.89(2), 
-0.90(3) and -0.92(3) meV for models a, b and c, re- 
spectively. The three values are identical within exper- 
imental error, and they excellently agree with J e ff de- 
termined from spin- wave experiments. They also agree 
with J e jf = —0.88 meV derived from the optical spec- 
troscopies applied to the Mn 2+ dimer excitations, see Ta- 
ble IIIII J e f f is obviously independent of the geometric 
size of the coupled magnetic ions and can therefore be 
regarded as a measure of the magnetic energy per Mn 2+ 
ion in the AFM state of CsMnBr 3 . 



C. Further insights from magnetic cluster excitations 

1. Exchange parameters from high-field magnetization steps 

Step-like features in high-field magnetization data re- 
sult from level crossings associated with the ground 
state of magnetic clusters and thereby provide in- 
formation about the exchange parameters. This 
is demonstrated here for the tetrameric nickel com- 
pound [Moi 2 028(OH)i 2 {Ni(H 2 0) 3 }4]-13H 2 0, henceforth 
abbreviated as {Ni4Moi 2 }. The magnetization is en- 
hanced from zero up to the saturation value of 8 \±b in 
steps of 2 hb at the fields 4.5, 8.9, 20.1, and 32 T as illus- 
trated by the differen tial magnetization data in Fig. fTTI 
(|Schnack et a/J . l2006f ). 

The four antiferromagnetically coupled Ni 2+ ions (s — 
1) in {Ni4Moi 2 } are arranged in a slightly distorted tetra- 
hedron, i.e., the Ni(2)-Ni(3) and Ni(2)-Ni(4) distances 
are slightly shorter than the other four Ni-Ni distances 
as shown in the insert of Fig. [12j thus the spin Hamilto- 



nian is described by 

H = -2 J (§1 • s 2 + §1 • s 3 + §1 • s 4 + s 3 • s 4 ) 

-2 J' (s 2 • s 3 + s 2 • s 4 ) , (66) 

which can be brought to diagonal form by choosing the 
spin quantum numbers according to the vector couplings 
S34 = s 3 + s 4 , S 234 = s 2 + S34, and S = §1 + S 234 
with < £34 < 25, |S 3 4 - s| < £234 < (#34 + s), and 
|5 f 234 — s\ < S < (S234 + S), respectively. The eigenvalues 
of Eq. [66] for si = s 2 = 83 = 54 = s are then given by 

E(Ss4 J 5234? S) = — J[S(S + 1) — 5234(5234 + 1) 

— 534(534 + 1)] — J' [5234 (5234 + 1) 

-5 3 4(5 34 + 1)- 5(5 + 1)]. (67) 

Figure [12] displays the energy levels £^(534,5234,5) 
normalized to J (assuming AFM coupling J < 0) as a 
function of the ratio x = J' / J . For x = 1 the energy 
levels are degenerate with respect to the total spin 5, 
and the energy splittings follow the Lande rule Eq. (|46|) . 
By applying a magnetic field the ground state changes 
in steps from 5 = to 5 = 4 for the field values corre- 
sponding to the maxima of the dM/dH data displayed 
in Fig. HU 

INS spectra measured for a polycrystalline sample of 
{M4M012} are shown in Fig.[l3l from which two ground- 
state transitions can be identified at ca. 0.5 and 1.7 meV 
(jNehrkorn et all [2010). The former is composed of two 
subbands at 0.4 and 0.6 meV attributed to a ZFS from 
magnetic anisotropy. An excited-state transition appears 
at 1.2 meV. From Fig. [12] we can readily conclude that the 
singlet 1 2, 1, 0) has to be the ground state. Moreover, the 
first-excited state has to be the triplet |2, 1, 1) centered 
at 0.5 meV, since transitions between 5 = states are 
not allowed. The observed splitting of the triplet |2, 1, 1) 
into two components can be ascribed to an axial single- 
ion anisotropy defined by Eq. ([TO]) , which has the effect to 
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FIG. 13 (Color online) Energy spectra of neutrons scattered 
from polycrystalline {N14M012} with an incident neutron en- 
ergy of 3.27 meV. The insert attributes the observed transi- 
tions to_th£j£w-eji^rgy_D_art of the splitting pattern. Adapted 
from iNehrkorn et a/1 120101 . 



split the states | £3462346') into the states 1 634 6234 SM). 
A least-squares fit to the energy spectra of Fig. [13] on the 
basis of the Hamiltonians Eq. ([66]) and ([TO]) converged to 
the parameters J = -0.25(2) meV, J' = -0.53(4) meV, 
D 0.22(5) meV ([Furrer et all l2010h which nicely re- 
produce the high-field magnetization data, see Fig. [TTJ 
The resulting low-energy splitting pattern is sketched in 
the insert of Fig. [13] 



2. Pressure dependence of exchange parameters 

By using external pressure the exchange parameters J 
can be determined for varying distance R between the 
magnetic ions. This is of importance for testing and 
improving theoretical models of the exchange interac- 
tion, where the distance usually enters in a straightfor- 
ward manner. A detailed understanding of the exchange 
interaction is, for instance, indispensable for the engi- 
neering of spintronics devices made of magnetic semi- 
conductors. In an effort to shed light on this issue 
the pressure dependence of J was investigated for an- 
tiferromagnetically coupled Mn 2+ dime rs in the semi- 
cond ucting compound Mno.02Zno.98Te ([Kolesnik et all 
120061 ). The corresponding energy- level scheme is indi- 
cated in Fig. [8] INS experiments performed for pressures 
of p = and p = 0.4 MPa gave evidence for an apprecia- 
ble pressure-induced upward shift of the observed dimer 
excitations as illustrated for the |2) — |1) transition with 
energy A\J\ in Fig. [14] The pressure-induced change of 
J amounts to dJ — —0.040(9) meV, accompanied by a 
0.49% decrease of the intra-dimer distance R, resulting 
in a linear distance dependence \dJ/dR\ = 1.8(4) meV/A 
for dR <C R. Similar INS experiments performed for 
CsMno.28Mgo.72Br3 gave evidence for a much stronger 
distance dependence of J with \dJ/dR\ = 3.6(3) meV/A 
(jStrassle et all 12004] ). 
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FIG. 14 (Color online) Pressure dependence of the |2) — )> 
|1) transition associated with Mn 2+ dimers in Mno.02Zno.9sTe 
measured by I NS. The lines denote Gaussian fits to the data. 
Adapted from iKolesnik et all 120061 . 



3. Doping dependence of exchange parameters 

It has been shown that a magnetic semiconductor can 
be converted by hol e doping from i ts intr insic AFM state 
to a ferromagnet ([Ferrand et all l200lf ). It is still an 
open question whether the holes are localized or itin- 
erant. The doping dependence of the exchange param- 
eters may shed light on this issue. Neutron spectro- 
scopic measurements were performed for single crystals 
of Mn^Zni-^Te, one with x = 0.05 and doped with P to 
a level of 5xl0 19 cm~ 3 , and another und oped reference 
sample with x = 0.02 ( Kepa et~al\ , l2003f ). The experi- 
ments were similar as those described in Sec. IIII.C.2l and 
gave evidence for a distinct doping-induced downward 
shift of the observed Mn 2+ dimer excitations as illus- 
trated for the 1 2) — >> |1) transition with energy 4\J\ in 
Fig. [T5] The doping-induced change of the exchange en- 
ergy J amounts to dJ = 0.013(3) meV, in reasonable 
agreement with dJ = 0.010 meV calculated from the 
RKKY interaction which indicates that the ferromag- 
netic exchange is mediated by weakly localized holes. 



4. Anisotropic exchange interactions 

Exchange anisotropy can generally be expected for 
type-L and type-J compounds where orbital degener- 
acy is present. This is exemplified here for the type- 
L compound KioICo^DsO^PWgOs^l^ODsO which 
contains a tetrameric Co 2+ cluster as sketched in Fig. [16] 
The combined action of spin-orbit and crystal-field inter- 
actions splits the 4 Ti single-ion ground state of the Co 2+ 
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FIG. 15 (Color online) Doping dependence of the |2) — >> |1) 
transition associated with Mn 2+ dimers in Mn^Zni-^Te mea- 
sured by INS. The lines denot e Gaussian fits to the data. 
Adapted from Kep a et al. 1 12001 
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FIG. 16 (Color online) Energy spectrum of neutrons scat- 
tered from Kio[Co4(D20)2(PW9 034)2]-20D 2 at T = 1.7 K. 
The inserts show a sketch of the tetrameric Co 2+ unit and 
the resulting energy level scheme, labeled with the dominant 
component \S12S34S, ±M) of the wave function. The highest 
states |1, 1, 1, 0) and |1, 1, 0, 0) located at 5.0 4 and 5.79 meV 
respectively, are not shown. Adapted from lOemente et all 
119971 . 



ions into six anisotropic Kramers doublets ( Carlin[[l986l ). 
By considering only the lowest single-ion level, an effec- 
tive spin Hamiltonian of the Co 2+ tetramer with 8 = 1/2 
for all ions can be written as 

H = — 2 J[si yX Ss yX + Si^ x S^ x + S2 yX S3,x + S2,x$4,x 
,y H~ s l,y s 4,y $2,y s 3,y s 2,y s 4,y 
+Ot(si,zh,z + h,zh,z + h.zh.z + S 2 ,zh,z)] 
-2 J f (si,xS2,x + h,yS2,y + OL's^zh^z)- (68) 

It turns out that for this particular system the eigen 
functions are approximately given by the spin func- 
tions | Si 2 S34 SM) constructed through the spin coupling 
scheme S i2 = §i + s 2 , S 34 = s 3 + s 4 , S = S12 + S 34 , with 
less than 1% S- mixing. 

INS experiments gave evidence for well defined tran- 
sitions associ ated with the Co 2+ tetramer as shown 
in Fig. d6] ([Clemente et all Il997h . The data analy- 
sis based on Eq. ([68]) provided the exchange parame- 
ters J 0.52 meV, a 2.4, J' 0.11 meV, and 
a' = 4.6, resulting in the energy level scheme indicated 
in Fig. [161 Both interactions J and J' are ferromag- 
netic, thus leading to an M = ±2 ground state, in agree- 
ment with magnetic suscept ibility and EPR experiments 
([Gomez-Garcia et all Il992[ ) . The exchange anisotropy is 
rather large as expected from the anisotropy of the Lande 
g matrix with components in the range 2.6-7.0 observed 
by EPR experiments. 



5. Higher-order single-ion anisotropics 

Anisotropy-induced ground-state level splittings are 
essential to understand the step-like magnetic hystere- 
sis curves and the related relaxation and spin reversal 
phenomena observed in the single-molecule magnets (see 



TABLE IV Anisotropy p arameters de termined for Fe4 
isomer AA) by EP R ((Bouwen et all l200lf ) and INS 
Amoretti et all l200ll ) experiments. 



D [/xeV] E [/xeV] B° A [jieV] B\ [/xeV] B\ [fieV] 

. _ 10 (4)xlO- 2 -0.(4)xl0 -4 

INS -25.3(2) 



-2.5(2) -1.5(3)xl0" 



Sec [Vj). EPR is the experimental tool of choice to de- 
termine ground- state level splittings, but because of the 
typically large ZFS parameters D high magnetic fields 
and/or high frequencies are needed to obtain sufficiently 
resolved spectra. INS experiments offer a valuable al- 
ternative in zero field which will be demonstrated here 
for the tetrameric iron compound Fe4(OCH 3 )6(dpm)6, 
or Fe4 in short, which has an S = 5 ground state and 
shows slow relaxa tion of the magnetization below 1 K 
( Barra et all Il999h . The four iron atoms lie exactly in 
a plane, with the inner Fe atom being i n the center of 
an is osceles triangle. EPR experiments ( Bouwen et all 
l200lh have shown that the single-ion anisotropics defined 
by Eqs. ([TP]) and (fTTj) are not sufficient to reproduce 
the observed signals in the ground-state multiplet, but 
higher-order anisotropy terms as given in Eq. ([T2]) are 
needed. The anisotropy parameters resulting from the 
EPR experiments are listed in Table IIVI 

It can be seen from Table IIVI that D is the dom- 
inant anisotropy parameter which splits the S = 5 
ground state into a sequence of five doublets (M = 
±5,±4,±3,±2,±1) and a singlet (M = 0) with ener- 
gies DM 2 according to Eq. (|T0)) . The other anisotropy 
parameters produce a slight mixing of the \M) states 
and/or give rise to small energy shifts. Nevertheless, the 
relevant selection rule in INS experiments is retained with 
good accuracy as AM = dzl, as demonstrated by the 
data displayed in Fig. [13 ([Amoretti et all l200lh . The 
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FIG. 17 (Color online) Energy spectrum of neutrons scattered 
from Fe 4 at T = 6.5 K . The transitions | ± M) -> \ ± M') 
within the 5 = 5 ground-stat e multiplet are marked for each 
observed peak. Adapted from lAmoretti et aZ.L l200ll . 
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FIG. 19 (Color online) (a) Magnetic field dependence (5||c) 
of the singlet-triplet excitations observed in SrCu2(BOs)2 by 
EPR (open dots) and INS (solid dots) experiments, (b) INS 
spectrum of Sr Cu 2 (BQ 3 )2 at Q = (1, 0,0) and B = QT ( B\\c). 
Adapted from iNoiiri et all Il999l and ICepas et all 120011 . 




FIG. 18 (Color online) Structural arrangement of the Cu 2+ 
dimers (red balls connected with black lines) in SrCu2(BOs)2 
including the Dzyaloshinski-Moriya interactions whose vec- 
tors are perpendicular to the (a,6)-plane. The arrows show 
the order of the spins in the expression d(si x Sj). Adapted 
from ICepas et all l200ll . 



decreasing energy spacing with decreasing energy trans- 
fer prevented the transition \M = ±1) — » \M = 0) to be 
resolved from the elastic line. The anisotropy parameters 
resulting from the INS experiments are listed in Table HV1 
The value of E is twice that obtained by the EPR mea- 
surements. Small discrepancies between the parameters 
determined by INS and EPR are frequently observed and 
are probably due to the high magnetic fields used by the 
latter technique, producing a rather large Zeeman split- 
ting so that a mixing of the ground and the excited spin 
multiplets cannot be neglected. 



dicular Cu 2+ dimers (s = 1/2) as illustrated in Fig. [181 
The gap associated with the singlet-triplet dimer excita- 
tions was determined by INS experiments to be 3.0 meV 
( Kagevama et all l2QQ0h . An almost perfect center of 
inversion at the middle of the dimer bonds forbids the 
Dzyaloshinski-Moriya (DM) interaction [Eq. J5J)] between 
the two spins of a dimer. However, each dimer is sepa- 
rated from the neighboring dimer by a BO3 unit, for 
which there is no center of inversion, so that the DM 
interaction is allowed between the next-nearest-neighbor 
(nnn) Cu 2+ spins as described by the Hamiltonian 



H = ^ ±de c ■ (si x Sj 



(69) 



where the sign depends on the bond (see Fig.[T8l) : and e c 
is the unit vector in the c direction ( Cepas et a/.Ll200lh . 
The effect of Eq. ((69]) is to split the AM = ±1 tran- 
sition associated with the singlet-triplet splitting into 
four branches under the application of a magnetic field 
parallel to the c axi s. This pred i ction was verified by 



EPR measurements (INoiiri et all 11999 ) and later con- 
firmed by INS experiments (|Cepas et a/.j200lf ) as shown 
in Fig. [T9l In addition to the four AM = ±1 branches, 
the field-independent AM = transition was observed 
in the INS measurements. The DM parameter resulting 
from these investigations turned out to be d ~ 0.18 meV, 
which roughly compares with the estimated value d = 
(Ag/g)J nnn w 0.5 meV. 



7. A novel tool for local structure determination 



6. Dzyaloshinski-Moriya interactions 

The compound SrCu2(BOs)2 is a two-dimensional 
spin-gap system with tetragonal unit cell. It consists of 
alternately stacked Sr and CUBO3 planes. The latter 
are characterized by a regular array of mutually perpen- 



Conventional crystallography is the standard tool for 
structure determination, and a periodic lattice is a pre- 
requisite for such studies. However, complex materi- 
als are often characterized by local deviations from per- 
fect periodicity which may be crucial to their properties. 
The most prominent bulk methods for local structure de- 
termination are x-ray absorption fine structure, nuclear 
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FIG. 20 (Color online) (a) Energy distribution of the \S = 
0,M = 0)-^|5 = 1,M = 1) Mn 2+ dimer transition observed 
for CsMno.1Mgo.9Br3 at T = 1.6 K. (b) Sketch of statistical 
distributions on Mn 2+ and Mg 2+ ions along the c axis, m 
is the number of Mn 2+ ions replacing Mg 2+ ions outside the 
centr al Mn 2+ dimer emb edded in the brown area. Adapted 
from iFurrer et all l2011bl . 



TABLE V Analysis of the INS spectrum of CsMno.1Mgo.9Br3 
given in Fig. l20l a). E m and denote the energy transfers 
and the normalized intensities of the individual bands m, re- 
spectively. The intensities I^ lc were calculated acc ording to 
a statistical model described in Fu rrer et al. I l2011bl . Relative 
error bars are given for the exchange couplings J m and for 
the local Mn-Mn distances R m - 
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1.762(2) 


0.081(22) 


0.119 


-0.8321(10) 


3.2311(3) 
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1.780(2) 


0.302(25) 


0.289 


-0.8408(10) 
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1.794(2) 


0.383(22) 
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-0.8477(10) 


3.2268(3) 
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1.811(2) 


0.173(22) 


0.193 


-0.8563(10) 


3.2244(3) 


4 


1.826(3) 


0.061(18) 


0.089 


-0.8637(14) 


3.2224(4) 



magnetic resonance, and atomic pair-distribution func- 
tion analysis. All these methods provide a spatial res- 
olution of typically 0.1 A, and their performance can 
hardly be improved. Magnetic cluster excitations are 
able to push the spatial resolution beyond the present 
limits through the dependence of the exchange coupling 
J on the interatomic distance R, which for most mate- 
rials is governed by the linear law dJ/dR — a as long 
as dR <C R. Modern spectroscopies measure exchange 
couplings with a precision of dJ/J « 0.01, thus spatial 
resolutions of dR ~ 0.01 A can be achieved as demon- 
strated below. 



We turn to Fig. [10] which displays Mn 2+ dimer ex- 
citations in CsMn^Mgi-^Bra. The total linewidths of 
the transitions A and B show an x-dependent increase 
beyond the instrumental energy resolution (FWHM 
= 55 /xeV) which was further investigated by high- 
resolution INS expe riments with FWHM = 15 fieV 
([Furrer et all . l2011bh . As shown in Fig. IMaV the en- 
ergy spectrum observed for the transition A exhibits 
marked deviations from a normal Gaussian distribution. 
It is best described by the superposition of five indi- 
vidual bands which correspond to specific Mn 2+ dimer 
configurations with particular exchange couplings J m 
(m = 0,1,2,3,4). The linear law dJ/dR = a was es- 
tablished for CsMn^ Mgi-^Bra with a = 3.6 meV/A 
([Strassle et all I2QQ4J ) . thus each of the five J m values 
can be associated with a particular local distance R m as 
listed in Table El 

How can the discrete nature of the local Mn-Mn dis- 
tances be explained? In Fig. l20Tb) different configura- 
tions along the Mn chain structure are sketched, where 
m is the number of peripheral Mn 2+ ions replacing the 
Mg 2+ ions. The introduction of additional Mn 2+ ions 
exerts some internal pressure within the chain, since the 
ionic radii of the Mn 2+ (high spin) and Mg 2 + ions are 
different with rMn =0.83 A> rug = 0.72 A (j Shannon! , 
Il976l h so that the atomic positions have to rearrange. 
In particular, the Mn-Mn bond distance of the central 
Mn 2+ dimer is expected to shorten gradually with in- 
creasing number of Mn 2+ ions as compared to the case 
m = 0. For any number m there is a myriad of structural 
configurations, resulting in a continuous distribution of 
local distances R m . This view, however, is in contrast 
to the observed energy spectrum displayed in Fig. l20Ta) 
which is clearly not continuous. In other words, the bond 
distance is not smoothly adjusted to its surrounding but 
locks in at a few specific values Rm. Obviously the re- 
alization of discrete local distances R m is governed by 
the number m of peripheral Mn 2+ ions and not by their 
specific arrangement in the chain. This surprising result 
is due to the one-dimensional character of the compound 
CsMn^Mgi-^Bra in which the mixed Mn T M gi_, T chains 
beha ve like a system of hard core particles ([Krivoglazl 
Il996h , In conclusion, the use of high-resolution spectro- 
scopies allows a rather direct determination of local inter- 
atomic distances in small magnetic clusters, in contrast 
to other techniques which usually have to be combined 
with simulations. 



IV. LARGE MAGNETIC CLUSTERS 
A. Introduction 

In the previous section, the scientific questions ad- 
dressed by the small magnetic clusters focused on demon- 
strating and elucidating the nature of the various basic 
interactions between the spin centers in condensed mat- 
ter systems. However, in large magnetic clusters, or the 
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molecular nanomagnets in the context of this review, the 
huge size of the Hilbert space makes a complete exper- 
imental characterization of the magnetic cluster excita- 
tions (usually) impossible. Accordingly, the fine details 
in the spin interactions such as exchange anisotropy are 
not detected, and the modeling of the data can with much 
success be based on spin models taking into account only 
the dominant interactions, which in most cases is Heisen- 
berg exchange. The spectroscopic experiments typically 
reveal the low-lying excitations or the low-energy sec- 
tor of these spin models, and one is hence naturally di- 
rected towards questions such as what is the nature of 
the ground state and elementary excitations. 

The key distinguishing feature as compared to the pre- 
vious section is the many-body structure of the wave- 
functions in the large magnetic clusters, and a main goal 
could be formulated as what novel quantum states are 
realized and which physical concepts allow us to ratio- 
nalize them. At this point the close ties to the field of, 
e.g., quantum spin systems becomes apparent, and in- 
deed methodologies developed there are often applied to 
molecular nanomagnets. Most of the examples presented 
in this section will elaborate on that. 

However, also distinguishing novel aspects come into 
play as a consequence of the fact that the molecular nano- 
magnets are not extended. For instance, phase transi- 
tions, either classical or quantum, are not possible in 
a strict sense. Conceptually most important, however, 
is that the wave vector q ceases to be a useful quan- 
tum number. One can actually expect that exactly those 
lattice topologies which cannot be expanded into an ex- 
tended lattice will exhibit the most interesting novel com- 
plex quantum states and magnetic phenomena. Research 
into this direction has just started, however, and only 
preliminary results are available at the time of the writ- 
ing of this manuscript. A further important point not 
emphasized yet is that the spins of the magnetic cen- 
ters in molecular nanomagnets are generally rather large, 
with Si = 3/2,2,5/2 being most often found, in con- 
trast to quantum spin system, where much focus is on 
spin- 1/2 systems. Typical metal ions would be Cr 3+ 
and Mn 4+ (s» = 3/2), Mn 3+ (s* = 2), and Mn 2+ and 
Fe 3+ (si = 5/2). The quantum-classical correspondence 
hence enters naturally in the discussion of the magnetic 
excitations in molecular nanomagnets. 

In the following those classes of molecules will be 
discussed for which considerable insight into the mag- 
netic cluster excitations has been obtai ned. Impor- 
tant classes such as odd-membered wheels dCador et all 
\2004 iHoshino et all l2009bl: lYao et all l200^ ferromag ; 
netic clusters (|C lenient e - Juan et all Il999t iLow et al 



20061: IStuiber et all 2011), discs ([Hoshino et all 12009a : 



Koizumi et al. 



20071 ). and others are not mentioned. 
An important subclass of molecular nanomagnets are 
the single- molecule magnets (SMMs), which have re- 
ceived the largest attention in the past and could be 
described as having given birth to the molecular nano- 
magnets as a research field. The above questions are also 



relevant, but the most interesting phenomena in SMMs, 
such as quantum tunneling of magnetization, are mainly 
related to magnetic anisotropy. The SMMs will be dis- 
cussed in Sec. El 



B. Theoretical description 

1. Spin Hamiltonian 

As mentioned before, experimental results on large 
magnetic clusters can often very well be reproduced by 
spin Hamiltonians which include only the most domi- 
nant terms. For clusters containing only type S and Q 
metal ions, on which we focus in this review, these are 
the HDVV Hamiltonian Eq. ([2]), the single- ion anisotropy 
Eqs. (Unj) and ([II]), and the Zeeman term Eq. ([H]). How- 
ever, in molecular nanomagnets the site symmetries of 
the individual spin centers in the cluster are very low, if 
they have symmetry at all. Accordingly, the single-ion 
anisotropy and g factors should in general be described 
as tensors. Also, in molecular nanomagnets often differ- 
ent kinds of metal centers are incorporated. This gives 
rise to the spin Hamiltonian 



H = -2 JjjSj • sj + 3j • Dj 

i<j i 
+/iB^Si gi B, 



(70) 



which in the following will be referred to as microscopic 
Hamiltonian in order to clearly distinguish it from effec- 
tive models which also appear. The dipole-dipole inter- 
action Eq. (J2J) has also to be included, but its effect on 
energy spectrum and magnetic behavior is very similar to 
that of the single-ion anisotropy term and may hence be 
lumped into the single- ion parameters. Experimental D 
and E values which were derived with the dipole-dipole 
interaction explicitly included in Eq. ([70]) will be indi- 
cated by a superscript " lig" . 

In contrast to the sites, the molecule itself may exhibit, 
or closely approximate, a high molecular symmetry. In 
fact, clusters with a particular symmetry are appealing 
for physical studies, and are thus preferred objects of 
investigations. The microscopic Hamiltonian simplifies 
then enormously and includes only few parameters. 

Usually the HDVV Hamiltonian dominates over the 
single-ion anisotropy, and the strong-exchange limit 
(Sec. III. Aft is an excellent starting point. Although the 
magnetic anisotropy cannot be ignored in the analysis of 
experimental data, the physics of interest in these cases is 
(usually) related to the Heisenberg interactions, and the 
discussion focuses on the corresponding Heisenberg spin 
models (notable exceptions are the SMMs discussed in 
Sec.[V|). The magnetic anisotropy may however also be so 
large that important effects appear which are not covered 
by the strong-exchange limit (S m ixing) ([Liviotti et all 
2002; IWaldmann and GiidelL 2005) or may need com- 
pletely different physical concepts for their description. 
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The examples selected below will demonstrate this point. 
It is added that from the values of magnetic parameters, 
such as J and D, it is usually not possible to infer a priori 
whether the strong-exchange limit is obeyed or not. The 
ratio D/J is generally small in molecular nanomagnets, 
and which case is realized needs a case-to-case analysis. 

The dimension of the Hilbert spaces encountered in 
large magnetic clusters poses a major obstacle, similar to 
that found in other areas such as quantum spin systems, 
and the same conceptual ideas are followed to tackle it. 
Indeed, essentially any technique developed in the con- 
text of quantum spin systems is also of interest for large 
magnetic clusters. However, some of them had been of 
particular use, and are mentioned next. 



2. Numerical techniques 

A most straight forward approach is to numerically 
solve the spin Hamiltonian for its energies and eigen 
functions, which is achieved in two steps, setting up the 
Hamiltonian matrix and then diagonalizing it. 

The major decision in the first step is the choice of the 
basis set, which could be the product states |{M^}) (with 
an obvious shorthand notation) or the spin functions 
\rjSM), where r] denotes the intermediate spin quan- 
tum numbers generated in a particular spin coupling 
scheme. The product states are most easily handled 
in computer code, yield sparse Hamiltonian matrices, 
and are eigenfunctions of S z which allows a block fac- 
torization for magnetic clusters with uniaxial symmetry. 
On the other hand, for the spin functions it is numer- 
ically demanding to calculate matrix elements (a num- 
ber of Wigner symbols need to be calculated) and the 
matrices are dense, but they have the intrinsic advan- 
tage of being eigenfunctions of the total spin operator S 
which results in a very effective block-factorization for 
the HDVV Hamiltonian. Interestingly, nearly all nu- 
merical work in the field of quantum spin systems has 
been based on product states; spin functions are rarely 
used. However, for large magnetic clusters diagonaliza- 
tion using spin function s has been u s ed wi t h much suc- 
cess dBaran et all l20Q8t iDelfs et ad Il993t iGuidi et ad 
l2QQ4t IWaldmann et ad Il999h . and efficient ITO based 
techniqu es have been develop e d for calculating matrix el- 
ements (IGatteschi and Pardii 1 19931) and employing spa- 



ements (Lxattescni ana rami, wyd) and employing spa- 
tial sy mmetries ([Schnalle and Schnacld . l2010HWaldmanni 

Complete information on the system is obtained by a 
full exact diagonali zation, and seve ral canned computer 
codes are available ( Bai et aZl . l2QQ0h . The largest dimen- 
sion of the Hilbert space which can be handled is ca. 
100 000 on a super computer, or about 15 000 on a (32 
bit) personal computer. If symmetries are systematically 
taken into account, quite large magnetic clusters can be 
treated on personal computers, e.g., a mixed valent Mn- 
[3 x 3] grid molecule with a Hilber t space dimension of 
4 860 000 (jWaldmann et ad l2QQ6ch . 



If full exact diagonalization is not possible, one may 
attempt to obtain the energies and eigen functions in 
a subspace. A first approach is to truncate the space of 
basis function, but the success of the procedure obviously 
depends on how well the selected basis states represent 
the s ought-after eigen functions ([ Schnalle and Schnackl 
l2QQ9h , 

A set of very efficient diagonalization methods is given 
by the sparse matrix diagonalization techniques, which 
allow an exact numerical calculation of a small number 
(~100) of selected ene rgies and /or eige n states, typically 
the low- lying states (jBai et ad 1200 0). In physics the 
most prominently used variant is the Lanczos method, 
while in chemistry the Jacobi-Davidson algorithm is more 
often applied. However, for the specific purpose of calcu- 
lating the low-temperature observables of large magnetic 
clusters, the simpler subspace-iteration techniques turn 
out to be quite powerful, since they provide both energies 
and eigen function even in the presence of degeneracies 
with very comparable convergence rates. 

These techniques employ an iterative process and work 
best for sparse matrices. Because of the latter it is most 
natural to use the product states, though spi n functions 
were applied in few cases ([Guidi et ad l2004h . The iter- 
ative process consists of repeatedly applying the Hamil- 
tonian matrix H to a vector x, 



= (H - <7l) • x„ 



(71) 



where a introduces a shift which allows to optimize the 
convergence rate, and the starting vector x may be a 
random vector. If more than one eigen pair is searched, 
the iteration is applied to a subspace of vectors X. For 
a practical algorithm, som e further improvements are 
suggested (|Bai et ad l2QQ0h . The Lanczos and Jacobi- 
Davidson algorithms are also based on matrix-vector 
multiplications H • x, but employ more sophisticated al- 
gorithms to extract the inf ormation on the eigen pairs 
from the generated vectors (|Bai et al 1 l2QQ0h . 

Besides these approaches a number of numerical tech- 
niques exist which aim at calculating observables directly 
without evaluating eigen pairs explicitly. Among these 
are e.g. Quantum Monte Carlo, Chebyscheff expansion, 
dynamic and finite- temperature Lanczos, transfer matrix 
and (dynamical) density matrix renormalization group 
methods. However, although promising, these methods 
have not yet been applied systematically to the analysis 
of large magnetic clusters as defin ed in this review, but 



few applications were reported (Engelhard t et al, 



Exler and Schnack, 2003; Schnack and Wendland 



2006; 



2010; 



Torbriigge and Schnack, 2007; Ummethum et al, 2012). 



More efforts in this direction would obviously be desir- 
able. 



3. Effective Hamiltonian techniques 

Another approach to describe the relevant low-energy 
excitations in a particular spin model is to replace the mi- 
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croscopic Hamiltonian by an effective spin Hamiltonian 
(mapping), which acts in a Hilbert space of (significantly) 
reduced dimension. It is emphasized that the states in 
the Hilbert space of the effective Hamiltonian do not have 
to be identical to those of the microscopic Hamiltonian. 
An effective Hamiltonian may be constructed from var- 
ious procedures, but the following simple technique had 
been particularly useful for rationalizing the magnetism 
in a number of large magnetic clusters. 

The method can be described as a first-order pertur- 
bation treatment, as integrating out a degree of freedom 
or as a mean-field argument and is guided by physical 
intuition. It starts with combining a subset of the spins 
Si in a collective spin S M = ^Z iefJL §i, where /i stands for 
the set of sites z, and selecting the sector of interest via 
the value of the quantum number associated to S M , 
which is usually the minimal or maximal value. Then it 
holds 



C7.S1 



(72) 



with a projection coefficient q, which depends on Si and 
. For the sector where all spins in the subset \i are fer- 
romagnetically aligned and assumes its maximal value 
S/j, = ^Z^i holds Ci = Si/Sy,. The effective Hamiltonian 
is then finally obtained by inserting Eq. ([72]) in the mi- 
croscopic Hamiltonian, which removes the spins of the 
subset /1 and replaces them by the collective spin S M . 

A typical situation where the above scheme could be 
applied is a cluster with one or few ferromagnetic ex- 
change couplings which are much larger than the other 
ones. Then the spins linked by a strong ferromagnetic 
coupling can be combined into one collective spin which 
replaces them in the spin Hamiltonian. Another situation 
is an antiferromagnetic bipartite lattice, where the spins 
on each sublattices A or B are ferromagnetically aligned 
with respect to each other in the ground state. This sug- 
gests to introduce two sublattice spins and S#, and 
to transform the Hamiltonian according to Eq. ([72]) . The 
subsections below will provide examples for the practical 
application of the scheme. 



4. Condensed matter techniques 

Sophisticated techniques for calculating properties of 
extended (quantum) spin systems have been developed. 
Some of them can directly be applied to large magnetic 
clusters as they work for any number N of spin centers, 
i.e., the thermodynamic limit is not assumed a priori but 
taken after completion of the calculation. However, in 
magnetic clusters translational invariance is not present, 
and it is hence a characteristic feature of the applica- 
tion of these techniques to clusters that they have to be 
adapted to work in real space and not momentum space 
as in extended systems, which may bring in novel as- 
pects. This area is still largely unexplored, but two of 
the simpler techniques in this class of methods have been 
applied with some success. They both elaborate on the 



observation that in molecular nanomagnets the spins Si 
are relatively large. 

Spin- wave theory (SWT), or the set of techniques em- 
braced by this acronym such as linear SWT, interacting 
SWT , or modified and f in ite- siz e SWT (IHirsch and T ang. 
19891: llvanov and Senl I2QQ4J: iRado and SuhF 1l963: 



Takahashil 11987k Zh ong and Sorellai Il993l ). is certainly 



a most successful theory in magnetism. All variants have 
in common that they start from a classical spin config- 
uration in the ground state, and expand around it by 
approximating the single-site spin operators §£. They 
are hence semi-classical in nature, and - technically - 
can be applied to any type of lattice, and in particu- 
lar to large magnetic clusters when formulated in real 
space ([Cepas and ZimanL l2QQ5h . It is necessary to dis- 
tinguish between SWTs for ferromagnetic and antiferro- 
magnetic clusters. For ferromagnetic clusters, SWT al- 
lows us to calculate the zero-temperature excitation spec- 
trum exactly for any arbitrary cluster. However, since 
SWT breaks spin rotational invariance it is conceptually 
inappropriate for antiferromagnetic systems with disor- 
dered ground states, such as one-dimensional (ID) sys- 
tems or magnetic clusters. Nevertheless, it is often found 
to produce signif i cant results e.g. for excitation energies 
(|lvanov and Senl . 120041 ) , but the reliability of the results 
should carefully be checked case-by-case. The fundamen- 
tal questions as regards the applicability of SWT to an- 
tiferromagnetic clusters shall not be addressed in this 
review, but results of some few case-by-case checks will 
be reported. For the technical imp lementation o f SWT 
for magnet i c clu sters we refer to ([Stuiber et al ., 2011; 
IWaldmannl . l2QQ7h . 

The large spins in molecular nanomagnets suggest 
an entirely classical description, where each spin oper- 
ator Si is replaced by a classical vector s$ of appro- 
priate length, or a unit vector times a prefactor 
Si = Si&i. A subtlety arises here as regards the appro- 
priate val ue of si, w hi ch can be argued to be best taken 
as Si = \/sAst 4- 1) (I Joyed. 1 19671 : iLuscombe and Lubani 
Il997h or Si = Si ([Fisherl . 1 19641 ). With this replacement 
the Hamiltonian becomes a classical functional, and the 
ground-state energy and spin configuration is obtained 
by minimization. Also, thermodynamic quantities can 
be calculated for quite large magnetic clusters of hun- 
dreds of sites, and it is here where the classical ap- 
proach was mostly applied to mole c ular n a nomagnets 



proacn was mostly a pplied to molec ular n anomagnets 
(iMuller et all l2001al ISchroder efall 12005a] : lYao et all 
l2006h . The calculation of dynamic quantities is also 
possi ble ([Luscombe and Lubani . Il997t ILuscombe et all 
1 19981 ) , but the discrete energies of the excitations in mag- 
netic cluste rs can of course not be recovered, though their 
intensities ([Mentrup et al. I, l2000h . Hence, the classical 
calculations are not actually suited to reproduce experi- 
mental spectroscopic data, but nevertheless can provide 
significant insight into the physics of the excitations in 
a large magneti c cluster through the quantum-classical 
correspondence ([Schroder et a/ll2005bl ), which is proba- 
bly the most important aspect. 
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C. Even-membered antiferromagnetic molecular wheels 

Molecular wheels are species in which the metal ions 
form almost perfect ring-like structures. The decanuclear 
wheel [Fei (OCH3)2o(0 2 CCH 2 Cl)io], or Fei in short, 
was the first wheel whose magnetic properti es were stud- 
ied ([Gatteschi et a/.l . ll994l : pTaft et all\l994 . Since then, 
because of their aesthetical appeal and unprecedented 
magnetism, the class of the molecular wheels has enor- 
mously grown and dozens of wheels with nuclearity rang- 
ing from N = 6 to N — 18 have been synthesized. 
In Fig. [21] the crystal structures of five representatives 
are displayed. The physics in the molecular wheels, and 
their relatives such as the modified wheels (Sec. IIV.D|) . 
turned out to be surprisingly rich. The presentation here 
is hence necessarily limited. 

The even-membered antiferromagnetic wheels have 
been the subject of intense research. In this section, 
homo-nuclear wheels will be discussed with S{ = s for 
all i. In view of the high molecular symmetry, which in 
hexa-nuclear wheels can be a crystallographic S6 sym- 
metry and in octa-nuclear wheels a C4 symmetry, the 
magnetism should be expected to be very well described 
by the generic spin Hamiltonian 



H 



(N-l ^ 
Si • + sat • §1 
i=i , 



s iz~ 3«t(«» + l) 



■HB9S-B. (73) 



Additional terms describing e.g. a variation of the 
exchange constants along the wheel, tilted single-site 
anisotropy tensors or higher-order terms should in prin- 
ciple also be present, but these are difficult to resolve 
in experiment. Sometimes a planar term E J2i(^ix ~ 
had to be added to Eq. ([73]) (vide infra), and evidence for 
weak Dzyaloshinski-Moriya in t eract ions were reported 
([Cinti et all I2QQ2I : lLante et aZ.L l2QQ9h . However, the lat- 
ter manifests itself only in very specific experiments, and 
is not further considered. The Hamiltonian Eq. ([73]) (plus 
sometimes an E term) generally provides a very good de- 
scription. 

The exchange coupling 2 J is typically on the order 
of few meV, and the single- ion anisotropy weak, \D\ < 
0.1 meV. In most cases it is of the easy- axis type (D < 0), 
with the exception of Feis, where D > (vide in- 
fra). The ratio D/J is small in the antiferromagnetic 
wheels, but the strong-exchange limit is not necessar- 
ily well obeyed. In fact, the effect of the anisotropy 
is rather determined by (Ns) 2 \D/J\ ([Chiolero and Loss! 
1 19981 ; IWaldmannL l2QQ2al ) , and antiferromagnetic wheels 
should be classified according to wether the anisotropy 
is "weak" or "strong" by the different physics in the two 
limits. 

The magnetic susceptibility data have conclusively 
pointed towards a 5 = ground state in the antifer- 
romagnetic wheels, which is intuitively anticipated from 



the expected alternating up and down spin configura- 
tion in the classical ground state. The first general in- 
sight into the excitation spectrum has been obtained from 
low-temperature magnetization curves on the Feio wheel, 
which are reproduced in Fig. [22] A sequence of magne- 
tization steps is seen as function of field, where at each 
step the magnetization changes by 2 \±b , and which occur 
in regular field intervals of ca. 4.2 T. This demonstrates 
that as function of field the ground state changes from 
5 = to 5 = 1 at the first step, 5 = 1 to 5 = 2 at 
the second step, and so on, as sketched in Fig. [23] From 
the fields of the magnetization steps or ground-state level 
crossings one can infer the zero -field energies of the states 
(IShapira and BindilattiL [2002), with the result 



E(S) = -S(S + 1). 



(74) 



The excitation spectrum consists hence of a set of spin 
multiplets with 5 = 0, 1, 2, ... whose energies satisfy the 
Lande rule Eq. ([46]) . see also Fig. [23] where A is the en- 
ergy gap between the singlet and triplet. Such a spectrum 
is also generated by a Heisenberg dimer Hab — AS^-Sb, 
with appropriate spins Sa and Sb- The Hamiltonian 
Hab appears hence as an effective spin Hamiltonian for 
the antiferromagnetic Heisenberg model on a ring. There 
is however some deeper physics in here, which will be 
mentioned in the next subsection. 

The role of magnetic anisotropy in antiferromag- 
netic whee l s was elucidated by EPR dPilawa et all , 
I1997L 120011 l2003t Ivan Slageren et all 120021) and mea- 
surem ent s of the magnetic torque (|Cornia et all 
1999a|IB; IWaldmann et al. , 1999). The first INS ex- 
periment on a wheel was performed on the clus- 
ter Na[Fefi{NfCH 2 CH 2 0),sifilCl, or NaFe 6 (Fig. ED 



([Waldmann et all 
produced in Fig. [24] 
transitions I and II, 



19991 ). The experimental data are re- 



and are characterized by two cold 
and a hot transition III. The in- 
terpretation of the observed transitions is straightfor- 
ward: Peaks I and II correspond to transitions from the 
5 = ground state to the zero-field-split 5 = 1 multiplet, 
and peak III to a transition from the |5 = 1,M = 0) 
level to the 5 = 2 multiplet, as indicated in Fig. [23] 
The INS energies yield directly 2 J = —2.0 meV and 
D = —0.038 meV, in excellent agreement with magnetic 
susceptibility and magnetic torque measurements, which 
puts NaFe6 in the "weak" anisotropy category. 

The following three subsections are organized by the 
"strength" of the magnetic anisotropy or parameter 
(Ns) 2 D/J. First, the physics in the "weak" anisotropy 
case, which is that of the L&^-band picture, is consid- 
ered followed by the "intermediate" anisotropy case. Fi- 
nally, wheels with " strong" anisotropy which may display 
quantum tunneling of the Neel vector are discussed. 
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1. Antiferromagnetic wheels with "weak" anisotropy and the 
L&^-band picture 

Most antiferromagnetic wheels fall into the weak 
anisotropy category. The anisotropy leads to a zero-field 
splitting of the spin multiplets, which is detected in e.g. 
INS experiments as was shown by the example of NaFe6, 
but the physics of the magnetic excitations is not affected 
by it. Hence, the appropriate model for the discussion of 
the physics is that of the antiferromagnetic Heisenberg 
ring (AFMHR) or Eq. (fSJ with D = 0. 

Initially, the wheels were regarded as models for ID 
antiferromagnetic chains with the implication that the 
physical concepts developed for them, which are char- 
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FIG. 22 (a) Magnetization vs field and (b) differential mag- 
netization dM/dB at a temperature of 0.6 K of the antiferro- 
magnetic wheel Feio. The data show the sequence of magne- 
tization steps, which were be tr aced here for fields up to 40 T. 
Adapted from Taft et all 11994 
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FIG. 23 (Color online) Sketch of the low-lying energies in 
even antiferromagnetic molecular wheels. The Heisenberg in- 
teractions give rise to spin multiplets S = 0, 1, 2, . . . which 
split in a magnetic field and lead to level crossings in the 
ground state, which changes from \S = 0,M = 0) to |1, — 1), 
1 1 , — 1 ) — > 1 2, —2), etc., and steps in the magnetization curve. 
The magnetic anisotropy splits the spin multiplets in zero 
field, as indicated to the left. Some allowed INS transitions 
are indicated by arrows. 
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FIG. 24 Neutron energy- loss spectrum of NaFe6 at the indi- 
cated temperatures (Q — 1.0 A -1 ). The transition i could be 
assigned to an impurity species in the sample. Adapted from 
Wald mann et all 1 19991 . 
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acterized by the strong quantum fluctuations in ID, 
should also apply to AFMHRs. However, the experi- 
mental and theoretical work has suggested a very dif- 
ferent picture of the excitations, which is denoted here 
as L&£?-band concept. This concept has its roots 
in Anderson's 1952 paper on antiferromagnetic spin 



(a) 



waves ( Anderson! , Il952fh and emerged in the course 
of several works (lA.ndersoni fl983: iBernu et all 1 1992 : 
Caneschi et al Jl996t IChioTero and LosslJl998l: iLhuilher 



2002: Schnack and Luban 



2000; Waldmann, 2002b) 



The L&£?-band concept is classical in nature, though 
some subtleties are present, and the meaning of "clas- 
sical" in this context is in fact not yet completely under- 
stood. It rather seems that in magnetic clusters of the 
size discussed here th e line between classical and quan- 
tum physics is blurred ([Konstantinidis et all l2Qllh . Fun- 
damental questions hence remain. However, from a prac- 
tical point of view the L&^-band concept seems to apply 
whenever the classical spin structure appears to describe 
well the ground state, and has allowed rationalizing the 
magnetism in a number of different classes of molecules. 
The molecular wheel [CrgFgPivie], or Crg (Fig. [21]) . has 
played an important role in this context, as it was the 
first system for which a detailed, and indeed complete, 
experimental demonstration of the concept was achieved 
by exploiting the full power of INS. 

Descriptions of t he L&£?-band concept are available 
( WaldmarmL l2005al) ; only major aspects are mentioned 
here. The key ingredient is the notion of rotational 
bands, where a rotational band is a sequence of spin 
multiplets with S = Smin,S m in + l,S m in + 2, ... whose 
energies increase according to the Lande rule E(S) oc 
S(S + 1). The analogy of this spectrum with that of 
a rigid rotator H = L 2 /(2J), where L is an angular mo- 
mentum and I the moment of inertia, is not coincidental. 
The L&E'-band concept can be summarized as: 

(1) In an energy- vs-S plot, the low-energy sector is 
characterized as a set of rotational modes; the lowest- 
lying mode called L band and a number tie of higher- 
lying modes called E band. The next higher-lying states 
are collectively denoted as quasi continuum. 

(2) The classification of the energy states in (1) is justi- 
fied by a selection rule: Spin transitions from the L band 
into the quasi continuum are forbidden. Hence, at low 
temperatures only INS transitions between states of the 
L band or from the L to the E band are allowed. 

Important ly, these points have to come together 

(Waldmann, 2002b). The presence of a rotational mode is 
by itself not sufficient to ensure the validity of the LhE- 
band concept and its consequences. In Fig. [25] are shown 
the energy spectrum and spin-spin auto-correlation func- 
tion C§(uj) of a N = 8, s = 3/2 AFMHR. The structure 
of the energies detailed in (1) is clearly seen. Also, the 
correlation function demonstrates, e.g., that transitions 
from the S = ground state into the quasi continuum do 
not occur. The above two points completely define the 
concept, but some further comments are at place: 

L band: The L band is intimately related to the occur- 
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FIG. 25 (Color online) (a) Simulated energy spectrum for 
a N = 8, s = 3/2 AFMHR with 2J -1.46 meV as for 
Crg. The low- lying levels form rotational bands. The L band 
can be rationalized as (quantized) rotation of the Neel vec- 
tor. The next higher-lying bands are denoted collectively as 
E band and may be regarded as (discrete) antiferromagnetic 
spin waves. Some INS transitions are indicated by arrows 
with associated labels. The numbers in brackets indicate 
the shift quantum number q for some spin multiplets. (b) 
Simulated spin-spin auto-correlation function C§(u) in units 
of 2\J\ for various temperatures. The L- and ^-band tran- 
sitions indicated in pa nel (a) a re identified. Adapted from 
Waldmann et al 1 120031 andlWaldmann, 2002& 



rence of long-range Neel order: For lattices which can be 
expanded to infinity, the presence of the L band in the 
spectrum is a necessary but not sufficient requirement 
for an ordered Nee l ground state in the infinite lattice 
(|Bernu et oZl Il992h , ft is hence also related to the cor- 
responding antiferromagnetic Goldstone mode, and the 
L-band states have spatial symmetries consistent with 
th e ordered Neel state. An excellent discussion is found 
in iLhuillierl l2002l Qualitatively, the L band reflects the 
rotational degree of freedom of the classical antiferromag- 
netic spin configuration in the ground state, or (quan- 
tized) rotation of the Neel vector. The L band is also 
known as the tower of states or quasi-degenerate joint 
states. 

E band: The E band is deeply related to spin waves: 
For lattices which can be expanded to infinity, the E band 
evolves into the familiar antiferromagnetic spin waves in 
the infinite lattice (where we deviate from the usual no- 
tation by not regarding the Goldstone mode or antifer- 
romagnetic Bragg peak as a spin wave, see Ref. 25 in 
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FIG. 26 (Color online) Neutron energy- loss spectrum 
of Crs at the indica ted temperatures. Adapted from 
Waldmann et al. I I2003L 



IWaldmarmL l2QQ7h . In finite clusters, the spin- wave spec- 
trum becomes of course discretized. It then consists of 
tie energies, and each of them corresponds to one of the 
tie rotational modes in the E band. Although originat- 
ing from extended lattices, the concept of spin waves can 
be carried over in an obvious manner to any lattice which 
obeys the L&E band c oncept, not only those ch aracter- 
ized by a wave vector ( Cepas and ZimanL l2QQ5h . Quali- 
tatively, the E band reflects the possible higher-energetic 
internal spin structures which results as excitations from 
the classical antiferromagnetic ground-state spin config- 
uration. 

In the AFMHR the cyclic symmetry gives rise to a shift 
quantum number q defined via the shift operator T\q) = 
e iq \q) with q = 0, ±27r/7V, tt by which the excitations 
can be classified. A complete theory is not existing but 
phenomenologically the energies of the L and E band for 
S > 1 can then be approximated by 

E(S, q) = \e(ir)S(S + 1) + e(q) - e(0) (75) 

with q = 0, 7r for the L band. e(q) can be regarded as 
the f inite-size vers i on of t he spin-wave dispersion rela- 
tion (Dr eiser et al. L l2010bh . Indeed, in the infinite chain 
q would become the wave vector, and e(q) would agre e 
with the familiar spin- wave dispersion ( Anderson! 1 19521 ). 
However, it is emphasized again that the L&E'-band con- 
cept is not limited to clusters with cyclic or a similar high 
symmetry, it can also be observed in clusters with differ- 
ent point group symmetries or no symmetry at all (vide 
infra). 

We turn now to discussing the molecular wheel Crg. 
In this molecule, eight Cr 3+ ions (s = 3/2) form an al- 
most perfect planar octagon. The system crystallizes 
in space group P42i2, and the molecule nominally ex- 
hibits C4 symmetry. However, disorder in some pivalate 
ligands and tBu groups is present, suggesting that the 
individual molecules are slightly distorted. INS data 
yielded the exchange parameter 2 J = —1.46(4) meV 
and a nisotropy D = —0.038(5) meV ( Carretta et all 
l2003bl ), consistent with thermodynamic and EPR results 
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FIG. 27 (Color online) Q dependence of the integrated neu- 
tron scattering intensity of Crs for the peak (a) L° and (b) 
Ei . Both peaks arise from transitions from the S — ground 
state into the S = 1 sector, but the involved final spin multi- 
plets differ by their shift quantum number q, as indicated in 
the panels. T he solid lines rep r esent the theoretical curves. 
Adapted from Waldma nn et all 2003. 



(van Slageren et al 1 12QQ2D . A variation of the exchange 
and anisotropy parameters along the wheel as allowed by 
a C4 symmetry was not detected. Evidence for a weak 
rhombic term with E/D « 0.11 was found. 

A careful analysis of the INS dat a provided a de- 
tailed picture of the excitations in Crg ( Wa ldmann et all 
2003). The experimental INS spectra, corrected for non- 
magnetic scattering, are compiled in Fig. [26j Compari- 
son with the theoretical result shown in Fig. [25] will be 
made. The L° transition, or transition from the S = 
ground state to the lowest 5 = 1 multiplet, is split into 
two close peaks because of the zero-field splitting from 
the anisotropy, exactly as discussed before for NaFe6- 
The splitting is however small demonstrating the " weak" 
anisotropy case in Crg and justifying an analysis in terms 
of the AFMHR model. This is further corroborated by 
comparing the experimental INS spectra to the theo- 
retical correlation function, which agree in any detail, 
demonstrating the validity of the above point (1). 

The INS data allowed also a detailed compari- 
son of the excitation intensities or oscillator strengths 
KAIITW^^IIA')! 2 , and very good agreement between ex- 
periment and theory was found. Furthermore, magnetic 
scattering intensity at energies higher than that corre- 
sponding to transition E2 in Fig. [25] was not observed, 
demonstrating the selection rule in point (2). 

Finally, also the Q dependence of the peak intensi- 
ties were analyzed. The experimental results for the L° 
and Ei transitions and the theoretical expectations are 
presented in Fig. [27] The good agreement is obvious. 
However, most important, the Q dependence provides 
a fingerprint of the underlying many-body structures of 
the quantum states involved in a transition. Hence, the 
observed different Q dependencies of the L° and E\ tran- 
sitions directly demonstrate the different physical nature 
of the excitations in the L band (rotations of the Neel 
vector) and the E band (spin waves). 

Having demonstrated the L&^-band picture for Crg, 
the question arises of how general it is. A conclusive 
answer is not yet available, but the results so far sug- 
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FIG. 28 (Color online) Physical properties of the antiferro- 
magnetic Heisenberg ring as a function of its size N and the 
spin length s. In the blue shaded area quantum fluctuation are 
weak and the spin structure is essentially "classical" . Here the 
L&^-band concept becomes valid and Has a good effective 
spin Hamiltonian. This "classical" regime is reached when 
the size N is significantly smaller than a characteristic corre- 
lation length £ , which roughly inc reases exponentially with s. 
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Adapted from Waldmann, 2002b. 



gest Fig. [28j The LhE band c oncept works the better 
the larger s and the smaller N ([Engelhardt et all 120061 : 
IWaldmannL l2QQ2bh . The crossover is qualitatively deter- 
mined by comparing the size N of the cluster to a length £ 
characterizing the decay of the antiferromagnetic correla- 
tions with distance, as indi cated in Fig. [2 8 } A precise def- 
inition of £ can be subtle ([Affleck et al. □19871 iHalda^el . 
1983a), but in any (disordered) system the correlations 
decay on a characteristic length scale, which is what is 
meant by £. The spin length enters then through the de- 
pendence of £ on s. which, motivated by Haldane ' s and 
other results (lAffleck et all Il987l iHaldand . Il983al lbh. is 
roughly exponential. Since for s = 1/2, £ is about 6 
sites, it is concluded that large s = 1/2 clusters are dis- 
tinguished from those with s > 1/2 in that the LhE 
band picture is never obeyed in them. 

The L&£?-band concept has some useful consequences, 
valid not only for the AFMHR. First, the L band can well 
be reproduced by an effective Hamiltonian which may be 
constructed along the lines outlined in Sec. IIV.B.31 This 
is carried out now for the example of the even antiferro- 
magnetic wheels. Their lattice is bipartite and introduc- 
ing two sublattices A and B is natural, where the spins 
on each sublattice are ferromagnetically aligned. Then an 
effective Hamiltonian is obtained by combining the spins 
on each sublattice to = XlieA &i and S# = X^e£^> 
and inserting = (2/N)Sa or = (2/7V)Sb, depending 
on the spin's sublattice, into the microscopic Hamilto- 
nian, which results in 



Hab = — 2jS^ • Sb, 



(76) 



with j = a\ J. For a\ one finds a\ = a^ s , where af B = 
4/iV. The two-sublattice Hamiltonian Hab produces 
obviously an energy spectrum according to Eq. EH im- 



mediately explaining the step-like magnetization curves 
(Fig. [22]) and L band in the energy spectrum. The ap- 
proximation can be improved to yield nearly exact re- 
sults if a\ is slightly c orrected to accoun t for the weak 
quantum fluctuations ([Waldmannl l2002al ) . For Crg one 



obtains a\ 



with af 71 = 0.5586. 



Secondly, the nature of the E band suggests applying 
SWT to the calculation of its energies. This route was 
explored in the last years, but the status is not yet clear 
and more work is needed. However, it seems that despite 
the conceptual problems of antiferromagnetic SWT the 
standard interacting SWTs do, at least for bipartite lat- 
tices, produce reasonable results for the energies of the 
E'-band states in the S = 1 sector. For instance, for the 
AFMHR the "dispersion relation" e(q) in Eq. ((75]) are 
obtained (vide infra). 



2. Antiferromagnetic wheels with " intermediate" anisotropy 

In antiferromagnetic wheels with substantial magnetic 
anisotropy the ground state and lowest excitation may 
better be described by quantum tunneling of the Neel 
vector (Sec. IIV.C.3|h but for the next-higher lying lev- 
els the L&^-band picture is still appropriate, though S 
mixing occurs as a novel feature. In the following the 
effects of a significant anisotropy on the excitations will 
be discussed by the example of CsFeg, as it is one of 
the best characterized large magnetic clusters and the 
generic Hamiltonian Eq. ([73]) has been confirmed with 
high precision. 

The chemical formula of CsFeg is Cs [Fes teas] CI. Eight 
Fe 3+ ions (s = 5/2) form an almost perfect octagon 
with a Cs + ion at the center (Fig. [21]). Depending 
on the synthesis, the system co-crystallizes with differ- 
ent sol vent molecules in spac e groups P4/n, Pna21, and 
P21/n ([Saalfrank et fl/J . ll997l) . and the molecule exhibits 
ideal or approximate C4 symmetry, but a dependence of 
the magnetic excitations on the solvent was experimen- 
tally not observed. CsFeg is a member of a family of 
wheels which are distinguished by the templating cen- 
tral alkaline ion, and a magn e to-str u ctural correlation 
was established ([Pilawa et all 120031 : IWaldmann et all 
120011 ). The magnetic excitations in CsFeg were stud- 
ied by low-temperature high-field torque magnetome- 
try, single-crystal high-frequency EPR at Q-band and 
190 GHz , single-crysta l nucl ear magnetic resonance 
(NMR) ([Schnelzer et all l2QQ7h . and sever al INS exper- 
iments c overing energies up to 25 meV (|Dreiser et al , 
2010a .b: Sch nelzer etaA I2QQ7I : IWaldmann et al. , 2006b, 



2005. 



200 lh . 



The data analysis could be accomplished by solving the 
microscopic Hamiltonian Eq. ([73]) using the numerical 
sparse matrix techniques described in Sec. IIV.B.2I How- 
ever, the Hilbert space dimension is 1 679 616 and the 
task can become time consuming. Hence, approximate 
schemes are desired, with the advantage that their accu- 
racy can always be checked by comparing with the exact 
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FIG. 29 (Color online) Neutron-energy loss spectrum in 
CsFes, after background correction, at T = 17 K, and 
best-fit curve calculate d from Eq. ([77]) . Adapted from 
IWaldmann et all l2006bl . 



results from Eq. (f73|) for some few cases. The lowest-level 
description built on the strong-exchange limit works well 
for, e.g., EPR experiments, but involves many parame- 
ters and furthermore misses important effects (for a de- 
tailed description of the strong-exchange limit appro ach 
w e refer to the book by iBencini and Gatteschil Il990l or 
to lDreiser et a/ll2010a ). However, if only the energy lev- 
els in the L-band sector - to use the language of the pre- 
vious section - are desired, a high-accuracy higher-level 
description is provided by the sublattice Hamiltonian ap- 
proach. Applying the ideas used for the HDVV term 
(Sec. HV.C.ip to the anisotropy (and Zeeman) term in 
the microscopic Hamiltonian yields the sublattice Hamil- 
tonian 



H AB = -2jS A -S B +d (s% z + S% tZ ) + gn B S • B,(77) 

with j = a\J and d = b\ D, where a\ = af B and 
&! = bf B [bf B = (2s - l)/(Ns - 1)]. The approxi- 
mation can again be improved by adjusting a\ and b\ 
to account for the weak quantum fluctuations, which re- 
sults in essentially exact energies a nd yields tr a nsition 
intensitie s to within 1 0% accuracy (IWaldmannL [2002a; 
IWaldmann et all . l2006bl ). For CsFe 8 , af 1 = 0.5536 and 
bf™ = 0.1870. The energy spectrum of CsFeg is nearly 
identical to that of Crg shown in Fig.[25j and the labeling 
of states and transitions are carried over. 

INS experiments allowed the observation of all L-band 
states up to S = 5 at an energy of 14.4 meV [transitions 
L°-L 4 with respect to Fig. [25]. The data could very ac- 
curately be simulated, Fig. [29l Also, the EPR transi- 
tions observed in angle-resolved high-frequency EPR in 
the field range of - 12 T could be described with high 
accuracy, Fig. [30l In the EPR experiments the L-band 
states up to S — 4 were probed. In the fits to the data 
only two free parameters, 2 J and D, were involved. The 
determined values are given in Table. IVII (INS #2 and 
HFEPR). 

The excitation spectrum in CsFeg was also probed for 
energies up to 25 meV in a high-energy INS experiment 




(b) 



-180 -150-120 -90 -60 -30 30 60 90 120 150 180 
Angle ^(deg) 



.52 6 

u_ 

CO 

4 



}: T = 5K 


T= 15K 




^^^^ 



30 60 
Angle /?(deg) 



-30 30 60 90 120 150 18 
Angle /?(deg) 



FIG. 30 (Color online) Experimental (crosses) and simulated 
(doted lines) EPR transitions as observed in angle-resolved 
EPR experiments on single-crystals of CsFes at frequencies 
of (a) Q band (35 GHz) ( T — 25 K) and (b) 19 GHz (T = 5 
and 15 K). Adapted from lDreiser et aELl2010al . 
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FIG. 31 (Color online) Neutron-energy loss spectrum in 
CsFe 8 (T = 5 K). (a) Q-sliced INS data, with curves off- 
set for clarity. Gray curves represent the INS data recorded 
at T = 58 K after Bose-correction. The solid line represents 
the best-fit si mulation of the INS sp ectrum, (b) S(Q,lj) plot. 
Adapted from lDreiser et al. I, l2010bl . 



([Dreiser et all . l2010bl ). Results are shown in Fig. EU At 
lower energies the transitions within the L band are again 
observed, however, at energies of ca. 7.5 and 10 meV two 
further cold transitions are detected, which can unam- 
biguously be related to the discrete spin- wave excitations 
Ei and expected in an octa-nuclear antiferromagnetic 
wheel. Above transition E2 no magnetic scattering inten- 
sity is observed, confirming the selection rule associated 
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TABLE VI Comparison of the magnetic parameters for CsFes 
obtained by different experimental techniques. 



Technique 2 J [meV] D [meV] 



Reference 



torque -1.90(10) -0.045(3) 

INS #1 -1.78(4) -0.048(1) 

INS #2 -1.80(2) -0.050(1) 

INS #3 -1.79(5) -0.050(7) 



Waldmann et a/., 2001 
Waldmann et aL 2005 
Waldmann et aL 2006b 
Dreiser et aL 2010b 



HFEPR -1.87(25) -0.0493(1) Dre iser et aL . 



to the L&^-band picture in Sec. IIV.CTTl The determined 

2 J and D values are listed in Table. ED (INS #3). 

The INS data were recorded on a non-deuterated poly- 
crystalline sample. This may explain the huge non- 
magnetic scattering in Fig. [3TJ which is typically ob- 
served in molecular nanomagnets above energies of ca. 2- 

3 meV. For CsFes it could however very well be accounted 
for using a Bose-correction analysis (Sec. IILB.l|h which 
allowed the unambiguous identification of the magnetic 
peaks, see Fig. l3Tla). 

Table IVII compiles the determined 2 J and D values, 
including those o btained by torqu e and h igh-re solution 
INS (INS #1) (IWaldmann et all I2QQ5L l2QQll) . The 
consistency is excellent, in particular considering the 
large range of energy scales probed in the experiments 
(-0.01 meV in 35 GHz EPR, -10 meV in high-energy 
INS). Efforts were made to infer the significance of fur- 
ther terms in the microscopic spin Hamiltonian not in- 
cluded in Eq. ([73]) . A J\ — J2 modulation of the ex- 
change constants along the wheel was found to be smaller 
than 20%, and the r hombic anisotropy to be negligible, 
E = 0.0000(3) meV ([Preiser et all l2010allbh . 

The determined 2 J value in the EPR experiment de- 
serves a comment. The EPR selection rules do not allow 
a direct observation of exchange splittings (Sec. IILB.3|) . 
Exchange constants may though be determined indi- 
rectly, through the temperature dependence of the EPR 
resonance intensities, which however in large clusters 
i s challenging, or through the ff-mix i ng mechanism 
dBarra et a/1 (20071: iLiviotti et a/1, 120021: iPiligkos et all . 
l2009t IWaldmann and Giidell . 120051 : IWilson et all 120061) . 
which was the case in CsFes. In the strong-exchange 
limit, the anisotropy splitting produces the "normal" 
zero-field splitting pattern, e.g., a D term in the micro- 
scopic spin Hamiltonian produces a zero-field splitting of 
the spin multiplet, which follows the M 2 behavior cor- 
responding to S 2 . However, if anisotropy is stronger, as 
compared to the exchange J, then the pattern is modified 
and deviations from M 2 occur. In perturbation theory 
this corresponds to higher-order terms (S 2 ) n with n > 1 
coming in, with weights oc (D/ J) n . Detecting these shifts 
in the zero-field splitting pattern allows the indirect de- 
termination of the strength of the exchange. The excel- 
lent agreement of the EPR 2 J value demonstrates hence 
that the generic Hamiltonian Eq. ([73]) does also describe 
the subtle ^-mixing effects very well in CsFes. 




FIG. 32 (Color online) 190 GHz EPR spectrum on a single 
crystal of CsFes. The intra- multiplet transition a is clearly 
observed. The inset to the left depicts the three lowest en- 
ergy levels and the allowed INS transitions I, II, and a. The 
inset to the right shows INS data recorded at T — 9.7 K in 
which the transition a should have be en observed if it wer e 
of appreciable intensity. Adapted from Dreise r et al , 2010a. 



The comparison of the INS and EPR experiments re- 
veal another interesting aspect of the excitations in an- 
tiferromagnetic wheels. According to the INS selection 
rules the transition \S = 1, M = 0) <R> \S = 1, M = ±1), 
or a henceforth, is allowed and should be detected at 
appropriate temperatures, yet it is not observed in INS 
experiments, albeit in EPR experiments (Fig. [32]) . It 
turns out that m£ra-multiplet transitions, such as a, are 
orders of magnitude weaker than m£er-multiplet tran- 
sitions because of the particular many-body structure 
of the wave functions, which is that of a bipartite lat- 
tice of two mesoscopically s ized spins on each subla ttice 
(S A = S B = 10 in CsFe 8 ) ([Waldmann et all l2QQ5h . By 
that reason the m£ra-multiplet transitions become too 
weak to be observed by INS. CsFes provides a convincing 
example, since here the transition a should have easily 
been detected by INS, see Fig. [321 Hence, the combined 
INS and EPR data demonstrate directly a hallmark fea- 
ture of antiferromagnetic wheels, namely the mesoscopic 
antiferromagnetic sublattice structure. 

Although anisotropy is appreciable in CsFes, the en- 
ergies of the spin-wave excitations E\ and E2 are actu- 
ally little affected. This, and fundamental interest, mo- 
tivated an analysis of the elementary excitations of the 
N = 8, s = 5/2 AFMHR model using different variants 
of SWT. All models provide predictions for the £?-band 
states or the "dispersion relation" e(q) in Eq. ([75]) , but 
only the last three yield also estimates for the singlet- 
triplet gap A (= ai|2J|) or L-band indeed. The findings 
are compared in Fig. [33] to the exact energies. Interest- 
ingly, all SWT models reproduce roughly the spin-wave 
excitation spectrum, which supports the notion that the 
L&£?-band concept is essentially classical. However, sig- 
nificant differences exist, and the ISWT and LSWT+A C 
models do best for the £?-band excitations. As regards 
the singlet-triplet gap, the LMSWT and IMSWT esti- 
mate it almost a factor of two too small. LSWT+A C 
comes closest, to within 7%. A similar analysis for the 
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FIG. 33 (Color online) Zero-temperature excitation spectrum 
of the N = 8, s = 5/2 AFMHR, as function of the shift quan- 
tum number q as calculated using exact numerical diagonal- 
ization (stars) and the indicated SWTs (lines) (in the latter q 
was assumed as continuous for clarity). LSWT = linear SWT, 
ISWT = interacting SWT, IMSWT = full-diagonalization in- 
teracting modified S WT, LSWT+A C = li near SWT with a 
shift. Adapted from Dreis er et al 1 l2010bl to which we also 
refer for details. 



larger Fe is wheel (N = 18, s = 5/ 2) confirmed the ob- 
servations ( Ummethum et a/ll2012l ). 



3. Antiferromagnetic wheels with "strong" anisotropy and 
quantum tunneling of the Neel vector 

The possibility of quantum tunneling of the Neel 
vector (QTNV) in antiferromagnetic materials at- 



tracted huge interest (Barbara and Chudnovskvl 1990; 



Gider et all 



1995 



Chudnovskv and Teiadal 11998 

Gunther and Barbara! . [l995fc E rive and Zaslavskiii 

1990; Shuvrko et aR I2007D . and initial attempts to 
establish QTNV concentrated on ferritin proteins 
(lAwschalom et all Il992t iGard . Il996b iGider et all 1996; 
iTeiadaL 11996^ The prediction that coherent QTNV 
might also be realized in antiferrom agnetic molecular 
whee ls with strong anisotropy ( Chiolero and Loss! 
1998) stimulated intense research. In this context the 
molecules Fei fSec. HVX]) . CsFe 8 fSec. HV.C.2|) . and Fei 8 
appeared as promising candidates (Fig. [21]). Reviews 
and follow-up articles on the QTN V scenario in antifer- 
romagnetic wheels are availa ble (Konstantin idis et all 
l201lULeuenberger et all 120031 : [Meier and Lossi l2001bh . 

The discussion in the following refers to Eq. ([73]) . Due 
to the anisotropy the spin functions \SM) are mixed 
strongly and S loses its significance as good quantum 
number. This changes also the dynamics of the Neel 
vector, defined as n = (S^ — Sb)/(Sa + Sb), which is 
then not anymore that of a rotation but tunneling. In 
the quantum spectrum this may be seen by the fact that 
the two lowest levels (ground state & first excitation) ap- 
proach each other but separate from the next higher- lying 
levels, as shown exemplarily in Fig. IMT d): A01 becomes 
smaller than e.g. A12. 

A semi-classical theory provided a clear description 
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FIG. 34 (Color online) (a,b) Classical ground-state spin con- 
figuration for a (a) N = 8, (b) N = 18 antiferromagnetic 
wheel with the Neel vector (long green arrow) pointing along 
(a) z, (b) y. (c) Shape of the potential y(n) in the case 
(B). The two tunneling paths from n = -\-y to — y are in- 
dicated. Simulated low-lying energies of a N = 8, s = 5/2 
antiferromagnetic wheel vs D/(2J) in zero field (M is then 
an exact quantum number) . Transitions I and II are observed 
in INS experiments (Fig. [35]). (e-g) Energy scheme as given 
by semiclassic al theory for (e) Feis, (f) C s Fes, and (g) Feip . 
Adapted from Konstantinidis et al, 2011; Waldmann et a/., 
20051. 120091 . 



([Chiolero and Loss] . Il998fc iMeier and Loss! l2Q01bh - De- 
pending on the magnetic parameters and strength and 
orientation of the magnetic field, several scenarios occur; 
we here focus on two: (A) D < and zero magnetic field 
and (B) D > and large magnetic fields with B along the 
x axis. In both cases, the Neel vector is strongly local- 
ized along two directions, namely n = d=z (A) or d=y (B). 
Classically, the ground state is then characterized by the 
spin configurations with Neel order sketched in Fig.lMTa) 
and (b), respectively, and rotation is hampered by an en- 
ergy barrier, corresponding to a potential V(n) with min- 
ima at the respective orientations [Fig. l34t c)]. However, 
quantum fluctuations allow for tunneling of the Neel vec- 
tor, which lifts the classical degeneracy and opens a tun- 
neling gap Aqtnv in the energy spectrum, which then 
corresponds to A01. 

The Neel- vector dynamics is characterized by the tun- 
neling action So/h, attempt frequency huo, barrier height 
AU, and tunneling amplitude Aq, 



So/h = Nsy/2\D/(2J)l 



(78) 
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AU 



sy/8\D(2J)\, 
Ns 2 \D\, 



s /h 

2tt 



exp(-So/h). 



(79) 
(80) 

(81) 



(a) 



The two cases (A) and (B) need to be distinguished 
now, which we do through a parameter c = 2 (A) and 
c — 1 (B), respectively. The ground state energy is then 
given as ffk^o, and the tunneling splitting as cAq [in 
the case (A) the tunneling gap has to our knowledge 
not yet been calculated, but numerical results suggest 
Aqtnv ~ 2Ao]. The semi-classical theory for QTNV 
becomes valid for large tunneling actions j^So/h ^> 1, 
which is equivalent to stating that the ground-state en- 
ergy is smaller than the barrier height, AU ^> fftc^o, or 
the Neel vector strongly localized, (0\h z / y \l) 2 — » 1, where 
the matrix element is given as (0\h z / y \l) 2 « l — c/(So/h), 
where z (A) or y (B) is the respective Neel vector com- 
ponent. 

As regards applying the semi-classical theory some 
points are worth noticing. First, it turned out that 
the J and D values as they appear in Eq. ([73]) or 



through 



and 



b qm 



in Eq. (|77|) should not be in- 
serted in the semi-classical formulae, but "corrected" J, 
D values as they are obtained by using a^ c = 4/N and 



6f = 2Ns 2 /[Ns(Ns + 2) - 3] (Waldm ann et all 1200 



(for a detailed discussion see lKonstantinidis et aZ.L 1201 lh 
Not doing so leads to e.g. significantly overestimated 
tunneling actions, which went unfortunately unnoticed 
in early works. Furthermore, the crossover from weak to 
strong anisotropy, or from rotation to tunneling of the 
Neel vector, is continuous and not abrupt [Fig. I34H)], 
and the QTNV scena rio is hence necessarily approximate 
([Santini et all l2QQ5h . This introduces some ambiguity, 
and different criteria for when QTNV is realized can be 
given, e.g., that the tunneling levels should fall below the 
top of the barrier or that the tunneling splitting should be 
exponentially small. It is noted in passing that since the 
semi-classical theory as it stands is an approximate the- 
ory, agreement with semi-classical theory is a sufficient 
but not necessary criterium for QTNV. 

The low-lying excitations in zero fie ld were mea- 
sured by INS for Fein, CsFeg, and Fei8 ( Santini et all 
120051 : IWaldmann et all 120051 l2QQ9h . Fei and CsFe 8 
were described before. The chemical formula of Fei8 
is [Fe 18 (pd) 12 (pdH) 12 (0 2 CEt) 6 (N0 3 )6](N03)6^MeCN 
« 48); the eightee n Fe 3+ (s = 5/2 ) ions are arranged 
in a cycle (Fig. [21]) ([King et a/ll2006[ ). The system crys- 
tallizes in space group R3- and the molecule exhibits crys- 
tallographic S6 symmetry. One nitrate and ca. 8MeCN 
solvent molecules are disordered. INS data for the three 
wheels are shown in Fig. [35j The lowest excitation I 
could clearly be detected as well as several higher excita- 
tions (for CsFes see Sec. IIV.C~2]h which facilitated a pre- 
cise determination of the magnetic parameters. In Feio 
a substantial rhombic term \E/D\ = 0.21 was found, and 
structural disorder had to be included in the analysis. In 
Feig a high-energy INS experiment evidenced a modula- 
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FIG. 35 (Color online) Neutron-energy loss spectra for (a) 
Feio [T = 2 K (blue), 5 K (red), 10 K (green)], (b) CsFe 8 
(T = 2.4 K), and (c) Fei 8 [T = 1.9 K (blue), 4.2 K (red)] 
showing transition I, which in the QTNV regime correspond 
to the Neel vector tunneling splitting Aq i . Ad apted from 
ISantini et all I2QQ5I : IWaldmann et all I2005L 120091 . 



tion of the exchange constan ts along the ring con sistent 
with the C3 symmetry (jUmmethum et all , l2Q12h . No- 
tably, the ratio of the excitation energies of transitions 
I and II decreases in the sequence of Feio, CsFeg, and 
Fei8- For better comparison of the wheels, the magnetic 
parameters reported in the original works were converted 
to 2j and d of Hab, Eq. ([77]) (in case of Feio the rhombic 
contribution was neglected and for Feig the appropriate 
averaged J was used) . The results are listed in Table EID 
For both CsFeg and Feio one finds d < or D < 0, 
while for Feig D > 0. Hence, in CsFeg and Feio [case 
(A)] the transition I observed by INS would directly 
correspond to the Neel vector tunneling gap if QTNV 
were realized in them. The tunneling gaps estimated by 
the semi-classical theory roughly agree with the observed 
gaps, but the other parameters compiled in Table [VII] in- 
dicate that QTNV is not well realized in Feio, and that 
CsFe 8 is borderline. For D > as in Fei 8 QTNV does 
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TABLE VII Experimental values and characteristic param- 
eters of QTNV for Fei s (c 1), CsFe 8 (c = 2) and Feio 
(c = 2). Adapted from ([Konstantinidis e£ aZ.L l201lh . 





Fei8 


CsFe 8 


Feio 


2 j [K] 


-5.1 


-11.1 


-6.31 


d [K] 


0.021 


-0.104 


-0.0276 


Aoi [K] 


not measured 


5.92 


3.83 


So/h 


5.90 


4.03 


3.42 


AU [K] 


22.2 


22.5 


9.25 


§hu [K] 


7.52 


22.3 


10.8 


cA [K] 


0.320 


5.08 


4.18 


(S /h)/c 


2.95 


1.01 


0.96 


(0\n z /y\l) 2 


0.83 


0.50 


0.42 


Aqtnv/AU 


0.014 


« 0.23 


« 0.45 
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FIG. 36 (Color online) Magnetic torque vs field for B nearly 
perpendicular to z (angle = 80.5°, T = 0.1 K), showing "wig- 
gles" due to the oscillations of the Neel- vector tunneling gap. 
The inset shows the simulated field-dependent torque and the 
tunneling gap as calculat ed quantum mechanicall y and semi- 
classically. Adapted from IWaldmann et all 120091 . 



not occur at low fields (|Chiolero and LossL Il998f h but 
in the high-field regime, which in Fei8 is reached above 
10.6 T [case (B)]. Here, QTNV is well realized in Fe i8 
as e.g. indicated by the exponentially small tunneling 
gap. The situation in the three wheels is probably most 
clearly demonstrated by the energy diagrams shown in 
Fig. IHe)- (g). 



The high-field Neel vector tunneling gap in Feig was 
not directly observed via spectroscopic techniques, but 
the field-dependent oscillations in the Neel vector tun- 
neling gap due to quantum interference were detected 
in low- temperature high-field magnetic torque measure- 
ments (jWaldmann et all 12009). In the semi-classical the- 
ory for case (B) the phase of the ground-state wave func- 
tion contains a topological term irNgfiBB /\4(2 J)|, which 
is proportional to the field. Hence, in increasing fields the 
phase is repeatedly tuned through destructive and con- 
structive interference, which gives rise to an oscillation 



of the tunneling splitting according to 

_Ngn B B\ 



A(B) = A QTNV 



sin 7T 



4|2J| J 



(82) 



Since the tunneling splitting also affects the ground-state 
energy, Eo(B) = e(B) + A(£?)/2, where e(B) is a smooth 
function, the oscillations can be detected by magneti- 
zation or torque measurements at low temperatures. In- 
deed, the observed "wiggles" in the torque curve (Fig.l36|) 
do directly correspond to the oscillations in the tunneling 
gap, which is demonstrated e.g. by comparing the numer- 
ically calculated curves for the torque and the tunneling 
gap (inset to Fig. [36]). The analysis also showed that the 
semi-classical theory yields highly accurate results, which 
underpins the notion of QTNV in Feis- 

It is added that wiggles in the torque as function of 
field can also occur due to a first-order mixing of the 
\S,M = -S) and \S + 1, M = -S - 1) spin levels at 
the field-induced level-crossing points by the magnetic 
anisotropy; an example for this ^-mixing mechanism will 
be presented below in Sec. HV.D.ll However, theoret- 
ically QTNV cannot be described by ^-mixing of two 
spin levels reflecting its different underlying physics, and 
also experimentally, the two mechanisms can unambigu- 
ously be distinguished from each other. For instance, 
the dependence on the angle between magnetic field and 
anisotropy axis z allows a clear-cut decision: In the S- 
mixing scenario the wiggles occur for both nearly par- 
allel and perpendicular fields, in contrast to the obser- 
vations in Feig, where ordinary staircase-like profiles are 
observed for parallel fields, as predicted by the QNVT 
scenario. 



D. " Modified" antiferromagnetic molecular wheels 



The ring topology considered in Sec. HV.CI can be var- 
ied in a number of ways, by "slight" modifications. The 
topologies which will be addressed here can be put into 
four categories. First, one of the magnetic metal ions in 
the ring, which carry spin s, is replaced by another mag- 
netic metal ion with different spin so ^ s. These clusters 
shall be denoted as "doped wheels", and the foreign so 
ion as impurity. Second, the cyclic boundary conditions 
are changed to open boundaries, e.g., by replacing one of 
the magnetic ions in the ring with a diamagnetic ion or 
by removing one metal center. These clusters shall be de- 
noted as "short chains". Third, a magnetic ion is added 
at the center of the ring with coupling paths such that 
the lattice remains bipartite. And fourth, a magnetic ion 
is added at the center with the coupling paths introduc- 
ing spin frustration, then called "discs". Representative 
examples are Cr^Ni, Cr6, Mn-[3x3], and a Fe7 molecule, 
shown in Fig. [371 

Much interest in such systems came from the theo- 
retical suggestions that an uncompensated spin intro- 
duced into an antiferromagnetic wheel e.g. by dop- 
ing may act as a tracer spin for the quantum tunnel- 
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(a) 



Mn-[3x3] 



FIG. 37 (Color online) Molecular structures of the doped 
wheel Cr 7 Ni: [ n Pr 2 NH 2 ] [Cr 7 NiF 8 Pivi 6 ], short chain Cr 6 : 
[R2NH 2 ]3[Cr 6 FiiPivio(H20)], Mn-[3x3] grid representing a 
wheel with magnetic center: [Mng(2POAP-2H)6], and disc 
Fe 7 : [Fe 7 03tea3Piv9(H 2 0)3]. 



ing of the Neel vector, which in this way could ex- 
perimentally be observed and possibly manipulated by 
EPR methods, which othe rwise would not be possi- 
ble (|Meier and Lossl bOQlal lbl). Furthermore, the excess 
spin may have S f = 1/2, which then sugg ests applica- 
tion o f t he cluster as quantu m bit (qubit) ([Meier et all 
l2QQ3al lbl: iTroiani et all . l20Q5h . Using such "mesoscopic 
spin- 1/2" clusters as "spin cluster qubits" could pro- 
vide advantages such as easier addressing. Significant 
progress in this direction have been made on the Cr 7 Ni 
wheel, e.g., long coherence times were observed, mag- 
netic coupling between two Cr 7 N i clusters introduced , 
and entanglement demonstrated ( Ardayan et all 120071 : 
ICandini et all l2010fc Pfimco et aZ.L l2QQ9h . Reviews are 
available ([Afironte et a/1 120Q5L l2006h - Since the physics 
is related mainly to the ground state, these exciting de- 
velopments will however not be discussed. We will fo- 
cus here on the basic question of the impact of topol- 
ogy on magnetic cluster excitations. The Mn-[3x3] grid 
molecule will be considered first, followed by the doped 
wheels and short antiferromagnetic chai ns. An abun- 
dance of discs have been synthesized ([Hoshino et all 
2009a), but detailed studies of the cluster excitations 
were to our knowledge not reported. 
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FIG. 38 (Color online) (a) Neutron-energy loss spectra of 
Mn[3x3] (2) at 2.5 K in two different energy ranges, (b) 
Simulated energy spectrum with anisotropy neglected. The 
observed INS transitions are indicated by arrows. The in- 
set sketches the observed transitions within the S = 5/2 
ground and first excited S = 7/2 multiplet. Adapted from 
iGuidi et all I2004 



most interesting from the perspective of cluster excita- 
tions. The molecule may be crystallized with different 
counter ions and solvents, yielding e.g. [Mng(2POAP- 
2H) 6 ](C104)6-3.75CH 3 CN11H 2 (1) or [Mn 9 (2POAP- 
2H) 6 ](N0 3 )6-H 2 (2). Nine Mn 2 + (s = 5/2) metal ions 
are arranged on the vertices of a 3x3 grid (Fig. [37]). They 
crystalize in space group C2/c, and the molecules exhibit 
a slightly distorted symmetry with the S4 symmetry 
axis (= z axis) perpendicular to the grid plane. Consid- 
ering the symmetry, the appropriate spin Hamiltonian 
for describing the magnetism reads 

H = -2J R ^Y] 8j • Sj + i + s 8 • si^ 
-2 J c (s 2 + s 4 + s 6 + s 8 ) • s 

8 

+D R *h + Dcslz + MbS ■ B, (83) 



1. The Mn-[3x3] grid molecule 



Molecular [nxm] grids attracted considerable inter- 
est in chemistry, for the employed preprogrammed self- 

and their physical prop ertie 
3, l2009t iRuben et all [200 
the m agnetic properties were re viewed in | Waldman: 
l2005ai The Mn-[3x3] grid (|Zhao et all l2000l ) is 



assembly synthes is strategy a n d thei r physical properties . 
For reviews see iDawe et all 120091 : IRuben et a/1 |2004 



where the spins are numbered as given in Fig. l37f c). In 
principle, the D values for the corner and edge Mn ions 
could be different, but this was found to not affect the 
magnetism significantly. 

From magnetization measurements, antiferromagnetic 
exchange interactions and a S = 5/2 ground state were 
inferred, which can be rationalized by the classical spin 
configuration, where corner and central spins point up 
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12 16 20 24 28 

Magnetic field (T) 

FIG. 39 (Color online) Magnetic torque vs field (upper curve) 
and the field derivative (lower curve) for Mn-[3x3] (1) with B 
nearly perpendicular to the magnetic z axis (angle = 80.5°, 
T — 0.1 K), showing the oscillatory torque behavior. The 
inset shows the energies of t he spin multiplets as ex tracted 
from the level- crossing fields (Wa ldmann et aEll201l[ ). 



and edge spins point down. The excitation spectrum up 
to en ergies of 4 meV was studied by INS ( Guidi et all 
2004). Some results are presented in Fig. l38Ta). The 



dimension of the Hilbert space is 10 077 696, and so- 
phisticated numerical approaches had to be developed 
for analyzing the data. Good agrement was obtained for 
2 J R = 2J C = -0.43 meV and D R = D c = -0.012 meV. 
The higher-lying excitations revealed a small deviation of 
the exchange constants from the S4 symmetry assumed in 
Eq. ([83]). The calculated energy spectrum with D R and 
Dq set to zero is presented in Fig. 155Tb). The two low- 
est transitions a and j3 stem from the zero-field splitting 
of the S = 5/2 ground state, as sketched in Fig. l38Tb). 
Peaks la, lb, and Ic correspond to transitions from the 
zero-field splitting levels of the ground multiplet to the 
next higher- lying S = 7/2 multiplet, and peaks E\ a and 
Eib go from the ground multiplet to the lowest S = 3/2 
multiplets. 

Besides that the INS data of such a large cluster as 
Mn-[3x3] were successfully interpreted, the inspection of 
the determined energy spectrum is interesting: As in an- 
tiferromagnetic wheels, the excitations may be classified 
as L and E band, or Neel- vector rotation and spin waves. 
This observation can be linked to the bipartite topology 
of the grid lattice, suggesting a "classical" spin structure. 

The L band h as been demonstr a ted in magnetic torque 
measurements ( Waldmann et all l2004f ). At very low 
temperatures the torque exhibits an oscillatory field de- 
pendence, Fig. [39j which can be related to a sequence 
of level crossings, where the ground state changes from 
S = 5/2 to S = 7/2, 7/2 — >• 9/2, and so one, similar to the 
situation in the wheels sketched in Fig. [23j In the exper- 
iments all states of the L band up to S = 23/2 were ob- 
served, and the energies of the multiplets as determined 
from the level-crossing fields do indeed follow the Lande 
pattern S(S + 1), as displayed in the inset to Fig. [39j 
Hence, the L band (or tower of states or quasi-degenerate 



joints states), which is the prec ursor to long - range Neel 
ordering in the infinite lattice ( Bernu et al. < Il992h has 
been experimentally demonstrated for the square-lattice 
topology. 

These findings suggest the construction of an effec- 
tive spin Hamiltonian for antiferromagnetic 3x3 grids. 
Applying the procedure used for the antiferromagnetic 
wheels to the eight metal ions on the periphery yields 

H abc = —^3r^a • &b — 2 JcSo • S# 



+d R (s% z + + D c s\ z 

+#//£?S • B, 



(84) 

with the sublattices A = {1, 3, 5, 7} and B = {2, 4, 6, 8}, 
and j R = 0.526J R and d R = 0A97D R for the Mn-[3x3] 
grid. Habc was indeed demonstrated to produce highly 
accurate results, a nd was found crucial in the analysis of 
experimental data dDatta et al. U2007l : IWaldn^m n. 2005b; 
IWaldmann et all l2004h ~ 

The oscillations in the torque signal or iginate from an 
interesting quantum-mixing m echanism ([Carretta et all 
l2003at IWaldmann et all 120041 ) . A magnetic anisotropy 
is generally expected to induce mixing of spin multiplets 
(S mixing), which often may be treated perturbatively, 
implying a "small" effect (but see also Sec. IIV.CT3l) . How- 
ever, if two states are close in energy, i.e., essentially de- 
generate, then standard (non-degenerate) perturbation 
theory will obviously break down, and the effect of the 
perturbation not be small or mixing of the states large. 
The anisotropy produces thus a strong mixing of the 
\S,M = -S) and \S + 1, M = -S - 1) states involved 
at a level crossing, and the ground state is described as 
a superposition 



\g)cxa(B)\S) + b(B)\S + l) 



(85) 



where the field-dependent a and b become equal at the 
level crossing (an obvious shorthand notation for the 
states was used; a 2 + b 2 = 1). Since states with dif- 
ferent total spin are mixed, the total spin will fluctuate 
strongly, also ca lled quantum oscillations of the total spin 
( Carretta et all l2003af ). This mixing is directly related 
to the oscillatory torque curve observed in experiment, 
which hence provides evidence for this phenomenon. 

S mixing should in principle also enable a direct de- 
tection of exchange splittings through EPR experiments, 
since the EPR selection rule AS = would not hold ex- 
actly. However, the mixing is usually not strong enough 
for such EPR transitions to gain sufficient intensity, but 
through this mechanism the S = 5/2 ^ S = 7/2 transi- 
tion could directly be observed in Mn-[3x3] as function of 
field in a mu lti- frequency single-crystal EPR experiment 
(jDatta et aZ.L l2QQ7h . 



2. "Doped" even- mem bered antiferromagnetic wheels 

A series of octa-nuclear hetero-nuclear wheels of gen- 
eral chemical formula [H 2 NR 2 ] [Cr 7 M'F 8 Pivi 6 ], or Cr 7 M' 
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in short, wit h e.g. M 



= C u 2+ Ni 2 +, or Mn 2 + were 
2QQ7MBaker et all l2011bl : 



synthesized (lAffronte et .. . , , _ _ 

lLarsen et all 2003a; Lave et all . 20051 ). and their excita- 



tions studied by different techniques, am ong which were 
low-t emperature torque mag netometrv (|Carretta et all 
l2005|) high-frequency EPR hPilizkoseUdl l2009h . and 
INS ([Baker et all l2011bl : ICaciuffo et all l2005f ). Also 



an analogous Fe 7 Mn cluster was investigated using INS 
( Guidi et all l2QQ7h . As mentioned before, the molecule 
Cr 7 Ni has attracted most interest, because of its poten- 
tial use in quantum information, and is focused on here. 

The synthesis strategy resulting in Cr 7 Ni is ex- 
tremely flexib le, and potent i ally m any derivatives of 
Cr 7 Ni exist ([Affronte et all l2007h . The material 
used in INS experiments was of chemical formula 
[H 2 N(C 2 D5)2][Cr 7 NiF 8 (0 2 (C 5 D 9 )i 6 ] and crystallizes in 
space group P4, without solvent molecules in the crystal 
lattice. Seven Cr 3+ (s = 3/2) ions and one Ni 2+ (s = 1) 
ion form a ring as shown in Fig. [371 The material has 
two advantageous properties for INS; it can be deuter- 
ated to a large extent and large single crystals can be 
grown. The appropriate generic spin Hamiltonian reads 



H = -2 J ^] Si • s i+1 - 2J r (si + s 7 ) • s 

i=l 
7 



where the spins Si to s 7 refer to the Cr 3+ ions, and so is 
the spin of the Ni 2+ ion. In view of the molecular sym- 
metry, also a rhombic term Ei(s^ x — sf y ) is expected 
and w as evidenced in EPR experiments ([Piligkos et all 
I2009D . but w as not resolved in magnetic torque and IN S 
experiments ( Caciuffo et a/.L [2005; Ca rretta et aD,|2005). 



The exchange interactions in Cr 7 Ni are antiferromag- 
netic and the ground state is S — 1/2, consistent with 
the expectation from the classical spin configuration with 
Neel order. The magnetic excitations were studied by 
two INS experiments, which shall be described in the 
following. In a first experiment th e excitations up to en - 
ergies of ~4 meV were detected ([Caciuffo et all [2005). 
A spectrum is presented in Fig. l40la). One large cold 
feature, with a double-peak structure, is observed at 
1.27 meV, and a hot peak at 2.08 meV. At higher en- 
ergies evidence for further magnetic scattering intensity 
is observed. The lowest feature can be associated to 
the transition from the S = 1/2 ground state to the 
lowest-lying S = 3/2 multiplet, and the double-peak 
structure of this feature to a zero-field splitting in the 
S = 3/2 level. The zero-field splitting (0.15 meV) is 
much smaller than the center of gravity of the S = 3/2 
multiplet (1.27 meV); anisotropy is hence weak in this 
molecule, which is of relevanc e for i ts potential appli- 
cation as qubit (iTroiani et all 120051 ). The 2.08 meV 
feature corresponds to an S = 3/2 — >• 5/2 transition. 
The data could well be simulated using Eq. ([86]h yielding 
2J = -1.46 meV, 2J' = -1.69 meV, D lig = -0.03 meV, 
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FIG. 40 (Color online) (a) Neutron-energy loss spectra on 
Cr 7 Ni at T = 2 K and 12 K. Circles represent the ex- 
perimental data; line s the simulation result. Adapted from 
Caciuffo et at. |, 120051 . (b) Simulated energy spectra with 
anisotropy neglected. Observed INS transitions are indicated. 



and D hg ' = —0.6 meV. The simulated energy spectrum 
with anisotropy neglected is shown in Fig. BUTb). The 
spectrum reveals again a L-band structure, which is con- 
firmed by t he INS experiment, and high-field torque mag- 
netometry ( Carretta et all |2005[ ). Also £?-band states 
are detected in the INS data. 



The unique advantage of the Cr 7 Ni system, that large 
deuterated single-crystals can be grown, allowed the di- 
rect experimental observation by INS of t he level-crossing 
beha vior as function of a magnetic field ( Carretta et all 
|2007[ ). In this experiment, a crystal of 0.4 g weight was 
measured at T = 66 mK with magnetic field applied in 
a range of - 11.5 T and an angle of = 50° with re- 
spect to the z axis. Experimental results are shown in 
Fig. [HI At the level crossing at B c = 10.5 T a small 
gap of 0.12 meV is observed, i.e., an avoided level cross- 
ing, which originates from the S mixing induced by the 
weak anisotropy effects. As explained in Sec. IIV.D.11 at 
the level crossing where states become almost degener- 
ate, the mixing effect will be strong, the wave function 
described by Eq. ([85]) with \a\ = |6|, and the total spin 
oscillates. This INS experiment hence provides direct ev- 
idence for the quantum oscillations of the total spin in 
Cr 7 Ni at the level crossing. 
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FIG. 41 (Color online) (a) High-resolution neutron-energy 
loss spectra on a crystal sample of CrrNi as function of mag- 
netic field (T = 66 mK, = 50°). (b) Simulated low-lying 
energy levels as function of magnetic field (lines) with all ob- 
serve d INS transition ene rgies included (squares). Adapted 
from ICarretta eioH 120071 . 



3. Short antiferromagnetic chains 

The synthesis strategy which allowed us to gener- 
ate the doped antiferromagnetic wheels in Sec. IIV.D.2I 
also afforded the generation of short antiferromagnetic 
chains, either by replacing one magnetic metal ion 
in the ring-like structure by a diamagnetic ion such 
as Zn 2+ or Cd 2+ , or modifying the syn thesis method 
to yi eld structures called "horseshoes" ([Affronte et all 
120071 ). The members of this family of clusters were 
studie d, besides the the rmodyna mic techniques, b y 
EPR flPiligkos et a!l.l2QQ7h. XMCD dGhirri aZl l2QQ9h. 
NMR (lAmjri et all boidt Micotti et aUl2006h. and INS 
dBaker et all 12011 at iBianchi et all 120091 : ICaciuffo et all 
I2QQ5I : lOchsenbein et all I2QQ81 l2QQ7h . Short antifer- 
romagnetic chains were also obtained through doping 
the Fei8 wheel (Sec. | IV.C.3|) with diamagnetic Ga 3+ 
([Henderson et all l2QQ8h . 

From the view point of the physics of the mag- 
netic excitations, comparison of short antiferromagnetic 
chains with even-membered antiferromagnetic wheels, 



or chains with open and periodic boundary condit ions, 
should prove interesting ( Ochsenbein et all , 120071 ), as 
well as comparing antiferromagnetic chains with even 
and odd number of metal centers ([Ochsenbein et all 
l2008f ). A large body of literature exists on one- 
dimensional q uantum spin chains, and also finite chains 
were studied ([Bogani eTo/1 l2QQ4t iDi Tusa et all . Il994t 

ap- 



1990). It 



Fujiwara et all 11998c lHagiwara et all 
pears natural that the physical pictures developed there 
can also be extended to the short antiferromagnetic 
chains considered here. However, there are indica- 
tions that this expectation is not fulfilled and the situa- 
tion in the antiferromagneti c chains much more involved 
( Konstantinidis et all l201lh . A definitive answer is not 
available at the moment. In the following we will de- 
scribe the molecular horseshoe Cr6, which represents a 
short antiferromagnetic chain with length N = 6, ex- 
hibiting S = ground states. 

The generic spin Hamiltonian for short antiferromag- 
netic chains is 



H 



N-l 



N 



i+1 



gfi B S'B, (87) 



with Si = s for all ions. From molecular symmetry also a 
rhombic term Ei(s^ x — s? ) and a modulation of the 
exchange coupling constant, in particular at the ends of 
the chain, may be present. The presence o f next- nearest- 
neigh bor exchange was also suggested ([Bianchi et all 
20091 ). However, these effects are considerably smaller 
than those due to the leading terms given in Eq. ([87]) . 
Various derivatives of C r^ have been synthesized 



and studied by INS dBaker et all 2011 at lLarsen et all 
l2003bl : lOchsenbein et all 120081 l2QQ7h . He re the cluster 
jNH2R l3[Cr6FnPivio(H20)1 is considered ([Larsen et all 
l2003bf ). It crystallizes in space group P2i/c, and the 
anion forms a string of six Cr 3+ (s = 3/2) ions, see 
Fig. [37l Using INS the magnetic excitat ion spectrum 
up to energies of 5 meV was studied ([Ochsenbein et all , 
I2007D : results are shown in Fig. l42l a). Three cold fea- 
ture I, IV and V are observed, and two hot transitions 
II and III. The data analysis yielded 2 J = —1.27 meV 
and D = in Eq. (|87|) . The model was later refined to 
2J -1.4 meV, 2J edge = -1.1 meV, D -0.028 meV 
and \E\ = 0.005 meV, where J edge ref ers to the coupling 
strengths of the outer coupling paths ([Ochsenbein et all 
I2008D . 

In Fig. 221(b) are shown the simulated energy spectrum 
for Cr6, using the parameters of the simplified model, and 
for comparison the energy spectrum of a hypothetical 
s = 3/2 antiferromagnetic wheel. The L&^-band pic- 
ture is immediately recognized in Cr6, and the observed 
INS transitions demonstrate the L- and i?-band states. 
In comparison to the wheel, however, the E band consists 
of two sub groups with an energy gap A± in the S = 1 
sector. This splitting can be associated to the formation 
of standing spin waves in the chain, as opposed to running 
waves in the wheel. The basic argument is simple, and 
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FIG. 42 (Color online) (a) Neutron-energy loss spectra on 
Cr6 . The solid lines represent simulations) . The inset sketches 
the formation of standing spin waves in a chain and the link 
to the gap in the spin-wave spectrum, (b) Simulated energy 
spectra for Cr6 and a hypothetical iV = 6, s = 3/2 antiferro- 
magnetic wheel. The observed INS transiti ons and the gap in 
the_E-band are indicated. Adapted from lOchsenbein et all 
I2007I . 



familiar from textbooks. In the hexa- nuclear wheel the 
spin waves with shift quantum number q and — g, corre- 
sponding to left and right running waves, are degenerate 
because of cyclic symmetry. However, the open boundary 
in the chain or missing link acts as a perturbation lifting 
the degeneracy, resulting in the formation of the symmet- 
ric and anti-symmetric linear combinations of the wave 
function \q) and | — q) : or standing waves indeed. The 
mechanism is sketched in Fig. l42l a). Hence the splitting 
of the £?-band transitions into peaks IV and V observed 
in experiment directly demonstrates the standing spin 
waves in the Cr6 antiferromagnetic chain. 

These conclus ions are supported by linear SWT 
([Anderson! . Il952h , which can approxi mate the energies of 
the i ? -band states in the S = 1 sector ([Ochsenbein et 
120071 : IWaldmannl . l2002bh . The result for the Cr6 horse- 
shoe and the same hypothetical TV = 6, 5 = 3/2 wheel as 
before is drawn in Fig. l43T a). clearly demonstrating the 
gap in the i?-band spectrum because of the open-chain 
topology. This motivated an analysis of the exact quan- 
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FIG. 43 (Color online) (a) Spin-wave dispersion as obtained 
in linear SWT for Cr6 and a hypothetical Af = 6, s = 3/2 
antiferromagnetic wheel, (b) Exact excitation energies in the 
S = 1 sector of Cr6 and comparison to the results of various 
variants of SWT. LSWT = linear SWT, ISWT = interacting 
SWT, MSWT = modified linear SWT, MISWT = modified 
interacting SWT SLSW T = spin-level SWT. Adapted from 
Ochsenb ein et all \ 2007 to which we also refer for details. 



turn spectrum in Cr6 by different variants of SWTs and 
a newly suggested spin- level SWT, wh ich adds quantum- 
corre ctions to Eq. ([76]) in first-order ( Ochsenbein et aJ] , 
|2007| ). Obviously interacting SWT does best in repro- 
ducing the exact energies, but the accuracy is a modest 
8%. The results of a similar analysis for the CsFeg and 
Feis antiferromagnetic wheels in Sec. IIV.CT2lar e recalled 
( Dreiser et al .l l2010btlUmmethum efa l. 2012). 

However, as mentioned before, the physics of the mag- 
netic excitations in short antiferromagnetic chains, al- 
though displaying a L&^-band structure in the energy 
spectrum, presen ts some subtletie s which are difficult to 
understand ( Konstantinidis et all l201ll ). In the eight 
membered short chain CrgZn, for instance, a detailed 
analysis of the wave functions indi cated a significant m ix- 
ing of the L- and £?-band states (Bia nchi et al. , 2009). 



E. Spin frustration in antiferromagnetic molecular clusters 

In the previous sections clusters with antiferromag- 
netic HDVV interactions on a bipartite lattice were dis- 
cussed, but studying quantum spin-frustration effects in 
large clusters is obviously also of high interest. In fact, 
since the possible geometrical arrangements of metal ions 
and ligand linkages are not restricted by the constraint 
of translational invariance in "zero-dimensional" clusters, 
competing interaction paths are almost always present in 
poly nuclear magnetic molecules, and bipartite magnetic 
molecules are rather the exception than the rule. How- 
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ever, the research on the magnetic excitations in spin- 
frustrated systems concentrated on few model systems. 

Regular antiferromagnetic spin triangles and lattices 
incorporating triangular units are most often considered 
in this context. The HDVV Hamiltonian of a triangle has 
been given in Eq. ([56]) , where for a regular triangle Si = s 
for all ions and J' = J. For s = 1/2, the energy spec- 
trum consists of a doublet of two S = 1/2 spin multiplets 
in the ground state, which transforms according to the 
irrep E of the D3 symmetry group, and a higher- lying 
S = 3/2 multiplet at energy §|2J|. The degeneracy of 
the two 5 = 1/2 multiplets is (often) considered as cri- 
terium for spin frustration. However, small deviations 
such as distortions of the triangle leading to J' ^ J or 
Dzyaloshinski-Moriya interactions lift the degeneracy in 
the ground state opening a gap A (Fig. |44|) . The ratio 
A/ 1 2 J I may here be considered as a figure of merit. 

This structure of low-lying energy levels has at- 
tracted huge interest for a variety of reasons. Spin 
frustration is one of them, but corresponding clus- 
ters are also attractive models for studying the 
Landau-Zener-Stiickelberg (LZS) transitions or d i ssipa- 
tion and decoherence in general ( Chiorescu et all feOQC 
Dobrovitski et all l2QQ0 l : iLandauTl 19321 : iLe 



(a) 



19871 ; IStuckelberdJl932UZenerl , 11932^ , or for applications 
as qubits ([Wernsdorfer et an, 120041), and Rabi oscilla- 
tions were indeed observed ([Bertaina et al l l2008f ). Many 
molecular tri-nuclear clusters are available and were stud- 
ied, also by INS, for their excitations. However, usually 
distortions are strong; the realization of regular triangles 
in bounded clusters is apparently not preferred by nature. 
The V15 molecule is one of the best realizations of a reg- 
ular triangle, according to its A/|2J|, and its magnetic 
excitations are discussed in the next section. 

Topologies with higher nuclearity, which also in- 
corporate triangular units and are of high symmetry, 
such as the cubeoctahedron or icosidodecahedron, 
have also b een intensely stu d ied, h owever, mostly 
by t h eory ( Cepas and Ziman, 2005; Konst antinidis , 



2009: IRousochatzakis et all [20081: ISchnackL 12010 : 
Schroder et al. I. l2005d) . Only few experimental spec- 



troscopic investigations of the cluster excitations are 
available. We will here focus on the magnetic Keplerate 
molecule Fe3o, which has become a main represen- 
tative in this field of research, also because of its 
relationshi p to the Kagome l a ttice ( " Kagome on a 
sphere") (IRousochatzakis et all I2QQ8I : ISchnackL l2010l : 



ISchnack et aZ.L l2QQlh . 

Many more molecules such as tetra- nuclear "butterfly" 
and hepta-nuclear disc-like molecules or odd-membered 
antiferromagnetic wheels and cycles have been synthe- 
sized and studied for spin frustration effects. 



1. The V15 molecule 

The cluster [Vi 5 As 6 O42(H 2 0)]K 6 -8H 2 O, or V15 in 
short, contains 15 V 4+ ions (s = 1/2) whose arrange- 
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FIG. 44 (Color online) (a) Sketch of the energy level scheme 
of an antiferromagnetic spin- 1/2 triangle, (b) Structure of the 
anion of V15. The thick green lines indicate the two hexagons, 
and the triangular area the spin triangle, (c) Sketch of the 
exchange interaction s in the V15 molecule. Adapted from 
iTarantul etaR 120071 and IChaboussant et all l2004al . 



ment can be described as two h exagons sandwiching a 
triangle, see Fig. I44fa) a nd (b) ([Gatteschi et al 
iMiiller and Doringl . Il988h . The system crystallizes in 
space group i?3c, and the molecule exhibits a nom- 
inal crystallographic D3 symmetry, however, a water 
molecule sits at the center of V15. Magnetic suscep- 
tibility demonstrated that the exchange interactions in 
V15 indicated in Fig. l44T b) are antiferromagnetic and 
strong [2Ji « 2J 4 « -13 meV, 2J 2 « 2J 3 « -26 meV, 
2J 5 « —70 meV], such that the hexagons are in a sin- 
glet state at temperatures below ca. 100 K and do not 
contribute magnetic moment. At low temperatures the 
magnetism can hence be described as that of a regular 
spin triangle with an effective antiferromagnetic interac - 
tion of 2 J = -0.211(2) meV (|Chaboussant et all l2002h . 

This picture was confirmed by detailed magnetiza- 
tion m easurements at low temperatures ( Tarantul et all 
l2007h . However, a small gap in the 5 = 1/2 ground- 
state doublet was dete cted. Magneti z ation data indi- 
cated A = 7(2) /ieV (IBarbara et all 120021) a n d low - 
frequency EPR A = 3 /xeV (jKaiivoshia et all l2QQ7h . 
while from INS a gap of A = 35(2) fieV was concluded 
( Chaboussant et all . 12002) . The origin of the gap has 



been discussed controversially, most often it has been as- 
sociated to Dzyaloshinski-Moriya interactions, but also 
to a distortion by e.g. the central water molecule. 

Evidence for the latter came from an INS experiment 
performed on a fully deuterated polycrystalline sample 
of V15 at ult ra- low temperatures and with magnetic 
fields applied (jChaboussant et~al\ , l2004al ) . The field de- 
pendence of the observed INS peaks and their assign- 
ment are presented in Figs. l45l a) and (b), respectively, 
confirming spectroscopically the expected energy spec- 
trum. The INS energies could excellently be fit and 
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FIG. 45 (a) Field dependence of the INS peak energies ob- 
served in V15 at T = 45 mK. Lines are fits to the data, (b) 
Derived energy spectrum and assignment of INS transitions, 
(c) Q dependence of the intensities of peaks III + IV + V in 
zero field, (d) Q dependence of the intensity of transitions III 
+ IV + V at T (circles) and transitio n I at 1 T (squares) . 
Lines are fits to the data. Adapted from IChaboussant et cdi 
l2004al . 




FIG. 46 (Color online) (a) Molecular structure of Fe3o- (b) 
Fe 3+ core forming an icosidodecahedron. The spin orienta- 
tions in the classical ground state are represented by the ar- 
rows, the colors red, green and blue refer to the sublattices 
A, B, and C, respectively, (c) Low-lying energy spectrum as 
function of total spin S as predicted in the rotational-band 
model Eq. ([88]), showing the L and E band, (d) Lowest en- 
ergy in each sector S z as calculated by DMRG (black squares) 
and comparison to the L band in the rotational-band model 
(red line). Both curves essentially superimpose confirming 
the S(S + 1) energy dependence of the lowest eigenstates in 
Fe3o- The inset shows the experimental magnetization curve. 
Pictures: (a , b,d) c ourtesy of J. Schnack, (c) adapted from 
iGarlea et aZll200d . 



yielded 2J = -0.212(2) meV and A = 27(3) /ieV. 
A careful analysis of the magnetic field and Q depen- 
dence of the INS peak intensities, shown in Figs. l45l c) 
and (d), allowed discriminating between Dzyaloshinski- 
Moriya interactions and distortions. It has been con- 
cluded that the gap A comes from a slightly distorted 
triangle, with exchange interactions 2J\2 = —0.21 meV, 
2 J23 = -0.23 meV, and 2 J 13 = -0.20 meV. 



2. The Fe3o Keplerate molecule 

The Fe 30 molecule [Mo 7 2Fe3o 0252(Mo20 7 (H 2 0))2- 
(Mo2O8H2(H2O))(CH 3 COO)i2(H2O)9i]-150H2O con- 
sists of 30 antiferromagnetically coupled Fe 3+ ions 
(s = 5/2) which are located at th e vertices of an 
icosidodecahedron (|Miiller et all Il999h . The molecular 
structure is displayed in Fig. l46Ta) and the iron metal 
core in Fig. HuT b). The spin arrangement clearly 
supports pronounced quantum-spin frustration effects, 
however, the large s = 5/2 s pins also suggest cl assical 
or semi-classical approaches ([Miiller et all l200"lbl ). The 
dimension of the Hilbert space in this molecule is a 
staggering 2.2 xlO 23 and understanding its magnetic 
excitation spectrum is obviously challenging. Fe3o has 



in fact become an ideal test ground for developing 
theoretical schemes and physical concepts. 

Classically, the antiferromagnetic ground state is char- 
acterized by three sublattices A, B , and C , and is highly 
degenerate ([Axenovich and Lubanl l2QQlh . Figure l46Tb) 
depicts one of the possible classical ground-state spin con- 
figurations. For the quantum spectrum a rotational-band 
model with a t h ree su blattice structure was conjectured 
( Schnack et all l200ll ), and based on ideas as those in 
Sec. IIV.B.3l a Hamiltonian was derived, 

2J - ^ \ 

Habc = — — [Sa • Sb + • Sc + Sc • Saj , (88) 

where S^, S#, Sc describe the three sublattices, with 
spin lengths Sa = Sb = Sc — 25. The predicted low- 
lying energy spectrum is presented in Fig. l46T c). and 
shows an L&^-band structure, similar as in the anti- 
ferromagnetic wheels (Sec. HV.C.lj) . However, the degen- 
eracies and spatial symmetry labels are of course differ- 
ent. For the existence of the L-band in Fe3o solid evi- 
dence came from the experimental magnetization curve, 
which increases linearly with magnetic field up to sat- 
uration at a critical field B c = 17.7 T consistent with 
a 5(5 + 1) energy dependence of the lowest states in 
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FIG. 47 (Color online) Neutron-energy loss spectrum od Fe3o 
at T — 65 mK. The solid line represents the estimated back- 
ground, and the such determined magnetic scattering is dis- 
played in the inset. Adapted from lGarlea et a/.Ll2006l . 
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FIG. 48 (Color online) (a) Low-temperature excitation spec- 
trum of Fe3o as predicted by linear SWT (LSWT), modified 
linear SWT (MSWT), and spin-level SWT (SLSWT). (b) Ex- 
perimental INS spectrum of Fe3o (points) and simulated spec- 
trum using SLSWT (line). Inset: Experimental minus simu- 
lated spectrum, providing evidence for magnetic scattering at 
ca. 0.3 meV which the SLSWT model does not account for. 
Adapted from Wal dmannl I2007I . 



each spin sector [inset to Fig. l46Td)] . From the crit- 
ical field the interaction strength has be estimated to 
2 J = —0.134 meV. The L band was also produced in a 
density matrix r enormalization group (DM RG) calcula- 
tion, Fig.HHJd) ([Exler and SchnackL 120031 ) . 

The higher-lying magnetic excitations in Fe3o were 
probed by INS experiments on a deuterated polycrys- 
talline sample. An experimental spectrum is displayed 
in Fig. |47j showing a broad magnetic feature in the en- 
ergy range 0.3-1.1 meV. The INS data could qualitatively 
be interpreted using Eq. (f88j) and associated to the E 
band. Its excitation energy is predicted to (26/5) \2J\ = 
0.67 meV [Fig. [461( c)] in rough agreement with the maxi- 
mum in the neutron scattering intensity (inset to Fig. l47j) . 
The detailed analysis yielded 2J = -0.108 meV. The 
width of the magnetic feature however remained unex- 
plained. 

In order to obtain a better insight the magnetic exci- 
tation spectrum was also calculated using modified lin- 



ear SWT (jCepas and ZimanL 12005). The predicted spec- 
trum is shown in Fig. l48t a). This theory suggests mag- 
netic scattering in the energy range of 3.5 1 2 J|-5.5|2 J| and 
hence accounts at least partially for the observed broad- 
ening. A novel spin-level SWT, which implements a first- 
order quantum correction to Eq. ([88]) . yielded an excita- 
tion spectrum in the range of 3.8|2J| -7.4|2J [Fig.EH(a)], 
and allowed to reproduce most of the observed magnetic 
scattering using 2J = —0.125 meV [Fig. l48lf b)]. consis- 
tent t o within 7% with the finding from the magnetic 
data ( Waldmannl . 120071 ). However, a peak- like magnetic 
scattering at ca. 0.3 meV is not reproduced. 

It was pointed out that these SWT techniques, when 
applied to three- (a nd more) pa r tite s ystems, have con- 
ceptual drawbacks (Wal dmannl . l2QQ7h . The spin- leve l 
SWT was later extended to higher-orders (jSchnalle et all . 
2009), further emphasizing this point. A fully satisfying 
underst anding of the magne tic excitations in Fe3o is still 
lacking (jGarlea et aZ.L l2QQ6h . 



V. SINGLE-MOLECULE MAGNETS 

A. Introduction 

The single-molecule magnets (SMMs) are a special 
subclass of the molecular nanomagnets. Some promi- 
nent examples are shown in Fig. [49l The SMMs are 
distinguished by exhibiting slow relaxation of the mag- 
netization or magnetic hysteresis at low temperatures, 
below a blocking temperature Tg. This phenomenon is 
not due to a long-range ordered magnetic ground state 
as in conventional magnets, but arises from an energy 
barrier for spin reversal at the molecular level. Fur- 
thermore, quantum phenomena such as quantum tun- 
neling of the magnetization are observed as character- 
istic steps in the magnetic hysteresis curves. These 
unique magnetic properties discovered about 15 years 
ago stimulated enormous research and allowed address- 
ing fundamental questions in quantum mechanics as 
well as suggested applications in information technol- 
ogy as classical or quantum bits or in molecular spin- 
tronics. Many authoritat ive reviews and book s were 
written, some o f them are Barbara and Guntherl . 
del Barco et all . 120051 : iBogani and Wernsdorferl . 
Christou et all l2000l: fchudnqvskv and Teiad; 



1999; 



2008; 



1998; 



Friedmann and Sarachikl, l2010l: iGatteschi anc Sessoli 



120031 : IGatteschi eTall . l2006bl : iLeuenberger et all . 120031 . 

The magnetic spectrum in the SMMs is in princi- 
ple also described by the microscopic spin Hamiltonian 
Eq. ([70|) . However, the most fascinating SMM phenom- 
ena such as quantum tunneling of the magnetization are 
related to the ground- state multiplet, and the anisotropy 
splittings in it. The effective spin Hamiltonian, which 
is known as "rigid-spin" or "giant-spin" or "single-spin" 
model, reads then 



H = H z + H± 



(89) 
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FIG. 49 (Color online) Molecular structures of three promi- 
nent SMMs, (a) Mm 2 : [Mni20i2(CH 3 COO)i6(H20)4], (c) 
Fe 8 : [Fe 8 02(OH)i2(tacn)6]Br8-9H 2 0, (d) Mn 6 : [Mn 6 2 (Et- 
sao) 6 (02CPh(Me) 2 )2(EtOH) 6 ]. Panel (b) shows the metal 
core in the Mni2 molecule, the exchange coupling paths, and 
the classical ground state spin configuration. 

= DS 2 Z + B° 4 d1(S) 

+E (S 2 X - sfj + Bldj(S) + BjOj(S), (90) 

where S is the total spin of the ground-state multi- 
plet. The terms were distinguished into H z and H± ac- 
cording to whether they commute with S z or not, and 
only the usually most relevant higher-order terms were 
listed [terms up to 6 th order were demons trated in high- 
precision experiments ([Barra et all 120071 ) , compare also 
with Sec. IIII.C.5] . In a magnetic field the Zeeman terms 
I^b9zS z B z and I~ib{9xS x B x + g y S y B y ) have to be added 
to H z and H±, respectively. In principle the giant-spin 
model can be derived perturbatively from Eq. ([TO]) . 

In a SMM the uniaxial anisotropy, which is the dom- 
inant zero-field splitting term in Eq. ([89]) . has to be 
of easy-axis type or D < 0. The energy levels \M) 
of the ground-state spin multiplet organize then into a 
parabolic band according to DM 2 , with the M = ±5 
levels being lowest and separated by an energy barrier of 
height U = \D\S 2 for integer and U = \D\(S 2 - 1/4) for 
half-integer S (see inset to Fig [50]). This energy barrier 
is responsible for the slow relaxation of the magnetiza- 
tion. The non-commuting terms contained in H± induce 
a mixing between the \M) levels, or quantum tunneling 
indeed. 

As regards spectroscopy, the major goal is the precise 
determination of the parameters in the giant-spin Hamil- 
tonian, where in particular the non-commuting terms in 
H± are of most interest, because of their relation to the 
quantum tunneling rates. As the relevant transitions are 
intra-multiplet with AS = many experimental tech- 




Energy transfer (meV) 

FIG. 50 (Color online) INS spectrum of Fe 8 at 9.6 K. Circles 
represent the experimental data, dashed lines Gaussian fits, 
and the solid line a simulation. The transitions are labeled 
according to | zb M) — > | db M'). The inset shows the parabolic 
energy spectrum with th e allowed INS transit ions indicated 
by arrows. Adapted from ICaciuffo et all , 1998 



niques can accomplish the task. EPR is one of them, 
and high-frequency hi gh-field EPR has indee d played a 
most important role (Ga tteschi et all l2006al ). INS has 
also been very valuable, adding the advantage of a zero- 
field spectroscopy. 

Figure [50] shows an INS spectrum recorded 
on a non-deuterated sample of the compound 
[Fe 8 02-(OH)i 2 ( t acn) 6 lBr8-9H 2 0, or Fe 8 henceforth 
( CaciufFo et all [1998). The molecule contains eight 
Fe 3+ ions (s = 5/2) in a butterfly- like arrangement, 
Fig. I49t c). It crystallizes in space grou p PI and exhibits 
appro ximate D 2 molecular symmetry ( Wiegh ardt et all 
1984). The frustrated AFM interactions in the molecule 
result in a S = 10 ground state. In the experimental 
INS data in Fig. [50] the allowed AM = ±1 transitions 
are observed, demonstrating the zero-field splitting in 
this spin multiplet. The "picket-fence" INS spectrum, 
which is typical for SMMs, results from the fact that 
the transition energies given by the dominant D term in 
Eq. dBU) vary as \D\(2M - 1) or linearly in M. Slight 
deviations from the regular pattern are however seen, 
which relate to the other terms in Eq. ([89]) . and precise 
values for its five p arameters could be determined 
( CaciufFo et all [l998). Besides the scientific impact of 
this work, it also showed that excellent INS data can 
be recorded on non-deuterated samples of molecular 
nanomagnets, which established the basis for many of 
the subsequent INS works in this area. 

Most of the SMMs which have been synthesized so far 
are based on 3d metal ions, and several of them were 
studied by INS, but mostly for the trans itions within 
the g r ound-state multiplet (for reviews see lBasler et all 
120031 : iBircher et all . l2006h . However, being related to 
anisotropy and not the exchange splittings, these stud- 
ies shall not be further discussed here; Feg may serve 
as a representative example. The potential of INS to 
also detect exchange splittings has been applied to only 
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very few SMMs such as Fe 4 (jCarretta et gfll2004bh , Fe 8 
( Carretta et g/.L l2006), Mni2, and Mn^, and derivates of 
them. In the following the situation in Mni2 and Mn6 
will be presented. 

It is added that in the last few years the focus shifted 
from 3d-based SMMs to 3d-4f or 4f SMMs (also 5d 
ions became of interest). Since 4f ions are type J and 
bring in significant anisotropy, while exchange interac- 
tions are weak, the analysis of experimental data is typ- 
ically more involved than in 3d-based type S/Q clusters. 
INS work on lanthanide containing SMMs has jus t be- 
gun ([Preiser et all 12012k iKlokishner et all l2QQ9bh . but 
promises exiting results in the future. 



B. The M1112 cluster 



The compound [MmsO^CHsCOOMHsO^], or 
Mni2 in short, was the first molecule for which SMM 
behavior and quantum tunneling of the magnetization 
was observed bel o w a blo cking temperature of ca. 3.5 K 
([Gatteschi et all l2006bf h and became the prototype 
SMM. The system crystallizes in space grou p J4, and the 
molecule exhibits 64 symmetry ([LisL |l980). It contains 
an inner tetrahedral core of four Mn 4+ ions (s — 3/2) and 
an outer ring of eight Mn 3+ ions (s = 2), Fig.l49la). The 
nearest-neighbor Heisenberg exchange interactions in the 
cluster result in a total spin S = 10 ground state, whose 
classical spin structure can be depicted as displayed 
in Fig. l49f b). The molecule's properties in its S = 10 
round-state mu l tiplet have been heavil y investigated 



Gatteschi etaU. l2QQ6bh, als o bv INS dBircher et al 



2004: Hennion et 



Sieber et al 
1999). 



19971: iMirebeau et all 11222 ; 



2005; Waldmann et al, 2006a; Zh ong et al 



However, understanding how the competing Heisen- 
berg interactions in the molecule give rise to the S = 10 
ground state is also of high interest, since through the S- 
mixing effects induced by the anisotropy the higher-lying 
spin multiplets can also significantly affect the S = 10 
ground-state multiple t and tu nneling rates ([Barra et all 
I2OO7I : ICarretta et all l2004af ). Several attempts were 
made to infer the values of the exchange constant Ji, 
J2, J3 a nd J4 [see Fig. B9Tb)1 with con t roversial con- 
clusions ( Hartmann-Boutron et all 1 19961: iRaghu et all 
120011 : iRegnault et all 120021 : ISessoli et a/.LIl993h 

The issue was targeted by high-energy INS 
experiments on a deut erated sample of Mni2 
( Chaboussant et~al\ , l2004bl ). INS spectra are shown in 
Fig.[5TJ A peak at ca. 1.2 meV is found, which originates 
from the | ± 10) —> \ ± 9) transition within t he S = 10 
groun d- state multiplet [D = —0.057 meV ( Bircher et all 
l2004f )]. Besides that, several inter-multiplet transitions 
reflecting the exchange splittings were observed in the 
energy range up 35 meV. From a detailed analysis of 
the experimental data the energy level scheme presented 
in Fig. [52] was derived. Interestingly, the anisotropy 
splittings in Mni2 are not significantly smaller than the 
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FIG. 51 Neutron-energy loss spectra of Mni2 recorded in 
two energy regimes. Data were recorded at the direct TOF 
spectrometer MARI at Rutherford Applet on Laboratory ISIS. 
Adapted from Chaboussant et al 2004& 



exchange splittings, as expected in the strong-exchange 
limit, which can be seen in Fig. l52l a) from the fact, that 
the lowest level of the lowest S = 9 multiplet falls below 
the top of the S = 10 ground- state multiplet. Hence, S 
mixing plays a relevant role, which was later confirmed in 
a high -precision EPR exp eriment on the S = 10 ground 
state (|Barra et all 120071 ). Using numerical large-scale 
calculations, the exchange coupling constants were 
determined to 2J\ = —5.79 meV, 2J2 = —5.33 meV, 
2 J 3 = -0.67 meV and 2 J 4 = -0.48 meV. 

The exchange couplings J3 and J4 are one order of 
magnitude smaller than J\ and J2, and spin frustration 
is hence small in the Mni2 molecule, which justifies the 
interpretation of the S = 10 ground state in terms of the 
classical spin configuration shown in Fig. H9Tb) and sug- 
gests a topological relation to AFM wheels ([Waldmannl 
l2005af ). In the context of Sec. IIV.C.1( this implies that 
the L&£^-band concept is also realized in the Mni2, which 
was underlined by a succe ssful linear SWT analysis o f the 
INS excitation spectrum ( Chaboussant et ~al\ . l2004bh . 



C. The Mn 6 clusters 

The SMM [Mn 6 02(Etsao)6(0 2 CPh(Me)2)2(EtOH) 6 ], 
or Mn6 henceforth, which has been synthesized 14 years 
after the initial di scovery of SMM behavior in Mni2 
(|Sessoli et aZ.L Il993h . was yet the first molecule found to 
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FIG. 52 Energy level scheme as derived from the experimen- 
tal INS data. Panel (a) shows a zoom into the lo w-energy 
spectrum. Adapted from IChaboussant et all l2004bl . 



exhibit a higher ani sotropy barrier than that in Mni2 
([Milios et all l2QQ7bh . It contains six Mn 3+ ions (s = 2), 
which are arranged into two triangles linked by oxygen 
atoms, see Fig. l49f d). It crystallizes in space group P2-| /n 
and th e molecule exhibits C2/7, symmetry (|Milios et all 
l2QQ7b[ ) . Several derivates of this molecule have been syn- 
thesized and a magneto-struc tural correlat i on for the ex- 
change couplings established (Mil ios et all l2007at ). 

Magnetic measurements demonstrated a S = 12 
ground state due to overall FM exchange interactions 
in the cluster. The zero-field splitting parameter in the 
S = 12 ground state was estimated to D « —0.054 meV, 
corresponding to an energy barrier of \D\S 2 ~ 7.7 meV. 
The exchange couplings were however determined to 
2 J « 0.4 meV, hence suggesting that in Mn6 S mixing is 
very strong and the strong-exchange limit or giant-spin 
model break down. 

This has been indeed convincingly confirmed in an 
INS experimen t on a non-deuterated sample of Mng 
dCarretta et a/.l , 12008). Experimental INS data recorded 



in two energy windows are presented in Fig. l53Ta). At 
low energies, at 1.1 meV and below, the "picket-fence" 
spectrum characteristic for the transitions within the 
ground-state multiplet of a SMM is observed, see inset to 
Fig. I53lfa). However, at already slightly higher energies 
(1.9 meV) the first inter- multiplet transition appear, fol- 
lowed by several more in the experimental energy range 
of 6 meV. This directly points to the fact, that several 
higher-lying spin multiplets are partially nested with the 
ground-state multiplet. 
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FIG. 53 (Color online) (a) Neutron-energy loss spectra in 
Mn6 recorded for two energy regimes up to 6 meV at the 
indicated temperatures (symbols). The lines represent simu- 
lations using Eq. (|91|) and the parameters given in the text, 
(b) Calculated energy spectrum as function of M. The color 
repres ents the value of the expectation value (S 2 ). Adapted 
from ICarretta etaR I2008I . 




FIG. 54 (Color online) The core of Mn6 with the labeling 
of the Mn 3+ s ites and the assumed coupling paths indicated. 
Adapted from Carretta et ol 2001. 



The data analysis was based on the Hamiltonian 
H = -2 > ^ J,,s z • s ? + > ' | Disl + BlOKsi)} , (91) 



i<j i 



where three exchange coupling paths Ji, J2, J3 were 
assumed as displayed in Fig. [54j The ligand-field pa- 
rameters for sites 1 and 1', 2 and 2', and 3 and 3' are 
identical by symmetry, those of sites 2, 3, 2', 3' were as- 
sumed to be equal, and B\ x jB\ 2 = D1/D2 employed. 
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The experimental data could be excellently reproduced 
with the best-fit parameters 2J X = 0.84(5) meV, 2J 2 = 
0.59(3) meV, 2J 3 = -0.01(1) meV, D 1 = -0.20(1) meV, 
D 2 = -0.76(2) meV, and B^ = -0.0010(3) meV. 
The resulting simulated energy spectrum is shown in 
Fig. l53Tb). which demonstrates that the energies of the 
\M) states of the S = 12 ground-state multiplet do 
not follow the generic parabolic M 2 curve but deviate 
strongly from it at higher energies, and that indeed sev- 
eral higher- lying S = 11 m ultiplets fall below t he top 
of the S — 12 multiplet. In lCarretta et al 1. 120081 it was 
further demonstrated that this unique structure of the 
energy spectrum has significant effects on the magnetic 
relaxation rates. 



VI. QUANTUM SPIN SYSTEMS 

A. Introduction 

Quantum spin systems have been attracting much at- 
tention due to numerous magnetic features w hich can- 
not b e interpreted by conventional spin models ( Sachdevl 
2008). In particular, classical magnetic phases such as 
ferromagnetism and Neel antiferromagnetism are pre- 
vented by (strong) quantum fluctuations, which are 
present in magnetic compounds if some of the following 
conditions are fulfilled: 

(i) The spin quantum number of the magnetic ions is 
low, i.e., Si = 1/2 or 1. 

(ii) The dimensionality of the magnetic system is low, 
i.e., d = 1 or 2. 

(hi) The connectivity of the network of magnetic ions 
(i.e., the number of spins to which each spin is coupled) 
is low. 

(iv) The interactions between the magnetic ions are 
geometrically frustrated. 

Representatives of the categories (i) and (ii) have 
been the subject of detailed investigations for a 
long time. The y include one- d i mensio n al magnets 
such as KCUF3 ( Hutchings et all Il979l ; lLake et a/1 
12005 ) and two-dimensional magnets such as La2Cu04 
(IColdea et a/]j200ll : [Headings et a/I . ^OlOtlVaknin et al 



(jKagevama et all l2000l: iKodama et all l2QQ2h 



1987), the latter being of tremendous interest due to 
the observation of doping-induced superconductivity 
([Bednorz and Mulier! 1 19861 ). More recently, novel 
materials have been synthesized which are formed by 
two and three magnetically coupled spi n chains, e.g. , 
two- l eg spin ladders such as SrC u9Q^ (lAzuma et al , 
19941) and fC^H 12 N) 2 CuBr 4 (jRiiegg et all 1 2008: 



Thielemann et al , 2009) as well as three-leg spin ladders 
such as Sr2Cu305 (jAzuma et o/lll994l ). Typical realiza- 
tions of the category (iii) are weakl y coupled dime r- based 
compounds i nclu ding KCUCI3 (Cavadi ni et all 
Kato et all Il998h . TlCuCl 3 (Oos awa et a/l 



Riieg g et al 



1999; 



Riieg g et al 



Riiegg et al 



20031). NHaCuCU (jKurmawan et al 
, 2004bh. BaCuSi 2 Q6 (I Jaime et at 
, 120071 : ISebastian et all [2005|), SrCu 2 (B0 3 ) 



1999; 



2002; 



2004; 



Cs 8 C r 2 Br 9 (jGrenier etaU, I2QQ4J: ~ Leuenberger et al, 
I l985h. S rc. Cr 2 Q8 (|Quintero-Castro et all I201C : 



and 
l2002h. 



Ba3Mn 2 0s ([Stone et al 



Trimer-b ased compounds 



Wang et all l2Qllh 

120081 : lUchida et att _ , . 

such as LaACuc»MoOi9 \0'm et a/1 l2QQ5h as well as the 
compound SrCu 2 (B03) 2 (characterized by an array of 
mutually perpendicular dimers) fulfill in addition the 
criterion (iv). 

According to the topic of the present work, we will 
consider here only compounds for which the presence 
of magnetic clusters is the most important ingredient 
to understand their quantum spin properties. The fol- 
lowing sections will therefore focus firstly on weakly 
interacting AFM dimer systems, for which quantum 
phase transitions involving gapless excitations and Bose- 
Einstein c ondensation hav e been observed most con- 
vincingly ( Giamarchi et all [2008). Secondly, the phe- 
nomenon of a spin-Peierls transition will be described 
that occurs in quasi one-dimensional (ID) antiferromag- 
nets due to the formation of spin pairs as a result 
of dimeriza t ion o f the regular array of magnetic ions 
(|Brav et all Il975l ) . Thirdly, the formation of magnetic 
polarons observed in transition-metal perovskites will be 
discussed which evolve upon hole doping and behave like 
magnetic clus t ers e mbedded in a non-magnetic matrix 
(iPhelan et a/.Ll2006h . 



B. Dimer-based antiferromagnets 

In dimer-based compounds the two magnetic ions 
are antiferromagnetically coupled according to the spin 
Hamiltonian Eq. ([44]) . In contrast to isolated dimer 
systems discussed in Sec. IIII1 the coupling between the 
dimers cannot be neglected, but the intra-dimer ex- 
change coupling Jo is larger than the coupling between 
the dimers J\ = Jo/7- As long as 7 ^> 1, the thermo- 
dynamic magnetic properties resemble those of isolated 
dimer systems, i.e., the dimer ground state \S, M) is a 
singlet |0,0), and so preserves full rotational invariance, 
unlike the Neel state of a square lattice antiferromagnet 
with 7 = 1. The singlet state is separated from the ex- 
cited triplet states |1, +1), |1, 0) and |1, — 1) by an energy 
gap A = 2 1 Jo I . This picture is confirmed by the magnetic 
susceptibility measu red for the dimer compound KCUCI3 
( Tanaka et a/.L[l997l ) as shown in Fig.l55la), which read- 
ily allows the determination of the gap energy, similar to 
the analysis of Fig. |2) However, the high-field magneti- 
zation data ([Oosawa et al l l2002h displayed in Fig.ESfb) 
are basically different from the case observed for iso- 
lated cluster systems. The latter exhibit sharp step- 
like enhancements of the magnetization (see Figs. [TT1 and 
[22]) , which are absent for KCUCI3. Instead, the mag- 
netization increases rather slowly above the critical field 
B c ^ 23 T, and the saturation moment of 1 /i#/Cu 2+ 
is reached at the saturation field B s w 53 T far above 
B c . This is due to the inter-dimer coupling Ji giving 
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FIG. 55 (a) Magnetic susceptibility measured for KCUCI3 
with B = 0.5 T oriente d perpendicular to th e cleavage plane 
(1,0,-2). Adapted from iTanaka et all Il997l . (b) Magnetiza- 
tion measured for KCUCI3 at T = 1.7 K with B perpendicular 
to the cleavage plane. The saturation field B s » 60 T was 
outsi de the available field range. Adapted from lOosawa et "all 
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FIG. 56 (Color online) Selected results obtained from INS 
experiments on TICUCI3. (a) Field dependence of the mag- 
netic excitation energies measured at Q = (0,4,0). The 
solid lines reflect a linear Z eeman model. The critical field 
is B c = 5.7 T. Adapted from lRiiegg et q/ll2Q0A (b) Splitting 
of the singlet-triplet excitation measured at Q = (—0.5, 0, 2) 
and T — 1.5 K. The asymme tric line shapes are ty pical reso- 
lution effects. Adapted from ICavadini et al. 2002. 



rise to energy dispersion of the singlet-triplet excitations 
which may be called triplons. The triplon dispersion 
can be calculated in the ra ndom phase approximation 
( Jensen and Mclntoshl Il99lh , since perturbation theory 
in I/7 is well defined. For the case of dimers forming a 
square lattice we find 

E M (k) = -2 J - 2 Ji [cos(k x d) + cos(k y d)] + q\ibBM, 

(92) 



where k = (k x ,k y ) is the quasiparticle wave vector, d 
the lattice constant, B an external magnetic field, and 
M = 0, d=l. In zero field the three triplons are degener- 
ate. With increasing field the triplons are split into three 
branches with M = +1,0,-1 as shown in Fig. l56t a). 
At a critical magnetic field B C1 the energy of the low- 
est triplet component intersects the ground-state 
singlet |0,0) and the ground state changes, thus B c is 
a quantum critical point separating a gapped spin-liquid 
state (B < B c ) from a field-induced magnetically ordered 
state (B > B c ). The triplet components |1,+1) can be 
regarded as bosons with hard-core repulsion, thus Bose- 
Einstein condensation occurs at the quantum critical 
point B C1 i.e., the gas of triplet bosons undergoes a phase 
transition into a novel condensate state with macroscopic 
occupation of the single-particle ground-state. 

The state of an individual dimer at position r is well 
approximated by a linear combination of the singlet |0, 0) 
and the triplet |1, +1), 

|*> = iz|0, 0) - vexp[z(k • r - E +1 t)]\l, +1), (93) 

where E+i is defined by Eq. ([92]) , an d the amplitudes u 
and v depe nd on the magnetic field ([Matsumoto et all 
20021 \2004 . The expectation values of the spin operator 
components of the dimer in the state Eq. ([93]) are 



<Si,*> = 



-(S2, y ) oc — cos(k-r -£+!*), (94) 



(S2,z) OC —. 



(96) 



uv 



~(s2,x) oc — sin(k • r - E +1 t), (95) 



The condensate at B > B c can be associated with the 
transverse order parameters (si, x ) and (si j2/ ) which are 
oppositely aligned at the dimer sites 1 and 2 due to the 
AFM dimer coupling Jo. The rotational symmetry 0(2) 
of the underlying spin Hamiltonian is spontaneously bro- 
ken for B > B c , giving rise to a dramatic change in the 
nature of the magnetic excitation spectrum. According 
to Eq. ([92]) , the mode associated with the lowest triplet 
|1,+1) exhibits a quadratic dispersion around the wave 
vector k = (ir/a,ir/a) and becomes gapless at the crit- 
ical field B c , but for B > B c it transforms in a fully 
isotropic system into a gapless Goldstone mode with a 
sound-like, linear dispersion 

EQl) = hs\k-k \ (97) 

where s is a velocity (Mat sumoto et all 120021 . l2004f ). 

One of the most widely studied dimer-based antiferro- 
magnets is the monoclinic compound TICUCI3 in which 
the Cu 2+ ions are arranged in centro-symmetric pairs. 
The intra-dimer coupling Jo = —2.7 meV dominates the 
inter-dimer couplings |Ji| < 1 meV. INS experiments 
confirmed the singlet-triplet nature of the magnetic ex- 
citations by apply ing an external mag netic field as visu- 
alized in Fig. Eg] dRiiegg et all l2003l i. For B > the 
singlet-triplet excitation is clearly split into three lines 
due to the Zeeman effect, with the central AM = 
line being twice as intense as the AM = ±1 side lines 
as predicted by Eq. JT5J). The triplons are well char- 
acterized by a three-dimensional extension of Eq. ([92]) 
with a quadratic dispersion around Q = (0,0,1) as 
shown in Fig. l57f a). The spin gap has a minimum value 
A « 0.70 meV at the zone center (0, 0, 1). With increas- 
ing field strength, the lowest triplet state |1, +1) reduces 
its energy continuously and overcomes the spin gap for 
a critical field B c = 5.7 T, see Fig. l56t a). Since a fi- 
nite number of triplet states |1, +1) is created at B C1 the 
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FIG. 57 (Color online) (a) Energy dispersion of the triplons 
observed in TIC11CI3 at T = 1.5 K and B = 5.5 T. 
The lines correspond to model expectations based on a 
three-dimensional ex tension of Eq. (|92[). Adapted from 
ICavadini et all 120021 . (b) Energy dispersion of the low-lying 
magnetic excitations observed in TIC11CI3 at different temper- 
atures, fields and pressures. The lines denote the linear be- 
havi or of the Goldstone mode a ccording; to Eq. (|97|) . Adapted 
from lRiiegg et q/.L I2Q03L l2004bL 



system undergoes a phase transition to a magnetically or- 
dered state. This condensation of triplet states at T « 
is therefore a prototype of a quantum phase transition. 

What is the experimental proof that TICUCI3 is Bose- 
Einstein condensed at B c ? Important evidence is pro- 
vided by the critical exponent <j) in the field dependence 
of the critical temperature, T C {B) oc (B c — B)v. The the- 
ory for a three-dimension al Bose gas predicts a universa l 
value (j) = 2/3 (see, e.g. iGiamarchi and Tsevlikl fl999). 
which was experimental ly confirmed , e.g., for the dimer 
compounds BaCuSi^ Oa ([Jaime et aZll2QQ4h and TlCuCl 3 
( Tanaka et all l2QQ7h in some temperature range. The 
ultimate proof for the Bose-Einstein condensation in 
TICUCI3, however, is offered by the properties of the mag- 
netic excitation spectrum above B c according to Eq. ([97)) . 
This theoretical prediction was observed for the first time 
by INS experiment s in TICUCI3 as shown in Fig. f57Tb) 
([Riiegg et all 120031 ) , thus the presence of a spin- wavelike 
mode with a linear dispersion is a convincing signal for 
the existence of the Bose-Einstein condensate. 

The singlet-triplet gap of TICUCI3 can also be closed 
by the application of hydrostatic pressure (thereby mod- 
ifying the parameter 7) which occurs at a critical pres- 
sure p c « 0.1 GPa. In contrast to the field induced 
case, all the triplet components can condense into the sin- 
glet ground state at the quantum critical point p c . The 
magnetic excitation spec trum has again th e nature of a 
gapless Goldstone mode ([Matsumoto et a/1 l2004f ) which 
was experimentally confirmed in th e pressure-induced or - 
dered phase as shown in Fig. l57T b) ([Riiegg et a/l l2004ah . 
In later INS experiments, it was demonstrated that only 
the longitudinal and one transverse triplet component 
soften at p C: whereas the other tr ansverse tripl e t com - 
ponent retains a finite gap at p c ([Riiegg et all l2Q08h . 
For p > Pc, the gap energies of both transverse compo- 
nents remain constant, whereas that of the longitudinal 
component gradually increases. The data could consis- 



tently be interpreted by using a mainly linear pressure 
dependence of the exchange parameters as well as a tiny 
exchange anisotropy. The gap energy of the longitudi- 
nal mode increases with the ordered magnetic moment 
above p c , with a fundamental ratio of y/2 between the 
gaps in the ordered and disordered states, thereby provid- 
ing a non - trivia l experimental test of the field theory 
([Sachdevl l201ll ). Such an amplitude mode is not present 
in a classical description of ordered magnets, but is a 
direct consequence of the underlying quantum criticality. 

Up to the present, the compound TICUCI3 has re- 
mained a prototype of a quantum antiferromagnet in 
which evidence for the Bose-Einstein condensation was 
given by many techniques. The concept of a Bose- 
Einstein condensation has been applied to some other 
dimer based compounds mentioned in Sec. PVT. Al as well, 
but often factors such as the large values of the exchange 
parameters as well as the presence of anisotropics violat- 
ing the rotational symmetry prevent the complete soft- 
ening of the lowest triplet component |1, +1) at B c . The 
influence of anisotropics, which can be d etermined rather 
preci sely by EPR techniques (s ee, e.g., iKolezhuk et all 
2004), was discussed elsewhere ([Giamarchi et all 12008). 

Dimerized antiferromagnetic chain systems are closely 
related to the ACUCI3 (A = K, Tl) compounds dis- 
cussed above. A well studied example is copper nitrate, 
Cu(N0 3 ) 2 -2.5(D 2 0), in which the spin- 1/2 Cu 2+ ions 
are arranged as chains of copper pairs. Each pair has 
a singlet ground state separated from a triplet at |2Jo| = 
0.44 meV. The weak interdimer coupling J\ = J0/7 with 
7^4 yields triplet excitations tha t propagate coher- 
ently along the chain (|Xu et al. ,2000). The experimental 
data displayed in Fig. [58] are well described by the one- 
dimensional variant of Eq. ([92]) with k = (^,0) where 
x denotes the chain direction. Above the singlet-triplet 
gap there is a second gap to the multimagnon contin- 
uum. Two-magno n bound states with S = 1 are visible in 
INS experiments ([Tennant et al., 120121: iTennnant et all 
2003 ) and have a dispe rsion (Schmid t and Uhrid . 12003: 



Uhrig and SchuTzl . Il996h 



E BS (k) = -2 J - ^[1 + 4cos 2 (M)], (98) 

as shown in Fig. [58] These 5 = 1 bound states exist 
only over the range \nir — k x d\ ^ 7r/3, where n is an odd 
integer. 



C. Spin-Peierls dimerization 

For a compound built up by identical atoms the elastic 
energy is lowest if the atoms are equally spaced. How- 
ever, there are compounds where the atoms move from 
an equally-spaced crystal to one in which the spacing al- 
ternates, i.e., the atoms form pairs. This is called dimer- 
ization (first proposed for one-dimensional systems by 
Peierls in the 1930s in a textbook entitled "Quantum 
theory of solids" ) , which is made possible through lower- 
ing of the free energy of the electronic subsystem by this 
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FIG. 58 (Color online) Excitation spectra of copper nitrate at 
T = 0.12 K. Panels (a) and (c) show background-subtracted 
two- and one-magnon INS data, respectively, while panels (b) 
and (d) show the s imulated T = spectra. Adapted from 
iTennant "oil . 12011 



maneuver. The early examples mainly included polymer- 
like organic materials, characterized by antiferromagnetic 
Heisenberg chains, being fir st discovered in t he com- 
pound TTF-CuS 4 C4(CF 3 )4 (iBrav et all Il975h an d af- 
terwards in MEM-f TCNQ^ (iHuizinga et all J l979l) and 
in DEM-(TCNQ) 2 (|Schwerdtfeger et qZlll982h . Here we 
exemplify the phenomenon of a spin- Peierls transi t ion fo r 
the linear-chain compound CuGeOs (|Nishi et all Il994h . 
followed by a discussion of the three-dimensional com- 
pound CeRu2Alio, in wh ich recently a dimer ization has 
been suggested to occur (Ro bert et al l l2Q10h , 

CuGeOs is a linear Cu 2+ (si = 1/2) chain compound 
crystallizing in the orthorhombic space group Pbmm 
with a = 4.81 A, b = 8.47 A and c = 2.941 A at room 
temperature (|Vollenkle et all Il967l ) . The coupling of the 
Cu 2+ ions positioned at (1/2,0,0) is strong along the c 
axis giving rise to physical properties typical of a one- 
dimensional magnet. The magn etic suscept i bility rapidly 
drops to zero below T c = 14 K (|Hase et a/.Ul993l ), which 
points to the opening of a finite energy gap associated 
with a singlet ground state. This can be understood in 
terms of lattice dimerization, i.e., the formation of Cu 2+ 
pairs resulting in a spin-Peierls ground state. Neutron 
diffraction experiments indeed gave evidence for Cu-Cu 
dimerization along the c axis below T c with an inter- 
atomic separation of 2.926 A, compared to 2.941 A in 
the high-temperature structure ( Hirota et a/.L ll994), ac- 
companied by shifts in the position of the 0(2) ions in the 
(a, b) plane as illustrated in Fig. [59j Since the coupling 
between two Cu 2+ ions in the c direction is mainly due to 
a superexchange interaction through 0(2), the dimeriza- 
tion is clearly driven by the 0(2) shifts, which results in 
two unequal and alternating exchange parameters along 
the c direction: 



Ji,2(T) = J c [l±5(T)]. 



(99) 
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FIG. 59 (Color online) Schematic representation of the struc- 
ture of CuGeOs. The arrows denote the AFM spin alignment 
of the Cu 2+ ions. The dashed lines mark the Cu-0(2)-Cu 
superexchange bonds, (a) T > T c = 14 K. (b) Copper dimer- 
ization for T < T c indicated by the dashed boxes. 
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FIG. 60 (a) Energy scan profiles of the magnetic excitation 
observed for CuGeOs at Q = (0, 1,0.52) for various temper- 
atures. The full and dashed lines denote Gaussian fits to the 
data obtained at T — 4 and 12 K, respectively, (b) Dispersion 
of the magnetic excitation observed in CuGe03 alonff the c 
direction at T — 4 K. Adapted from Nishi et all 11994 



The existence of the spin-Peierls gap wa s verified by 
INS experiments as shown in Fig. l60Ta) (|Nishi et all 
Il994f ). At T = 4 K a sharp peak corresponding to the 
singlet-triplet transition appears at an energy transfer 
of 2.5 meV for Q = (0,1,0.52) which is close to the 
Brillouin zone center (0,1,0.5). With increasing tem- 
perature the peak moves to lower energies and broad- 
ens (T = 12 K), while at T = T c = 14 K the energy 
gap vanishes as expected. The dispersion of the singlet- 
triplet transition is shown in Fig. l60f b). The energy at 
the zone boundary (z = 1/2) yields directly the intra- 
chain exchange parameter J c = 10.4 meV from the for- 
mula 7rJ c /2 = 16.3 meV derived by Pes C loizeaux and 
Pearson (jDes Cloizeaux and Pearson! Il962l ) . This infor- 
mation is useful to derive the value of S(T) in Eq. ([99]) . 
The mean-field theory gives 



S(T) 



ACQ 

PJc 



(100) 



with p = 1.637 (|Brav et aLl . ll975|) . Substituting A(0) = 
2.1 meV by extrapolation of the gap energy A(4iT) = 
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FIG. 61 (Color online) Schematic representation of the Ce- 
Al zig-zag chains in CeRu2Alio. The Ce 3+ dimerization for 
T < To is indicated by the dashed boxes. 



2.5 meV [see Fig. lfjUT a)] to zero temperature yields 6(0) = 
0.12 and Ji(0)/J 2 (0) = 1.27 from Eq. flggj). 

Recently the ternary compound CeRu2Alio has at- 
tracted much attention because of a phase transition tak- 
ing place at Tp = 27 K whose origin remained unclear 
initially ( Strvdoml l2009f ). CeRu2Alio crystallizes in the 
orthorhombic space group Cmcm in which the Ce 3+ ions 
are separated from each other by an exceptionally large 
distance of 5.2 A, so that the interpretation of To as a 
magnetic phase transition has to be discarded. Alterna- 
tive mechanisms such as a charge-ordered state as well 
as spin- d ensity- wave formation a l so have serious short - 
comings (iMatsumara et all 120091 : iNishioka eTaH l2QQ9h . 
( Tanida et all l2010al lbl) suggested the formation of Ce 
dimers within the (a, c) plane, bearing some similarities 
to a spin-Peierls transition. Neutron scat tering; experi- 
ment s confirmed the latter interpretation (Robert et all 
l2010f ). Neutron diffraction measurements gave evidence 
for a displacement of the Al atoms below To , resulting in 
a spin-Peierls transition along the one-dimensional Ce- 
Al zig-zag chains accompanied by a dimerization of the 
Ce 3+ ions as sketched in Fig. [611 

The final proof of the Ce 3+ dimerization is provided by 
the results of INS expe riments performed fo r CeRu2 Alio 
as shown in Fig. l62fa) ([Robert et a/ll2010l ). At temper- 
atures T < To there is no spectral weight at low energies. 
There is a well defined peak at 8 meV corresponding to 
the gap energy. With increasing temperature the gap 
shifts to lower energies. The peak at 8 meV can there- 
fore be identified as the singlet-triplet excitation associ- 
ated with the formation of Ce 3+ dimers below Tq which 
is supported by the oscillatory Q dependence of the in- 
tensities according to the dimer cross-section Eq. (3.9) as 
illustrated in Fig. l62f b). From the peak position and Eq. 
(3.3) the effective intra-dimer exchange coupling results 
as J e ff = —4 meV, which was con sidered to be unreal- 
istically large ( Robert et all l2010l ) . However, this value 
is based on a truncated basis with s = 1/2. The true 
exchange coupling between the Ce 3+ ions results from 
scaling with the de Gennes factor £ = (g — l) 2 j(j + 1), 
where j is the total angular momentum. For Ce 3+ with 
g = 6/7 and j = 5/2 we have £ = 0.18, and for an 
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FIG. 62 (Color online) (a) Temperature dependence of the 
energy spectra of neutrons scattered from CeRu2Alio for 
Q = 1.5 A -1 . The lines denote least-squares fits, (b) Q de- 
pendence of the intensity of the peak at 8 meV. The line cor- 
responds to th e dimer cross-section Eq. (3.9) with R — 5.2 A. 
Adapted from lRobert et qZll201cL 



s = 1/2 system (g = 2, j = 1/2) we find f = 0.75, thus 
the correction by the de Gennes factor reduces the true 
exchange coupling to Jce-Ce ~ — 1 meV. 



D. Polarons 

In complex materials, competing interactions can lead 
to the spontaneous formation of nano-sized regions of a 
different phase. If an additional charge is introduced into 
the material by doping, a fermionic quasiparticle called 
polaron can be formed. The resulting lattice polarization 
and deformation acts as a potential well that decreases 
the mobility of the charge. Polarons have spin, though 
two close-by polarons are spinless. The latter is called 
a bipolaron whose existence was the driving idea be- 
hind the discovery of high-temper ature superconductiv- 
ity in the copper-oxide perovskites (|Bednorz and Munerl 
1986). In the meantime the existence of charge-ordered 
stripes was postulated (IZaanen and Gunnarssonl Il989f ) 



and experimentally verified ([Lucarelli et a/.L 120031 ). The 



stripes are superconducting regions, separated by AFM 
regions which act as Josephson junctions by the proxim- 
ity effect. 

The existence of polarons is the key to un- 
derstand the rich phase diagrams of the giant- 
magnet qresistive man g anese and cobalt perovskit e s upo n 
doping ( Phelan et all 120061 : ISalamon and Jaimel l200lf ). 
This is demonstrated here for the hole-doped lanthanum 
cobaltates of type Lai-^Sr^CoOs. The ground state 
of the parent compound LaCoOs is nonmagnetic, cor- 
responding to a low-spin state of Co 3+ ions with Si = 0. 
It was widely believed that the addition of each hole into 
pristine LaCoOs through the substitution of a Sr 2+ ion 
for the La 3+ ion creates a Co 4+ ion in the lattice which 
has a nonzero value of Si in any spin-state configura- 
tion, thereby inducing a magnetic moment in the system. 
However, already lightly doped cobaltates with x « 0.002 
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FIG. 63 Energy spectra of neutrons scattered from 
Lao.998Sro.oo2Co03 at T = 1.5 and 10 K. The open circles 
correspond to the nonmagnetic reference compound LaCo03 . 
For clarity, the intensities of the T = 10 K data are shifted b y 
400 neutron counts. Adapted from [Podlesnva k et al. , 2008. 



give rise to an order of m agnitude larger magnetic sus- 
ceptibility than expected ([Yamaguchi et a/l Tl996). It 
was proposed that the holes introduced by Sr doping do 
not remain localized at the nearby Co site, but each hole 
is distributed among several neighboring Co sites, leaving 
the latter in the intermediate Co 3+ state (with Si = 1) 
and thereby forming a multi-site magnetic polaron. Such 
spin-state polarons behave like ferromagnetic nanoparti- 
cles with a very large total spin in an insulating nonmag- 
netic matrix. 

The existence of spin-state polarons was confirmed 
by inelastic neutron scattering experiments as shown 
in Fig. [63j A magnetic excitation is observed for 
Lao.998Sro.oo2Co03 at an energy transfer of 0.75 meV, 
which is absent for the undop ed parent compound 
LaCo0 3 (iPodlesnvak et all |2008[ ). The ground state of 
Co 3+ in the intermediate spin state is an orbitally degen- 
erate triplet which is split by a small trigonal ligand field 
into a singlet and a doublet. The transition between these 
two levels is the source of the peak observed at 0.75 meV. 
The peak intensity diminishes with increasing tempera- 
ture, in agreement with the Boltzmann population factor 
Eq. (|^2j) for a singlet-doublet transition. The Q depen- 
dence of the intensity of the observed excitation exhibits 
a clear oscillatory behavior as shown in Fig. [62J which re- 
flects the size as well as the shape of the polaron through 
the structure factor. The neutron cross-section of a clus- 
ter comprising TV magnetic ions can be approximated for 
polycrystalline material by an extension of Eq. (|3T|) : 



d 2 a 
dVtduo 



N 



cx F\Q) [<S||T W ( Si )||S') 2 +2 



i<j = l 



sm(QRij) 
QRij 



x(S\\f^(s t )\\S')(S'\\fV(s 3 )\\S)}, (101) 

For the special case of a AS = transition (which is 
relevant in the present context), the reduced matrix ele- 
ments can be factorized and set to 1. The lines in Fig. [6^1 




FIG. 64 (Color online) Q dependence of the intensity of the 
transition observed at 0.75 meV in Lao.99sSro.oo2Co03. The 
insert sketches different types of Co multimers. The lines are 
the result of s tructure factor calculatio ns based on Eq. (|10ip . 
Adapted from IPodlesnvak et a/1 [2008. 



correspond to calculated cross sections for different Co 
clusters sketched in the insert. We clearly see that the 
Q dependence of the cross section is an unambiguous 
fingerprint of the geometry of the multimers; in partic- 
ular, the data observed for the 0.75 meV transition in 
Lao.998Sro.oo2Co03 are perfectly explained by the scat- 
tering from an octahedrally shaped Co heptamer. The 
total moment of this heptamer, consisting formally of one 
central Co 4+ ion (si = 1/2) and six surrounding Co 3+ 
ions (si = 1), is 13 in good agreement with th e mag- 
netic susceptibility data ([Yamaguchi et a/lll996l ). 

The result of the above experiment gives a clear micro- 
scopic explanation why hole doping of as little as 0.2% 
dramatically affects the overall magnetic properties of 
the entire system, i.e., the magnetic susceptibility is an 
order of magnitude larger than expected. Additional 
charge carriers increase the number of such spin-state 
polarons, which form a percolative network resulting in 
a metallic state with long-range ferromag netic order at 
the c ritical Sr concentration x c = 0.18 iPhelan et all 
I2006D . The formation of spin-state polarons may be a 
common mechanism present in other Mn- and Co-based 
oxide perovskites as well. 



VII. CONCLUSIONS 

In this review we have attempted to provide an im- 
pression of the physics which one encounters in magnetic 
spin clusters, and their relevance to a variety of different 
physical systems. The presentation had obviously been 
limited to examples whose magnetic excitations have al- 
ready been very well or comparatively well understood. 
Many of the important scientific questions which are un- 
der current research or will possibly become of relevance 
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in future have however not been addressed, though some 
of them were indicated at few places in the text. We wish 
to conclude here by picking them up again. 

The small magnetic clusters have been demonstrated 
to be ideal experimental systems for studying the basic 
mechanisms of the magnetic interactions between spins 
and the underlying physical principles. The research on 
them started at least six decades ago and many funda- 
mental questions could be addressed (Sec. IIII[) . However, 
despite this long history, the research field is not yet 
exhausted, but continuously increases due to the ongo- 
ing improvements of the experimental equipment. For 
instance, for many fascinating extended magnetic com- 
pounds the corrections to the HDVV model such as the 
weaker anisotropic and/or higher-order exchange inter- 
actions have to be known with accuracy as they can be 
crucial for understanding their phase diagrams. However, 
the higher-order interactions are intrinsically hidden in 
the commonly applied analyzes of spin-wave dispersion 
relations, but become accessible by studying the cluster 
excitation spectra in the related diluted materials; an ex- 
ample was given in Sec. IIII.B.6I 

For most magnetic systems the conventional HDVV 
model, possibly including the weaker corrections to it, is 
perfectly appropriate, which is not surprising for mag- 
netic ions of type S and type Q. The basic assumptions 
underlying the HDVV picture are phenomenological, but 
there was little clear-cut experimental evidence in the 
past for an eventual failure of this model. The situation is 
quite different for clusters with magnetic type L and type 
J ions, and the exchange in general has to be described by 
more involved interaction terms. An obvious extension of 
the HDVV coupling is to introduce multipole interactions 
based on standard tensor operator techniques or to re- 
place the exchange parameter by an exchange tensor 
J(mi^mj\m r ^m f j)^ with mi and m- being the spin quan- 
tum numbers of the initial and final ionic states, respec- 
tively; the latter formalism was verified in the dimeric 
Ho 3+ compound Cs 3 Ho 2 Br 9 f Furrer et a/1 Il990h . This 
topic will certainly become more important in the future, 
e.g., magnetic molecules in corporating 4f me t al ions are 
currently intensely studied ([Klokishner et all l2009al ) . 

The observed existence of three-spin interactions in 
magnetic systems discussed in Sec. IIII.B.5I is of partic- 
ular relevance in the general context of many-body in- 
teractions, which already fascinated Kopernikus (1473- 
1543) and Kepler (1571-1630) in their investigations 
of the mutual gravitational interaction between three 
planets. There has been ample experimental evi- 
dence for many-body interactions in the past, no- 
tably in studi e s of u ltra-cold gases of alkaline atoms 
(IBiichler et all 120071), the stabi lity of molecules like 
ozone (IZhaunerchvk et all 2007), fo ur-atom exchange 
in 3 He ([McMahan and WilkinsL Il975h. c hiralitv in mag- 
netic compounds (|Grigoriev et al. , 2005), ring-exchange 



(Coldea et 



^ , l200lf) and three -body correlations of vor- 
tex states (jMenon et all l2006h in high-temperature su- 
perconductors, or domain- wall fluctuations in (anti-) fer- 



romagnets ( Shpvrko et a/l 120071 ). Novel quantum phases 
with intriguing physical properties can arise from many- 
body interactions which are usually described on the ba- 
sis of pair potentials G(r, t). However, it would be highly 
desirabl e to extend the a nalytic tools beyond van Hove's 
theory rtvan Hovel Il954f ) to include higher-order corre- 
lation functions in both space and time, for wh i ch ap - 
proximate solutions exist, e.g., for G(r , r',t) (|Hul 120071 ), 
G(r.t.t') flWord and Trammelll . Il98ll) and G(r,r',t,t') 
( van Zon and Schofieldl l200lf ). Adequate experimen- 
tal techniques should be developed as well which pro- 
vide a direct access to higher-order correlation functions. 
Indeed, novel neutron scattering techniques were sug- 
geste d for this purpose which include neutron interferom- 
etry dRauch and Sudal Il998h and spin-echo techniques 
(jGrigoriev et all\2004 ). 

Besides the small magnetic clusters the ground state 
and excitation spectrum in the large magnetic clusters, 
as they have been called in this review, has emerged as an 
attractive research field in the last 15 years. The scien- 
tific questions encountered in them are ultimately related 
to the possibility of many-body (quantum) effects, with 
links to the field of quantum spin systems. However, 
as compared to the system sizes considered in the latter 
area, in which one-, two-, and three-dimensional lattices 
of interacting quantum spins are considered, the sizes of 
the " large magnetic clusters" have to be considered small 
as they are (far) away from the finite-size scaling regime, 
in which the magnetic properties start to resemble those 
in the thermodynamic limit N — >• 00 (exceptions are spin- 
1/2 clusters because of the short correlation length for 
Si = 1/2, see Sec. IIV.CI and Fig. [28]). In this sense the 
large magnetic clusters should be called zero dimensional 
and belong to the class of mesoscopic systems. 

A typical consequence of that could be described as the 
"loss" of the wave vector as a good quantum number, 
since translational symmetry or a finite-size version of 
it such as the cyclic symmetry in wheels is generally not 
present in the large clusters. On the one hand this implies 
that concepts which are developed for extended lattices 
may have to be adapted or interpreted in novel ways if 
applied to the large clusters. The application of spin- 
wave theory to the antiferromagnetic wheels and short 
chains (Sees. HV.C.ll and IIV.D.2|) represents a showcase. 

For the antiferromagnetic Heisenberg rings or molecu- 
lar wheels the cyclic symmetry allows introducing a shift 
quantum number q which emerges into a wave vector in 
the thermodynamic limit. Therefore the available liter- 
ature results for the spin- wave dispersion relation in the 
antiferromagnetic chain, as they have been deriv e d from 
the various spin- wave theories (|lvanov and Senl [2004) , 
can be taken over directly with the wave vector replaced 
by the discrete values of q. As was demonstrated in 
Sec. IIV.C.2I for the CsFes molecular wheel this yields 
the energies in the E band with some accuracy. How- 
ever, besides this success, fundamental issues are emi- 
nent. For instance, spin- wave theories do a priori violate 
the spin rotational invariance of the HDVV Hamiltonian, 
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and their applicability to large clusters with their disor- 
dered ground states is hence fundamentally flawed, yet 
they can produce reasonable energies and one may ask 
why. Also, since the wave vector becomes discretized, 
magnetic excitations with long wave lengths do not ex- 
ist in molecular wheels (and large magnetic clusters in 
general), which raises the question of whether the excita- 
tions in the E band, albeit their energies can be derived 
by spin-wave theory, should actually be interpreted as 
spin waves. This is a sensible question, an d the notion of 
" clus ter spin waves" may be introduced ( Stuiber et all 
l201ll ). We here adopted a pragmatic view and called 
any excitation which is obviously related to a spin-wave 
energy a spin wave. 

The situation becomes even more involvled in the short 
antiferromagnetic chains. At the level of linear spin- wave 
theory the literature result for the finite antiferromag- 
netic chain can again directly be carried over, but this 
does not work for the more sophisticated spin- wave the- 
ories such as interacting spin- wave theory since the open 
boundaries result in site dependent corrections. That is, 
these theories have to be formulated in real space and 
not momentum space. Furthermore, a description of the 
excitations in classical terms, which is at the heart of the 
L&£?-band concept and motivated applying spin-wave 
theories, is actually not obvious in the short chains as the 
wave functions do not significantly overlap with the semi- 
classical configurations. Surprisingly, their ground state 
is in fact neither in the classical nor the quantum regime 
and the intuitively clear distinction between these t wo 
regimes becomes blurred ( Konstantinidis et a/ll201lh . 

The loss of the wave vector on the other hand is also 
indicative of the possibility of lattice topologies in large 
clusters which are not possible in extended systems. In 
fact, most magnetic molecules which were synthesized 
have a complex topology which cannot systematically 
be expanded into an infinite lattice. The V15 or Mni2 
molecules presented in Sees. IIV.E.1I and IV. Bl are ex- 
amples. However, for obvious reasons, the research on 
the many-body aspects of magnetic cluster excitations 
concentrated on lattice topologies or magnetic molecules 
with an "appealing" symmetry, but the most interesting 
magnetic phenomena may be overlooked this way. In this 
sense only the simplest systems were studied so far, yet 
the understanding of their excitations can be challenging 
and the case of Fe3o establishes a dramatic example. The 
large number of lattice topologies available through the 
molecular nanomagnets, which thanks to the productiv- 
ity of chemists will certainly grow further, presents obvi- 
ously a wide area for future research. 

Except of t wo c ases ([Carretta et all 120071 ; 
Waldm ann et al. I, l2QQ6al ) (see also Sec. HV.E.ip the 
INS works on molecular nanomagnets were using pow- 
der/polycrystalline samples. Experiments on single 
crystals would offer new possibilities to unravel the 
many-body physics in the cluster excitations, e.g., 
the spin-pair correlation function can be mapped out 
more directly than possible with powder samples. 



Single-crystal INS experiments have become state of 
the art in the last decade for inorganic compounds. 
The current drastic improvements in the INS technique 
make such experiments possible now also for molecular 
nanomagnets (which are challenging because of the small 
number of magnetic centers as compared to the many 
non-magnetic ligand atoms). Single-crystal INS work on 
molecular nanomagnets will certainly been seen more 
often in the next years. 

Magnetic cluster systems have also been considered 
as promising units for quantum computing. Much ef- 
fort went into "mesoscopic spin- 1/2 clusters" or "anti- 
ferromagnetic cluster qubits", i.e., large clusters which 
exhibit a total spin S = 1/2 ground state and may 
be used as qubits at low temperatures, where only 
the gr ound state is thermally populated ([Meier et all 
2QQ3al lbl). Such clusters may present advantages over 
atomic scale qubits, such as easier addressing and read- 
out. Significant progress were made, e.g., coherence 
times long enough for quantum computations could be 
achieved in the Cr^Ni "doped" antiferromagnetic wheel 
([Ardavan et all l2QQ7h and two molecular " qubits" could 
be linked such as to provide entanglement between 
them suggesting the possib i lity of two-qubit op erations 
([Affronte et a/.U2QQ5l l2006t ICandini et aZl l2Ql(jh . How- 
ever, the multi- level structure of the low- lying energies 
as it can be provided only by magnetic clusters has also 
been explored. The particular structure of the ground- 
state spin multiplet in the single-molecule magnets was 
theoretically shown to allow for an implementation of 
Grover's search al gorithm or to build dense an d efficient 
memory devices ([Leuenberger and Lossl l2QQlh . Further- 
more, by accessing also higher- lying levels nearby to 
the ground state Rabi oscillations could be observed in 
the V 15 molecule and the Fe4 single- mol ecule magnet 



([Bertaina et all 120081 : ISchlegel et aR 120081 ). Lastly, the 
spin frustration or orbital degeneracy in the ground state 
of regular spin triangles has been theoretically shown to 
allow for a coupling of the spin degree of freedom to 
electric fields and currents (spin-electric currents), such 
that the spin qubit can be manipulated through cur- 
rents supplied to a cluster by e.g. scanning tunneling 
microscope (STM) techniques ( Georgeot and MilaL 12010b 
iLehmann etaH 120071 : iTrif et a/.U 2008l). These~examples 
indicate that the cluster excitations as present in mag- 
netic clusters may provide novel quantum computation 
schemes, and exciting results in this direction can be ex- 
pected in the future. For instance, the many-body nature 
of the large magnetic clusters suggests the possibility of 
decoherence free subspaces in them, which however has 
to our knowledge not yet been explored. 

Magnetic clusters play also a role in the emerging field 
of molecular spintronics (IBogani and Wernsdorferl 120081 : 
iRochal l2005t ISanvito and Roc ha. 2006). Significant suc- 
cess has first been obtained with molecules containing 
one metal ion, which have been incorporated into break 
junctions or contact leads produced by electro migration. 
For instance, the Kondo effect could be observed in the 
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mole cule [Co(terpy(CH 2 )5-SH) 2 ] 2+ this way (|Park et all 
2002). However, also the current- volt age characteristics 
of magnetic clusters of two metal centers or even the 
single- mol ecule magnet Mn i2 have been measured in ex- 
periment (|HeerscheL [2006). Man y further schemes are 
curre ntly under investigation ( Bogani and Wernsdorferl 
120081 ). As in the context of quantum computation it 
can also be envisioned for molecular spintronics that the 
unique energy level structures or complex many-body 
states provided by magnetic cluster systems will open 
novel opportunities and allow for functionalities not cur- 
rently considered. 

The molecular nanomagnets represent obviously a 
valuable resource of (large) magnetic clusters, but in 
the recent years not only the chemical route towards 
large clusters has been advanced, but several beauti- 
ful exam ples of artificially eng ineered spin clusters also 
emerged (Jamneal a et al .,2001). Here ensembles of mag- 
netic metal ions experiencing nearest-neighboring ex- 
change interactions were fabricated directly on surfaces 
using STM techniques, and their magnetic excitations 
measured through recording the current- volt age curves. 
For instance, short antiferromagnetic chains of Mn 2+ 
metal ions with chain lengths of N = 1 — 10 were pro- 
duced and the magnetic ground stat e as well as first ex- 
cited excitation observed this way (|Hirjibehedin et a/1 
2006). Apparently in both areas, that of the molecular 
nanomagnets and the artificially engineered spin clus- 
ters, one faces very similar scientific questions, yet the 
technological challenges in putting them forward in real- 
world applications are complementary. In comparison 
with the artificially formed clusters, the molecular nano- 
magnets are available in an abundance of different lattice 
topologies assuming complex many-body states and their 
excitations can be studied in any detail using powerful 
experimental techniques such as inelastic neutron scat- 
tering. However, the artificial nanostructures are avail- 
able directly on surfaces and have already been proven 
to maintain their magnetic properties and their function 
to be addressable on the surface. The molecular nano- 
magnets and the artificially engineered spin systems are 
hence complementary in the sense that the advantages of 
each class may help the ot her to overcome its problems 
(Konsta ntinidis et al. 1. 120111 ). 

Furthermore, magnetic cluster excitations are funda- 
mental in a variety of other systems. For instance, the 
clusters may be linked together by a network of weak 
magnetic m£er-cluster exchange interactions, such that 
the cluster excitations may "travel" through the net- 
work and become dispersive. This can lead to fascinating 
novel quantum phases, and the Bose-Einstein conden- 
sation in the dimer-based compounds, which was pre- 
sented in Sec. lVI.B| is a striking example. However, such 
networks may obviously not only be built from dimers, 
but also trimers or tetramers, which can additionally in- 
troduce spin frustration and possess unprecedented be- 
havior. At this point a further link between molecular 
nanomagnets and quantum spin systems emerges, as not 



only small clusters may be incorporated into the network 
but also large clusters or molecular nanomagnets, which 
should be expected to result in an interesting interplay 
between the complex quantum states realized within a 
cluster and the complex phases generated in extended lat- 
tices. Indeed, the synthesis of extended exchange-coupled 
networks of, e.g., CU3 molecular unit s or single-molecule 
magnets have a l ready been reported pvanov et all , l2010l : 
iMivasaka eTaH 120061 : iMorimoto etaR l2009h .~All in all, 
the competition between the quantum states in a mag- 
netic cluster and the cooperativity introduced through 
an extended network of z?7ier- cluster interactions should 
continue to be an attractive playground for research. 

Lastly, clusters of magnetic ions are relevant in bi- 
ology as they constitute the active sites in many met- 
alloproteins. They are ubiquitous in living matter and 
contain sites with one to eight metal atoms, sometimes 
with multiple occurrence of the smaller clusters in the 
same protein molecule. The primary function of, e.g., 
the iron-sulfur clusters lies in mediating one-electron re- 
dox processes and as such they are integral components of 
respiratory and photosynthetic electron transfer chains. 
So far the understanding of the electronic ground- and 
excited-state properties relied on magnetic susceptibil- 
ity, magnetic circular dichroism, EPR, Mossbauer and 
resonance Raman scattering measurements, which pro- 
vided information on the anisotropy of the g factor, the 
valence state of the iron atoms, the total spin quantum 
number of the ground state, and sometimes the intra- 
cluster exchange interactions (Bei nert et aZ.L [l997h . The 
latter are often determined from experiments on model 
systems as demonstrated, e.g., by magnetic susceptibil- 
ity measurements of a cubane-type Fe4S4 cluster, giving 
J = —18 meV for the exchange para meter of t he tw o 
semi-independent (Fe2S2) 2+ dimers (|Yoo et all [l997). 
To our knowledge, the very informative INS technique 
has not yet been applied to the study of metallopro- 
teins, since prohibitively large deuterated samples with 
volumes of typically 1 cm 3 would be needed. However, 
with the advent of third-generation neutron sources and 
novel beam- focusing techniques, reducing the requested 
sample size by several orders of magnitude, such experi- 
ments will become feasible in the near future. 

As so often, improvements in the experimental tech- 
niques are a main driving force for scientific progress, and 
it is finally mentioned that the development of the instru- 
mentation used for studying magnetic cluster excitations 
has been impressive in the last decade, and should be 
expected to be so also in the future. Neutron scattering 
can serve as a representative example. Neutron facilities 
such as the high-flux reactor HFR at the Institute Laue- 
Langevin (ILL) in Grenoble (France) and the spallation 
neutron source ISIS at the Rutherford Appleton Labora- 
tory in Didcot (U.K.) which both presently undergo sig- 
nificant upgrades as well as the successful commission- 
ing or planing of new third-generation neutron sources 
such as the Spallation Neutron Source (SNS) at the Oak 
Ridge National Laboratory (USA), the spallation neu- 
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tron source at J-PARC in Tokai (Japan), and the Eu- 
ropean Spallation Source (ESS) in Lund (Sweden) open 
truly exciting perspectives to unravel the many unsolved 
questions in the area of magnetic cluster excitations. 
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